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Partially Plastic Thick-Walled 
Cylinder Theory 


By M. C. STEELE,' URBANA, ILL. 


thick-walled 
A quanti- 


Previous theories for partially plastic 
cylinders under internal pressure are reviewed. 
tative comparison is given for (a) compressibility versus 
incompressibility of material, and (6) von Mises’ versus 
Tresca’s theory of failure. The former reveals that deflec- 
tions at the outside and bore surfaces agree closely, al- 
though considerable percentage differences may be found 
Large differences (except 
two 


in the axial stresses and strains. 
for the axial stress) are found in comparing the 


theories of failure. Based on the comparison and available 


experimental evidence, a theory is presented in closed 
form to include the Hencky stress-strain relations, in- 
compressibility, and Ludwik’s strain-hardening function. 


NOMENCLATURE 


The following nomenclature is used in the paper: 


do; de = principal stress increment; principal strain in- 
crement 
€ = principal stress; principal strain 
= principal stresses and strains in hoop, radial, and 
axial directions, respectively 
= elastic component of strain; plastic component 
of strain 
= effective stress; effective strain 
= tensile yield stress 
= yield shear stress in tension 
= maximum shear stress; maximum shear strain 
= maximum shear strain at yield 
= internal radius 
= variable radius 
= plastic-elastic boundary radius 
= outside radius 
, aX, B = plastic flow functions 
P = internal pressure 
G; } = modulus of rigidity; elasticity 
m = strain-hardening parameter 
v = Poisson's ratio 
sur- 


(€@),/4, (€@)bore = circumferential strains at outside and bore 


faces, respectively 


INTRODUCTION 


The partial overstraining of thick-walled cylinders by internal 


1 Research Assistant Professor, Department of Theoretical and 
Applied Mechanics, University of [Illinois 

Contributed by the Applied Mechanics Division and presented 
at the Annual Meeting, Atlantic City, N. J.. November 25-30, 1951, 
of Tae American Soctrety or Mecuanicat ENGINEERS 

Discussion of this paper should be addressed to the Secretary 
ASME, 29 West 39th Street, New York, N. Y., and will be accepted 
until July 10, 1952, for publication at a later date Discussion re- 
ceived after the closing date will be returned 

Note Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those 
of the Society Manuscript received by ASME Applied Mechanics 
Division, June J1, 1951. Paper No. 51—A-25 


uniform pressure has claimed the attention of several prominent 
Table 1 lists some of the contributions and illus- 
The 


matical solutions require numerical procedures for their comple- 


investigators 


trates differences in basic assumptions, rigorous mathe- 


tion. In general, the solutions are so complex that little practical 
use may be made of them by design engineers. 

The mathematician’s approach to the problem under con- 
sideration is fundamentally different from that of the engineer. 
The former is interested in an accurate knowledge of the stresses 
and strains in terms of the applied load or loads for all points in 
an idealized body. In general, the engineer seeks an estimation 
of the variation of factors with applied loads, which are of signifi- 
cance only in the functional design of the part. The rigorous 
solution due to Hill, Lee, ? shows that in Nadai’s 
(2) approximate solution, the axial stresses are in error by as much 


and Tupper (1 


as 667/; per cent. This factor may or may not be of interest to 


the engineer. For materials exhibiting flat-topped (well-defined 
yield 
hardening, the engineer 
In the case of the thick-walled cylinder, interest lies in 


stress-strain curves, or small amounts of strain- 


is interested in deformations rather than 


point 


stresses 
deflections of bore and outside surfaces for specific internal pres- 
sures. An accurate knowledge of the stresses would be required 
only if rupture theories were to be considered before the deflec- 
tions become the governing feature in design. The axsal stress 
is small, in general, compared to the other two principal stresses, 
and a considerable percentage error in the axial stresses may be 
tolerated if their effect on the bore and outside deflections is 
small. 

The necessity for accurate solutions to the thick-walled cylin- 
der is not denied. Indeed, history has shown that present-day 
elastic mechanics of materials could not have been so soundly 
formulated without the solutions of the theory of elasticity on 
which to base it. It is believed that a similar evolution must 
take place in the case of plasticity if the advantages it has to offer 
are to be understood and used widely by industry. The thick- 
walled cylinder under internal pressure has been widely treated 
by the methods of the mathematical theory of plasticity, and 
the author feels particularly fortunate in having these solutions at 
his disposal for an attempt at formulating the problem on a 
“plastic mechanics of materials’’ basis. 

It is thus the purpose of this paper to examine the causes lead- 
ing to complexities in the analyses listed in Table 1, and by a 
partial elimination of them, to present a theory which, it is hoped, 
will be of use to design engineers 

JUTLINE OF SOLUTIONS 

Table 1 will be better understood by outlining the require- 
ments for a solution for the stresses and strains in a partially 
subjected to internal 


plastic thick-walled cylinder pressure 


They are as follows: 
(a) Radial Equilibrium 


Aye, 


Oy 


{1} 


2 Numbers in parentheses refer to the Bibliography at the end of 
the paper 
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TABLE 1 
Plastic 
stress- 
strain 

Author law 
Prandtl- 


PARTIALLY 


PLASTIC THICK-WALLED-CYLINDER THEORIES 


Compressibility 
Plastic Elastic 
domain domain 


Flow 
condition 


I'ype of 


Axia! boundary solution 


Comp Comp Tresca a=0 Numerical 


JUNI 


1952 


Hill, Lee, and Tupper 
( 


Comp Comp 


) 
Hill, Lee, and Tupper 


1 
MacGregor, Coffin, Hencky Comp Comp 
and Fisher (12) 
Sokolovsky (16) 
Sokolovsky (16) 
Sokolovsky (16) 


Hencky Comp Comp 
Hencky omp Comp 
Hencky comp Somp 


Hencky 
Prandtl- 
Reuss 


Hencky 


Sokolovsky 

Hodge and 
(10) 

Hodge and 
(0) 


nxcomp Comp 
‘omp Comp 


16) 
White 
White 


comp Comp 


Allen and Sopwith Hencky ‘omp Comp 
11) 


Hencky Incomp Comp 


Cook (17) 


Nadai (2) Hencky Incomp Comp 


Steele (13 Hencky Incomp Comp 


Steele (theory of this Hencky Incomp 


paper) 


* Inclusion of strain-hardening 
defining theory of failure 


Comp or 
incomp 


Tresca Closed ends Num? 


von Mises Open ends 
von Mises 
von Mises® 
von Mises 


von Mises® 
von Mises 


von Mises 


Kither « 0, or Closed 
closed ends, or open form 
ends 

Closed ends 


Tresca 


Closed 
form? 


lresca with 
upper and 
lower yield 
von Mises Closed 
form 
Closed 
forn 


€: 0 


Either « 0, or 
closed ends, or open 
ends 

Either « = 0, 
closed ends, or open 
ends 


Tresca 


Trescae or Closed 


foru 


nonasterisked solutions assume material to flow according to functic 


Solution dependent on experimental measurements which Cook made on closed-ended cylinders 


Straure Con patib lity 


hither 


"ke 
/ o,y dy = 0 (open ends 
/ 1 


~) p 


o,y dy = — (closed ends 


It should be noted that a comparison of the theories in Table 

is hampered by inconsistency among authors in the selected 
axial boundary. Choosing the axial strain as zero is convenient 
mathematically and is popular for the more rigorous solutions. 
It is a condition, however, not found in practice where cylinders 
have either open or closed ends. Solutions for ¢, = 0 and closed 
ends do not differ greatly, but the difference is larger for open- 


ended cylinders 
(d) Plane Sections to Remain Plane 
€, is independent of 7 


Elastic-Plastic 
implies continuous €g and € 


e) Continuity Across Interface. Continuous 


Equilibrium de- 
Conditions at the elastic-plastic boundary 


displacement ™ 
mands continuous @,. 
if incompressibility of 


may be deduced from these. For example, 


plastic material only is assumed, then ¢€, cannot be continuous 
Also, if the tensile stress-strain relation is discontinuous (e.g.. 
upper and lower yield points), then og must be discontinuous when 


Tresca’s yield condition is used. 


(f) Plastic Stress-Strain Relations. One of two plastic stress- 
strain relations replaces the well-known Hooke’s law for elastic 
The Prandtl-Reuss (3) relations provide the basis for 


They take 


bodies. 
the 
the form 


“incremental strain” or ‘‘flow’’ type of theory 


Ede, = do; vd (a: + + did 
(plus two similar equations for the other two principal st 
increments ). 

The Hencky (4 


formation” theory. 


relations describe the “‘total strain’ or 


They may be writter 


(plus two similar equations for the other two prin ipal strains 


The basic differences in Equations [5] and [6] 
It is well established that the 


have been de- 
scribed adequately elsew here ( 4) 
incremental theory is the more correct, and that if the principal 
axes of stress remain fixed in direction, and the principal stress 
ratios are constant, the Prandtl-Reuss and Hencky relations 
vield the same result. Frequently, even though the conditions 
on stress axes and ratios are not strictly maintained, solutions 
using Equation [6] vary only a small amount from those solved 
by the Prandtl-Reuss relations. This is the case for partially 
plastic thick-walled cylinders subjected to internal uniform 
pressure only, and the small inaccuracy in using the Hencky re- 
lations is compensated for by certain advantages of mathe- 
matical convenience. 

Vaterial The total strain in 


an overstrained material is composed of an elastic and a plastic 


(g) Compressibility of Plastic 


component; thus 


It is observed from experiments on metals that there is no change 


in volume due to the plastic component of strain, Thus volume 


changes during plastic deformation are elastic; hence 


1 — 2 
&4+@2+6 = (a, + Ge + G:) 


E 





STEELE 


If the stresses and strains obey Equation [8], the material is said 
to be compressible. If the plastic component of strain is pre- 
dominant 


&+6@7+6 20 9) 


When Equation [9] is used, the material is said to be incom- 


pressible. In the cylinder problem under consideration the plastic 


material is restrained by an outer elastic hoop The plastic 
strains are thus comparable in magnitude to the elastic ones, and 
errors are introduced if Equation [9] is used. The inclusion of 
compressibility, however, considerably complicates a solution. 
It is seen from Table 1 that there is only one solution in closed 
form observing compressibility of the plastic material. 

(h) Flou The 


(maximum shear stress) and von Mises (maximum energy of 


Condition theories of failure due to Tresca 


distortion ) are well known They are usually stated as 


a [10] 


where 


and 


a; o:)? +(e 


{11} 


All the cylinder theories incorporating strain-hardening make use 
The “effective 


stress’’ 


of a generalization of Equation [11 


and “‘effective strain” are defined as 


for all combinations of stress 
As 
Hencky stress-strain relations, however, this 


and it is assumed that ¢ = f (é 
Hill (5) has pointed out the limitations of this assumption. 
for the case of the 
strain-hardening function has approximate application in many 
instances and in particular to the problem under consideration. It 
is not common to make use of an extension of Equation [10] for 
Ludwik (6) in 1909 stated that the 
maximum shearing stress as a function of the maximum shearing 


strain-hardening material. 


strain described the subsequent flow of a ductile material in the 


strain-hardening range. Little use has been made of this 


as- 
sumption. Tests on low-carbon steel thin cylinders* under com- 
bined stress have shown that von Mises’ theory closely describes 
failure and subsequent flow condition. In steel with a drop in 
stress at the yield point, Morrison (7) states that initiation of 
vield is due to a critical maximum shear stress but that plastic 


flow continues in accordance with von Mises’ criterion 
Experimental work by Steele and Young (8) on low-carbon 
steel, and Macrae (9 


under internal fluid pressure, contradict the findings of the pre- 


on nickel-gun-steel thick-walled cylinders 
vious paragraph from thin-walled cylinder tests. Strains were 
found to agree closely with theories based on the Tresca condition 
and were in definite disagreement with von Mises’ criterion. It 
is of importance that experimental results lay on the side of the 
strain curves computed using Tresca’s condition, remote from the 
curves based on von Mises’ theory. Steele and Young attributed 
this to the mechanism of flow in steels containing a flat-topped 
They showed that the theo- 


retical assumption that a partially plastic thick cylinder consists 


portion at yield in a tension test 


of two regions, both containing isotropic material and separated 
by a circular boundary, is incorrect. By polishing the ends of 
cylinders it was shown that the overstrained portion consists 
mostly of elastic material separated by a small number of thin, 


spiral, wedge regions of plastic material 


3 Bibliography (2), chapter 17 


PARTIALLY PLASTIC THICK-WALLED CYLINDER THEORY 


A comparison will be effected in the next section of this paper 
among certain of the theories listed in Table 1 Previous com- 
parisons, with one exception (10), have involved several varia- 
bles, e.g., axial boundary, compressibility, theory of failure, 
and so forth, and the effects of the individual variable have been 
masked. It is proposed to adopt a standard set of conditions and 
examine the effeet of compressibility and theory of failure on the 
significant stresses and strains, In addition, the effect of making 
Poisson's ratio 0.5 rather than the normal value (0.3 for steel) is 
discussed for theories which assume in ompressibility of the 
plastic material 

Comparison oF THEoRIES* 


Hodge and White (10) have 
of the Reuss stress-strain relations and the Hencky relations for 


given a quantitative comparison 
cylinders in plane strain, and composed of compressible material 
obeying von Misés’ flow condition. Their results show that the 
two theories yield almost identical values for the stresses and 
strains. It is realized that this cannot be generalized strictly 
to the cases of open and closed cylinders No solution is known 
to the author which includes the incremental theory, compressi- 
bility, and von Mises’ flow condition for closed or open-ended 
cylinders. Incremental solutions using Tresca’s theory all yield 
the same values for hoop and radial stresses since the equilibrium 
equations may be integrated directly (for zero strain-hardening); 
the solution for these stresses is statically determinate. It is 
difficult to estimate the errors involved in the computed strains 
by using the Hencky stress-strain relations for open-ended cylin- 
ders rather than the Prandtl-Reuss equations, but there seems at 
present no alternative to the selection of the former relations for 
the basis of a practic al solution 
A quantitative comparison follows to examine the effect of (a 

compressibility, and (5) theory of failure. The basis of the com- 


parison is as follows 


1 Hencky (total strain) plastic stress-strain relations 
2 Material exhibits idealized, flat-topped tensile stress-strain 
curve 
Open-ended cylinders 


=2 


3 
4 Diameter ratio (k 


5 Stresses plotted for cylinders vielded halfway through the 
wall(m = 1.5). 

The comparison is effected with reference to three theories, 
viz., those of (a) Allen and Sopwith (11), (6) MaeGregor, Coffin, 
Fisher (12), (13) 


theories may be obtained from Table 1 


these 
Allen 


and Sopwith’s, and Steele’s work are compared, the effect of com- 


Information 
It is seen that if 


and and (c) Steele on 


pressibility may be estimated with all other variables standard- 
ized, Similarly, 
MacGregor, Coffin 
stresses and strains for the von Mises and Tresca flow condition 


if Allen and Sopwith’'s theory is compared with 
and Fisher's solution, the differences in the 


are found 

Figs 1, 2, and 3 show plots for these three theories of the prin- 
( ipal stresses and strains at the outside and bore surfaces 

With regard to compressibility versus incompressible (plastic 
it is seen that the only differ- 


In Steele’s solution a disconti- 


material) solutions for the stresses 
ence lies in the axial direction 
nuity arises at the plastic-elastic boundary when Poisson's ratio 


for elastic material is taken as 0.3, since 
l 2v 


on the elastic side, and 


€e + 


is included in a PhD thesis sub- 
of Glasgow 


* The comparison presented here 
mitted by the author 
April, 1951 (see Bibliography 13 


to the University Scotland, in 
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RADIAI MFERENTIAI STRESSES FOR 


tes oF (a) ALLEN AND Sopwits (v = : (b) MacGreGor 
= ns Fi } Bore CrirncUMFERE) 
IN, AND FisHer (+ 0.3 c) STEr 0.3 ann 0.5 
ALLEN AND SoPWITH 


PIsHEeR (¥ 


on the plastic sice his would be con 


urate knowledge of the axial stresses is re« 
pointed out in the introduction, tl 
bore surfaces of the cylinder are 

that 


neer, and Figs. 2 and 3 show 
ire small. Indeed the plots 


or Allen and Sopwith’s, and Steek 


tiated on the se ile used The largest 


noticed in the axial strain which has 
ferential strain. 

To eliminate the physically unrealist 
stress and bore circumferential strain, it has been suggest« 
Nadai® and Prager (14) that the elastic material should be 
sidered incompressible also, Le., V The stresses and 
strains for vy = 0.5 are shown plotted in Figs 2, and 3, and it is 
seen that the curves show closer agreement with the more correct 
theory of Allen and Sopwith 
axial stress and strain, but the outside 
strains vary at the most by 3 per cent 
once the plastic region has ext« 
wall thickness or greater Thus it is co 
based on incompressibility in both the pl 
will be of value to the engineer, su is in ted in factors 
which are not affected seriously b dic issumptions 
It should be emphasized that the f going Is m ‘ ral state- 
ment applicable to all problems in pla it t one suitable for 
problem 


a plastic mechanics-of-mater! ils s 


Theories of failure are compared | ylutior len and 
Sopwith, and MaeGregor, 


lassed as rigorously corre 
Fig. 2 rsipe CIRCUMFERENTIAL (TENSILE) AND AxtaL (Com assed as rigorou 

PRESSIVE STRAINS FOR THEORIES OF a) ALLEN AND SOPWITH 
¥ 0.3 b) MacGrecor, Corrin, ano Fisner (¥ = 0.3); ( 


STEELE 0.3 AND US 
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large diff ‘ mput tresses and strains® (with the 


exceptior ind decision as to which theory to 


use depends on the 


xperimental behavior of the cylinder mate 


rial 


STRESS 


It was pointed out in the introduction that tests performed 
on thick-walled cylinders by Steele and Young, and Macrae, gave 
results at variance ith von Mises’ theory 
ted to the mecha 


on Mises’ theory may be 


und this was attribu- 
ding. It is further reasoned that 


SHEAR 


operative within the wedge regions 





of overstrain, but t tegrated effect of such regions on observed 
strains at the outsid d bore surfaces is more closely in agree SHEAR STRAIN 
with Tres« \ Evidently some new theor must bie. 4) Linear Srrarn-Harpenine 
” devised t nate | ver-all effects of wedge regions ot 


nent 


DIaGRaM 
*verstrain in an elastic mats il, but so far this has proved ex- 
ceedingly difficult 

Most gun steels have short flat-topped portion to the tensil 
stress-strain curve before strain-hardening becomes operative 
It is suggested that the flat-topped portion gives rise to the wedg 
type of yielding in members in which 
Further, for 


a stress gradient 

geometri i Hnlial 

materials and exhi 

tion in «a tension tes amount of strain-hardening withi 

wedges will be dependent on the length of the flat portion and 

the strain-hardening characteristics of the material. Experi 
have shown (Steele and Young) that for mild steel, whicl 

has a long flat portion, compared to the strain at criti 

linder outside and b circumferential strains may 


eri 


ribed approximately by Tresca’s theory (zer 


IL is Suggest 1 i 1 


o strain-harde 
pPiastic mec hanics-of-ma 


esca’s theory as extended t 
xisting in a cylinder of strain-harden- 
y in closed form is given in the next 


simplifying assumption of Incompr 


material 


20] 


20! gives 
Dis a 
equation 


satisfies the 


Mises-Tresca « 


the assumptior laws agree for 
follows the cust g treatr t for theories 
However, the two flow iditior 
equally well so stres 
for example, t 

torsion) were 
theories for tl P i \ r r ind muous I boundary 
consideration probabl ‘ 1b 7 continuou tic bound 
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Hence in elastic domain (n < where 


From Equation [13] 


‘ L oe 
and in plastic domain (1 < y 


in the plastic domain. Substituting for €@, 8, and o@ 
] ’ log ( ) ’ 


ms () 


P, hence 


m) log n? 25) 
The strain at the outside diameter is obtained from the elastic 
equations 
juation in the elastic region from the elastic « quations 
Kee = 4 The end condition gives 
substituting for og and a, 


~} p 
o,y dy 0 (open en ° closed ends) 
(Eee Qn j ' 9 


Assuming open ends 


*nJ3 E m 43 + 
n 
jdy = 8 
J / 2(1+yp 


From which 


From Equation | Summarizing the domain 
€@ € 


Ce Oo, 


Substituting for y and 7 from Equations [18] and [19], respectivel 


n 


2G (1 
Eee 


In the plastic domain 


2G 


In particular, the bore circumferential strain is given by 


E€@)vo 1 Ee 
=(1 +» ; 


s 9 


28 
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where 





Lc 


— TONS 


hog n* 


oe 6; 


Curves for ¢,/s, o/s and @,/s are given in Fig. 5 for various 
values of the strain-hardening parameter m 
Figs. 6 and 7 show plots of 
Fes E€@ wore 


and 
8 


respectively, for various values of m and for 1 
CONCLUSIONS 


Previous partially plastic thick-walled-cylinder theories are 
highly mathematical and not suitable for use in design offices A 
search for assumptions leading to a simplified theory for design 
may be made from the point of view that engineers are interested 


in certain significant factors—for the most part in deflections of 
it Pic. 5 Raptat, CircUMPerenTIAL, AND AXIAL STRESSES FOR 


outside and bore surfaces for specific applied pressures is 
Various Values or Srrain-HargpentInc PxramMerer mn 


shown that on the assumption of incompressibility of both elastic 
ind plastic materials, a solution for the significant factors in de- 


OPCN 


K-2 


MFT RENTAL STRAINS FOR varmioUs 
STRAIN HAROCNING PARAMETER, m 











Cte)a,* 





a 
E fe )qoer 
6 


hie. 6 OvutTsipe CIRCUMFERENTIAL Tensite) anp Axiat (Com- 
PRESSIVE) STRAINS ror Various VaLvues or Strrarn-HarpENING Fic. 7 Bore Circumrerentiat Strains ror VaRtous VALUES OF 
PARAMETER m STRAIN-HARDENING PARAMETER m 
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sign is in close agreement with a more rigorous mathematical 


theory based 


mild-steel thick-walled 


upon compressibility. Experimental work on 


cylinders under internal fluid pressure 


disagrees with von Mises’ (maximum-energy-of-distortion) theory 
of failure and more closely 


supports that due to Tresca (maxi- 


shear stress On the basis of the foregoing conclusions 


Hencky stress-strain relations 


mum 


an extension ol i’s theory to cover strain-hardening, incom- 


pressibility, and is considered a 


sound basis on which to furnish a “‘plastic’’ mechanics-of-ma- 


terials solution. Such a theory is presented in closed form, i. 


f stresses and strains applied pressures 


the values o 


without resort to Jen , rical integratior 
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On the Stresses in a Notched Strip 


B 


In a previous paper by the author (1),* a theoretical solu- 
tion for a notched strip under longitudinal tension is 


The 


system of linear equations. 


infinite 
ot 


It seems, how- 


given. result demands the solution of an 
A considerable amount 
labor is involved in solving such a system. 
ever, that the labor can be diminished by adapting to the 


In 


this paper such a process is described and then applied to 


solution a process known as the promotion of rank. 


solve the problem of a notched strip under transverse bend- 
ing. The solution of this problem seems also to be new. 
The numerical results obtained are compared graphicalls 


with the experimental results available. 


tension along 
by the 


tion in qué 


present 


suthor tion is 


tructed b ven biharmonk 


which represents 
tl widtl 


ions, in a 


en stress in unnotched strip of 


the 
of the 


1 Same 


ch of two 1armonic functions gives no stress at 


angential 


F the 


rs of 


infinity str uidition, gives no normal or t 


i} 


stresses along bot pt at the cents 


remaining boundary 
satisfied 


Thi 


The stresses 


notches which con 
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rank consists in replacing them by two The coefficients pox, pes are determined successively by 


process ol promotion ol 
and we’, ws of ranks equal the set of k-pairs of equations 


sets of new tunctions We, Ws, 
that is to say, such that rr, derived from them } 
. P2t,oPo,0 Pr 
20, cos 46 and r@, begins with terms Doro 
i zoo, 
For the sake of completeness, 


to their suffixes, 
begins with terms in cos 
in sin 48, sin 68, , respectively 

we shall write wo, wo’ in place of Xo, x We form we. we’ by the , a 
Por ol 


rule 


where p q ire coefficients to be adjusted so that rr, de- 
rived from ws, @:' may contain no term independent of @ and at 


the same time r#, may contain no term in sin 26. These condi- 
tions lead to two pairs of relations 
Similarly, gee on, ¢ ire determined by 


so that they determine p Poo’ and gro, 4 simultaneously 


6, derived from w: as 
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The coefficients 
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Coon, ¢ oa respectively In general, we 
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In the previous paper function 


structed in the form 


Now, by means of the two se 
structed in the following form inste 


Also, we write the stre 


For similar reasons, the factor 6?7' is introduced to render the 


less. Besides 


parametric coefficients ¢ ind D., dimensior 
and r@, at the edge of the lower notch, derived from the basic 


sin 280 13 stress function Yoo, will be denoted, in gener al, between the r ange 


0 w/2 and #2. by 
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dy, sin 280 20 
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Note that xo is here tacitly assumed to be even in xz. The 
boundary condition of no traction at the edge of the lower notch 


requires that rr; = r@, = 0 when derived from x. It leads to 


bo + Cobo.o + 
be + Cobo.» + 


Doboo’ = 0 


Dobo2' + Cobo. + Dobra 


Dexbar2e') = 0 


+ Dords, 


The parametric coefficients ( and Des can now be determined 
successively pair by pair 

It is noted that the foregoing process, as it stands, is applicabk 
to cases in which the strip is notched symmetrically by two semi- 
circular notches of equal radii under the action of a basic stress 


For, 


condition of no traction at the edge of 


function which is even in z but may be either even or odd in y 


in such cases, the boundar 
satisfied automatically if the condi- 


the upper notch will be same 


tion at the edge of the lower notch is satisfied 


It is noted further that by this process additional parametric 
coefficients, if required, can be computed without affecting pre- 


vious computations. This property is sometimes of distinct ad- 


vantage in numerical 
When the 


a straightforward matter to calculate the stress at any point in 


omputations 


parametri oefficients are thus evaluated it is then 


the strip 


TRANSVERSE-BENDING PROBLEM 
The foregoing process of promotion of rank will be applied to 


solve the transverse-bending problem for the notched strip as 


shown in Fig. 1. The basic stress function of an unnotched strip 
under a constant transverse-bending moment 6*7' in the plane of 


the strip is given by 
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that this stress system is even in z but odd in y. The two 
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Table 1 reproduces the values of 
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TABLE 2 (continued 
2k = 6 

Coefficient box.2. 
07960 
76868 4.95592 
346409 18. 60187 
62763 0.96173 762 
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rhe other limiting value 3.07 at A = 0 agrees with the corre- 
I ig 2 shows the results 


sponding value in the tension problem 
graphically, where the experimental values available are also 


plotted for comparison (6, 7 
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An Alternative Formulation of the 
Laws of Mechanics 


By H. M. TRENT, 

The author proposes a set of laws of mechanics which 
are duals of those of Newton. These laws are based upon 
the concept of the continuity of transmitted forces around 
closed loops. They are particularly useful in those engi- 
neering problems where only knowledge of forces is of in- 
terest, since the laws lead directly to equations in which 
forces are the dependent variables. 


years have elapsed SINnCt 
Dur 


§IWO hundred and -three 
Newton enunciated his celebrated laws of motion. 


sixty 


ing this time, many ol the world’s greatest scientists have 
studied the consequences of these laws and have applied them to 
innumerable macroscopic systems. In view of these well-known 
facts it may seem heresy to suggest that the laws of mechanics 
can be formulated in an alternative fashion. On the other hand 
as will be seen, the proposed reformulation makes use of hints 
writers in the 


appearing in books and arti les by a number o 


field of mechanics, among whom may be mentioned D’Alembert, 


Appell, and Firestone Such a reformulation seems worth at- 


tempting for two practical reasons: (a) because dynamical equa 
tions formulated in a different order can be particularly useful in 
the engineering sciences, and (6) because a new point of view on 


an old subject helps to clear up obscure points. 
NEWTONIAN MECHANICS 


It is well known that Newton's laws in the usually stated form, 
due actually to Euler, apply to the dynamics of particles; a par- 
ticle being defined as a body in which rotations and deformations 
are insignificant in the problem under consideration. Studies 
of the invariant properti+s of the laws show that if they are valid 


for one set of reference axes, they are valid for a second set of 


axes if, and only if, there is no relative acceleration between the 
Hence Newton's laws are valid only for a very special 
The 


demonstration that any given choice of reference axes is indeed 


two sets, 


set of reference axes which are often called “inertial axes 


inertial presents one of the most difficult problems in Newtonian 
mechanics even to this day. Be that as it may, any acceleration 
of a particle measured relative to a set of inertial axes is called an 
absolute acceleration This concept of an absolute accelera- 
tion plays a principal role in the laws of classical mechanics 
Let a mechanical system be defined as a collection of particles 
which, in general, interact on each other. We admit the possi- 
bility that excitations may be applied to the system by external 
devices. One of the great triumphs of Newtonian mechanics is 
that the absolute motions of the particles within a system can 
Head, Applied Mathematics Branch, Mechanics Division, Naval 
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be computed from only a knowledge of the externally applied 


forces and the nature of the interactions. On the other hand, if 
a know ledge of the internal forces is desired, the motions of each 
particle either must be computed or be known by some other 
in many prac tical cases, a knowledge of inte rnal 


that 


means Since, 


interest, it seems unfortunate these 
forces cannot be calculated directly from a knowledge of external 


This deficiency will 


forces 18 of primary 


excitations and the natures of interactions 
be removed by the methods presented in the following discussion 
they what 


nodal in character since specily 


This characteristic 


Newton's laws are 


happens at a point in the system is most 
easily understood in terms of D’Alembert's principle. Let a 
parti le of mass m be acted on by n-forces, f,, fe, f f Then 


by the second law of motion 


where a is the absolute acceleration of the particle 
With D’Alembert, define a force i, 30 that fy 


have 


then - 


f + ty = () 


The force f,, is usually called an “inertial force’ and we have the 


principle, due to D’Alembert, that the sum of all forces including 
the inertial force acting on a particle is equal to zero Ap 


so well known that they 


point 
plications of this principle are need no 
further comment. The important point is the analogy between 
this principle and Kirchhoff's current law which may be written 
the sum of all currents ent ring 4 Junction point ts zero 

Let us pursue this analogy further. Those familiar with mod- 


ern electric network analysis are aware that a knowledge of the 
applied excitations and the nature of each element in a system 
together with Kirchhoff's current law is sufficient for an analysis 
This 


simultaneous equations in which voltage drops are the dependent 


se-called nodal analysis yields a set of 


of the network 


variables. Should a knowledge of the currents flowing through 


elements be desired, a second calculation must be made Sup- 


pose, however, primary interest centers on the currents flowing 


elements Is it necessary first to compute voltage 


The there is an alternative 
cedure—the this kind 
made of Kirchhoff's voltage law which states that the sum of the 


through 


drops? answer is “no,” for pro- 


mesh analysis."’ In of analysis use is 
voltage drops around any closed path is equal to zero. This is a 
circuital concept and one which focuses attention on the con- 
tinuity of currents in a network. Dependent variables are now 
circuital currents and voltage drops are computed, if at all, as 
a second step. The important point for our purpose is that a 
shift from a nodal to a circuital basis actually changes the nature 
of the dependent variables used in the formulation 

ourselves, by 


Newton's 


Returning to the mechanical system, we ask 


analogy, the following question If we shift from 


laws, which are nodal in character, to a circuital basis, will the 


order in which motion and forces are computed be reversed? 
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under of this The | mand side 
D’Alembert forces is assume 


r expression 0 


ft the continuity of 
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schemati 
ing discussion sho 
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device like 
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that the applied force 
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, ] ’Alembert for 
tion of D’Alembert \ 


fp, which we shall 


f the partick 
gy with the assum; 
uivalent to a force 

5] then becomes 


t Equatior 
We now se 


ma is ¢ 


assumptions 


a transmitted gain have mect 

D’Alembert force is direction 
ption of D’Alembert 

Three postulates 


ind u Is taken as 
of the continuity of 


ple torces, 


Hereafter, it will assumed that the 


real transmitted force. It has one quality analogous to that o 
Maxwell displacement current " continuity of force 
through a , the displacement current assures con measurable 
tinuit current through a condense lied unless tl 
essential difference | 

asurable 


one 


tually 
ep ui i 
inferred onl 


with interactions be- 
“he first t 

n) external torces 
laws of motion 


Let there be 


ind let F, represent one of these forces 





AN ALTERNATIVE FORMULATION THE LAWS OF MECHANICS 
tal nature of the formulation ly r : ependent variables 
stem, it is necessary to have three é | I fores nd Assi them to i 


iditions 1) It is necessary to know the — loops i lection ua mace 


1 
iture of ipplied to the system. Typical cases int t appe c hese be 
ure t - me either a transmitted force or s I h rt the force 
motior i1ust be possible to express the relative spring ; force transmitte 
oss a interact 1 terms of the transmit 1 f A chow ) ‘ th these inte 
of a particle rela 
ihed in terms of the t - * next step the torn ition 18 to w 

easity that any sum of 1 f the re y ‘ ies around ¢ 


p is equal to zero must be 


where € is ompressibility of the spring and B is the 
of the co nt of ty of the dashpot 


are substitutes ition 11] the follov 


Inertial 
Reference 


forces are show! lotted lit \ ve beer 
ease in reading the drawin practical embodimer 
system might be a tion-type vibrator in whict 
centric mass, mM, 1 tform being vibrated, 
prescribed motior the itform relative to the 
and the platfor: 
spring which i 
The diagram Lin ometrically independent loops 


or meshes, which means th he equations for the system must 
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With a little experience it is easy to go from the schematic dia- 
gram directly to equations like Equations [13] by simple in- 
spection 


CONCLUSION 


It is clear that the circuital type of formulation has led di- 
rectly to dynamical relations involving only excitations and trans- 
mitted forces. Again it is emphasized that no direct specifications 


of an inertial reference are required 


only an 
motions be 


admission that 
desired, they can be 
obtained by inserting the values of f, and f, into Equations [12 


one exists. Should relative 


Two remarks are in order in this connection: (1) Any velocity 


ola particle relative to the inertial reference can be computed 


only to within an additive constant since the initial values of the 


indefinite integr ils are not known in general On the other hand 
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absolute accelerations can be computed exactly, and if observed 
values of acceleration with respect to arbitrarily selected co- 
ordinates do not agree with the computed values, then indeed the 
co-ordinates selected are not inertial. (2) The computed values 
of any relative velocity across an interaction will be observed in 
practice if it is referred to co-ordinates which are not rotating 
This is a much less restrictive requirement. 
In conclusion, then, there are alternative ways in which a me 


chanical system can be analyzed—one based on Newton's laws, 


which leads naturally to equations in kinematical variables, the 
’ 
| 


other, a circuital type of analysis which leads to equations in force 
variables. The most efficient one to use in any given case de 

pends upon the nature of the system and on the results desired 
It does appear that the circuital method is the natural one to use 


when interest is focused primarily on internal forces 





A Variational Principle for the Mesh-Type 


Analysis of a M 
By R. A. TOUPIN,' 


The mesh-type equations for a dynamical system con- 
sisting of N-mass points are derived from the Newtonian 
form of the equations of motion. The appropriate La- 
grangian for a variational principle is established for these 
equations. 

Ss 


Trent,? 


mechanical 


by methods 
the 
interest is primarily in the solution for the internal forces 
the 


The usual procedure 


pointed out there are alternative 


of analyzing a system. Frequently 


of a mechanical system, the actual motion of components 


being of secondary concern however, 
is to solve the problem first for the detailed motion of the system 
and then to compute the internal forces by means of some func- 
kinematical variables and the 
The Newtonian 


equations of motion leads directly to this procedure 


tional relation between the inter- 
the 


It is pos- 


nal forces of the system formulation of 
sible, however, to derive a set of equations, sufficient in number, 
to determine the internal 


kinematical solution 


forces directly without recourse to the 


These equations follow from the Newto- 


nian form of the equations of motion. Further, it is possible to 
give an analytical development of the new system quite analo- 


gous to the ordinary « xposition in analytical mechanics. 


Tue Mesu Equations 
Consider a dynamical system consisting of N-mass points sub- 


ject only to interaction forces Let f* be the ith component of the 


force acting ob m du 
» “ 


to the presence of m The Newtonian equa- 


tions of motion for such a system are the 3.N equations 


where 


the f* 


re 


Ordinarily, are some specified functions of £', the relative 
wo 


position vector of the two particles, or in the ge neral case, a func- 
tion of E* and the time explicitly 


ue 
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Newton's third law requires that 


mo mm 
The usual procedure is to substitute the Functions [3 
9 


into [1] and 


and obtain differential equations for the z The mesh equa 


ms 


tions for the system can be obtained as follows 


Integrating Equation [1 vields 


once 


The constant of integration is fixed by the initial conditions 


we define 


Then Equation [5 


Now form the differen combining Eq 


uation 


7] two at a time 


There are 3N(N quantities J‘, but due to their antisym- 


“ 


metric property, only BN(N 1)/2 We shall 
We now 


may be inverted to yield equa- 


are independent 
take as the independent variables those having yu 
the 


form 


Pp 
assume that Functions [3 


tions of the 


We do not discuss the general class of functions for which this is 
But the 
Inversion 
R 


possible for usual 
the 


u 


functions occurring in mechanics, 


9 readily obtained 


is Substituting Equation 
It 
should be pointed out that here we have considered a system of 
to 


into vields the mesh equations for the systen also 


mass points subject only interaction forces. Classical me- 


chanics separates the treatment of external and internal 
of such a 


forces 
system In the formalism developed here, external 
forces can be introduced by choosing the system so as to include 
The h 


then be regarded as infinite in the final equations 
If we define the scalar function 7 


the sources of the external fields mass of su 4 source can 


it 
given by 


is easy to verify that the left-hand side of Equation [8] is 
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The Newtonian equations of motion are not covariant under the 
i 


1 equations take the form transformation, Equation [22], unless 


aT 0 
ol* in which case the transformation is called Galilean. It is inter- 
esting to note that Equations [8] do not require the sec ond of the 
The relative velocities play the role of generalized forces in Conditions [23) for their covariance They remain valid in a 
Lagrangian mechanics. If the forces are conservative, we may co-ordinate system that is moving with arbitrary translator) 
lefine a potential function W such that motion relative to an inertial system. It is this property that 
illows the following inte rpretation of the mie sh type equations 
on Substitution of Equations [7] into [8] yields the identical re- 
itions 
24 


of Legendre’s dual 


isterisks 


substitutions 


iv define the function 4 ilogous 


lissipation function 


become 


ne 


function of J 


Lagrangian fun 


With this 


seen with the 


The equations daisct 

uding the interpre tion of t V proy ies Of Impulse 
ind = foree It I t deriva 1 given here wil 
further clarify r connectior ith t wtonian equations 
f motion for a 
zation. The author 
Ft Died onts — tion of this work and 
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The Limit Design of a Transverselv 


— 


Loaded Square Grid 


BY JACQUES HEYMAN 


An earlier paper’ discussed the derivation of a break- 
down criterion for beams subjected to combined bending 
and torsion. The present paper deals with the limit de- 
sign of grids, formed by two sets of parallel beams inter- 
secting at right angles, subjected to transverse loading at 
the joints. This form of loading introduces both bending 
and twisting moments in the beams, and modes of failure 
under these combined actions are investigated. Exact 
solutions are determined for some simple grids, but gen- 
eral methods are demonstrated which lead to upper and 


lower bounds oh the collapse load. 


Space Frames 


S in the previous paper,’ the type of space frame to be in 
vestigated has members which lie all in the same plane 
Thus bend- 


ixes lie perpendicular to the plane and sl 


all loads acting perpendicularly to this plane 
ing moments whose ear 


forces in the plane are zero. Any member of the frame is then 
acted upon by shear forces parallel to the applied loads and by 
two moments whose axes lie 


VM and a torque 7 
iusing breakdown at a hinge 


in the plane, that is, a bending 


moment It is assumed that the shear force 
plays no part in « this usually will be 
the case, but attention would have to be paid to example Ss in 


which exceptionally large shear forces were encountered 


The customary structural assumption is made that deforma- 


tions are small compared with the over-all dimensions of the 


frame, so that it is admissible to satisfy the equilibrium conditions 

in the undeformed rather than in the deformed configuration 
It was shown in the previous paper’ that a plastic hinge would 
V and a 


be formed in a box beam carrying a bending moment 


T when 


torque 


pure bend- 


Mis the lim 


ing Later in the present 


full plastic moment) in 
the slightly 


where t moment 


paper, more general 


eriterion 
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PROVIDENCI 


will be used, where @ is a (known) number and corresponds to the 


shape of cross section employed. In the work immediately fol 
lowing, however, @ wil! be taken as unity in order to simplify the 


calculations, so that 
vy? 


rhe methods proj l are in no this restriction 
on @ 
no twisting if ding moment V at any 


For plane 


section must satisfy 


und | quation orresponding re 


space fr 
V,? 


with other work,* ® the 


moment M> will be 


weight pe 


wr “ 
with limit taken 


tional to VM 


so th a frame consisting of 


terial consumptior 1) to some seale by 


\ 
> Vl, 


limit moment and length 


Tue Square Grip 


s square opening of side length L as shown in Fig. 1 


a a 


This opening mig rrespond to tion | en bulkheads 
and de yuilding. The 


probiem is to cove ) ng 


KS in 
with a grid of beams in order to 
support a | niformly distributed over the area. For the 
sake of simy t it ! ” taken as 0 


entrated loads of y r t the jou the geri If 


quai con 
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there are n equally spaced beams spanning in either direction, 
then there will be n? joints, so that 


] 


In addition, the elementary span length will be denoted by J, so 
that 


g 


The general problem will be approached through the investiga- 
tion of very simple grids consisting of a smal! number of beams. 
For all the grids under discussion, it is assumed that the equally 
spaced parallel beams have equal limit moments which may, 
moments of the beams 


from the limit 


This restriction, as will 


however, be different 
spanning at right angles ” seen, does 
not affect the generality of the proposed methods 


SinGLe Beam Spannine Eacu Way 


Consider a grid formed of two beams of section moments M 
ind M, arranged symmetrically in the opening and intersecting 
it the center at right angles The ends A and B of the beams are 
encastré, so that the limit reached at those 
rigid, allowing bending moment 
beam to the other The 


It is easy to see that the mode of 


moments may be 


points, and the joint O is perfectly 
ind torque to be transmitted from one 
P acts at the joint O 


il grid corresponds to the mechanism of 


load p 
of this symmetri 


failure 


one degree of freedom shown in Fig. 2. The eight hinges (shown 


M> 


2 


Fie. 2. CoLttapse Mope or Sincie-Beam Grit 


as dots) are all formed under pure bending, so that the problem is 
in fact one of plane failure in this instance. 
By considering the vertical equilibrium 


VM, M, 
4 +4 
l l 


P= 


This value for p may also be derived from the principle of virtual 
work. If the deflection of the center O of the grid is 6, then the 
angle change at each of the eight hinges is 6//. Now the plastic 
work done in a hinge is equal to the angle change multiplied by 
the limit moment at that hinge, and equating this to the work 


5 Ff) 
ps am( + 4M (’) 


The material consumption 


done by the load Pp 


[10] 


which leads directly to Equation [9 
in the grid is given from Equation [6] by 


M, + Ms [11] 


WeQL 


Hence, using | quation (9 and remembering that L = 2/ 
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WwW l L , 
= yi, 12] 

- 8 
The ratio of W/P is a measure of the material consumption per 
unit load, and it will be noticed that this expression is independent 
of M, and M,. This means that considering the problem as one 
of design, any values of M, and M; so chosen that the grid fails 
under a unit load will lead to a material consumption that is con- 
stant. In particular, if M@, = M, = M 


l 
My = PL 
16 


Two Beams SpanninG Eacn Way 


2 beam grid shown in Fig. 3 
The loads p = '/,P act at 


The mode of failure of the 2 x 
again introduces no twisting moments. 
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be 
— +45—+— 


hail 


A 


| 
; 


CT'wo-Beam Grip 


the joints O, and hinges are formed at the two ends of the beams 
AO, BO. (The beams forming the square OOOO are also sub- 
Vertical equilibrium or 


jected to the limit moments M, and M, 


virtual work gives 


1.€ 
The material consumption is given by 


W = 2L(M, + M, 


and in particular, if M 

V 
It will be noticed that the material consumption per unit load is 
reduced in the ratio °/» when compared with the previous case 


Turee Beams Spannine Eacu Way 


. 3 beam grid, twisting moments are introduced it 
breakdown 


Witha3 
the collapse, and use must be made not only of the 
criterion (Equation [3] ), but also of the orthogonal flow relation- 
ship which may be derived from it.* If the pl 
bending angle at a hinge is denoted by 8, and the plastic chang: 


istic change in 


in twist at the same hinge by @, then the flow relationship corre- 


sponding to Equation [3] is 
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It will be useful at this stage to derive expressions for the moment the internal and external work is always greater than the true 

and torque in terms of the changes in angle at a hinge. Solving collapse load. This kinematic principle gives a valuable method 

Equations [3] and [18] simultaneously for determining an upper bound on the collapse load A mecha 
1 } : PI } 


nism (of one or more degrees of freedom) is constructed by the in 
sertion of yield hinges at arbitrary locations of the structure, and 
the work done by the external loads during a small displacement 
equated to the work dissipated in the hinges. The load deter- 
mined by this equation 1s then greater than the true collapse 
load. Thus, for the example cited, it is known from Equation 
In addition, the total plastic work done at the hinge is (4/8 + !22] that the true collapse load p, is less than */s(Mo/l) (or equal 
76), since M and T remain constant during collapse, so that to this value, if the assumed mode happens to be correct 
In order to find a lower bound for p,, the frame may be exam 
Plastic work at hinge = My 0/6? + 8? 20} ined statically under the load p */(Mo/l). The various 
torques, shears, and bending moments are shown in Fig 5 where 
The three-beam grid is shown in Fig. 4. The elementary span 


length is / = '/,L, and the load at each of the joints is p = '/,P +. 
It may be shown that minimum material consumption is ac hieved 
for such a grid when beams of equal section span in either direc- 


@! 

> yale 

tion. For this reason it will be assumed that WM, = M. = Mg, so “._®) rs ™ 
yrne 


that Equations [3 and [18] apply at every hinge in the grid ~\ 


cm ) 
Owing to symmetry, it will be necessary to consider only the por- => : 8 Me 
) ’ => 4 
ga ts 3 i 


tion of the grid enclose the dashed line in Fig. 4 
2 _-; 
:.» i i? cs fort 
a (Ko-™p ~ 
“Sad 3 t 2 Me =) ~ 
» (0) 4* ” 5 ; 
ae la 


I > Staticat Anacysis or Tuoree-Beam Grip 


the single redundancy of the central portion ot the frame has 
been introduced as the unknown bending moment M,. It will be 
atelLe seen that the moment at O is 


11 
( “M 2M») 


while the moment at the end C of the beam CD is 


4 Turee-Beam Grip 


An exact analysis may be made for this rather simple grid, and 
the collapse load is given by 


Mo 


» = 229 9 2 
} eae 21 (5. + Mo) 


However, it is possible to find close bounds on this value p, of the 


Suppose the yield criterion is exceeded at both these critical 
collapse load by the use of two principles, which perhaps are best — ° segue on . _ — 


‘Hustrated by ex umpk s of thei appli sth points In the same ratio k, that is 
Suppose collapse of the grid occurs in a way similar to the pre- 

vious case, that is, beams AD and BC develop hinges at both ends 

while the central portion of the grid behaves as a rigid body 


-(") 


Then, by virtual work 


V 2 = k?M,? 


These simultaneous equations are satisfied by M, = Mo, and 
this value of VV, gives the bending moment distribution shown in 
Fig. 6, from which the “worst” bending moment (Mo) occurs 
at the two critical points considered 

Now, suppose all the values in Fig. 8 to be multiplied by the 
I:quation [22) was derived from the fact that, during collapse, the — factor 5; a bending-moment distribution results that is an 
work done by the external forces on the structure is equal to equilibrium with a load of 
the work dissipated in the hinges. Now this is true only if the as- : 
sumed mode of collapse is the mode that actually will occur; if ; (") M SM 

; 


= 24) 
l ” t 


the assumed mode is not correct, it has been shown by Drucker, a 
Ciree » am *rager® 7 R » § ‘te “i b quat r 
ireenberg, and Prager, that the load determined by equating and which nowhere violates the yield condition. Such a distribu 

“On Limit Design of Beams and Frames," by H. J. Greenberg tion has been termed statically admissible, and it has been 
and W. Prager, Proceedings of the ASCE, Separate No. 59, February shown,*? that the load in equilibrium with such a distribution is 
1951 


The Safety Fact ton Blleatio-Plestic Body in Plane Strain,” always less than the true collapse load. The principle thus gives a 
1e@ Sale actor an astic- astec ane Strs 


by D. C. Drucker, H. J. Greenberg, and W. Prager, JouRNAL oF method for determining a lower bound on the collapse load; if, 
Apptiep Mecuanics, Trans. ASME, vol. 73, 1951, pp. 371 378 for any given load, a distribution of moment, torque, and shear 
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under the load p, it will be found, in general, that it is impossible to 


satisfy equilibrium conditions at the joint , the total load at each 


joint being either in excess of or less than p 


extra (positive or negative forces are req 1 
order to produce the originally assumed collapse 
appre 


The significance of these extra forces st 
vctually performing a calculation beam grid 
Suppose as a first tria the nt 


d O are given by 


is, joints ADC and BCO lie 


central joint has been tak 


hanism of 


2.135 
was estimate 
ier to improve these bounds, extra yi ld hinges must be in- for equilibrium, the torques 
At first sight it would seem be equal, the bending moments must 


the two critical points 
n Fig. 7 (drawn looking occur woth points (Equation [3 


hanism shown 7 
uu endent de e ratio B/@ the same at both point 
Table 1 satisfy this condition appr 
emphasized that this extra stati il condition is 
the applic ation ol the method 
ment in the value for the upper 


und D that total 0.19, 
The total work done 
ing this to the work done by 


ianism is ysis of the fram 
voth of Equa- oint gives 
al equilib- £ le, appears that a load 
joint D to produce the defile 


only one 


he collapse nearly double tl 
lution which issumed defle 
equations are 


values of the 


ms of the one 
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tal load at joint C is less than the actual load, so 


that the deflection should be increased. 


sane Way, the 
A better bound for the 
collapse load should be obtained by making these changes in de- 
flection 

In work nore my ited examples, it is found that the 
rs a marked resemblance to a 


out-of-t 


ijusting ce ions be 


process of a 
relaxation process 1 a reduction in the lance 


force at one jou ! ‘ iereased errors 
joints. In fact 
tion of the defi 
deflection of joint ¢ 


Suppose ther } lor a sex 


neighboring 


1 simple trial shows that reduc- 
int D immediately requires a reduction 
latter 


in the instead of an increase at the 


joint. ond trial, deflections 


0.60, 6 1.00 


ts ADC lie on a straigh 


iis gives the value of the 


nd the yints indicate that the 


deflecti 00 smal] deflection of C is too 
large. The chat m made and the process repeated, but a 
lower bound oul ermi on the collapse load before 


work small changes are made in the 


further 
bending momet fe s, an equilibrium sys- 
tem may be cor ‘ the gris vd e load given by 
examination that these 
of the 


L.OL6 


Equation [29 is found on changes 
assumed 
Thus a 


1 by reducing the load in the ratio 


cause the y nditior ” Violatec me 


hinges, the lation being in the ratio 


lower bound may be determine 


1/1.016. leadi o the continued inequality 


] per cent 


certainly close enough for any practical 


and these bounds are 
purpose 
Before dealing with the analysis of more complicated grids, a 


summary may be made of the proposed procedure, based upon 


the results obtained in the foregoing: 


1 Insert yield at a large number of points in the grid 
producing a me sm of many degrees of freedom. The hinges 


should be place t all the sections at which i ictual 


hinges might o 


ollapse. 
but reasonable) def 


rmine the corresponding changes in angle at 


the work dissipated in the hinges to the 


work done th nal loads gives a value of the load which is 


mn excess of t load 


3 t l t out-ol-balance forces at ea 


necessary assumed deflections out-of- 
balance for¢ tion as the actual load at a 


same direst 


joint estimated as too large nd 


lorce wi mai it 
it the upper bound is a good estimate of the col- 
Adjustments are 


lapse made 


made it ind twisting moments in 


4 TRANSVERSELY 
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order to produce an equilibrium system, from which a lower bound 


may be determined 


inalysis be made in the 


It is not recommended that a ical 


first few trials, as this usually involves more labor than does the 


determination of a good upper bound 


More Compuicatep Grips 


information of the possible saving it 


In order to obtain some 
material when a grid formed of a large number of beams is used 


grids were analyzed, composed of four, five, and nine beams span- 


ning in either direction. The section moments of the two sets of 


beams were again made equal, and the procedure outlined in the 


foregoing was used to determine bounds on the collapse load and 


hence on the materia) consumption per unit load. These dimen- 


vaterial consumptions are plotted in Fig. 9. Bounds 


SIODIESS 


were obtained for the nine-beam grid to within +5 per cent, 


which is probably close enough for any design purpose. However 


were still and 


The 


hinges chosen to determine the upper bound is 


the out-of- large at this stave, 


quite 
the next step to red he upper bound was clesrly defined 
arrangem 
the grid in Fig. 10 

tly reduced as the number 


show! quarter of 


The rate of saving of material is gre 
of beams is increased, as may be seen from Fig. 9; a full investiga- 


tion of the economy effected by using a large number of beams 


nit mat 
mot 


ion fo 
beam grid. 
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would have to be based upon such factors as the cost of fabricating 
For ex- 
ample, the use of a nine-beam rather than a three-beam grid would 


a complicated grid and the increased costs of erection. 


give a saving of about 25 per cent in raw material; however, this 
economy might be completely lost in the increased cost of welding 
81 as against 9 joints 


More GENERAL BREAKDOWN CRITERION 
In conclusion of the present work, the more general criterion 


M? + aT? = M,? 


[32] 
will be applied to both three and four-beam grids. In order to 
obtain some measure of the material consumption, it should be 
noted from Equation [32] that the variation of @ does not affect 
the value of the moment to cause breakdown when no torque is 
present. On the other hand, the value of the breakdown torque 
in pure twisting is given by 


Mo 
Va 


[33] 


Thus the variation of @ implies that the torque to cause break- 
down is varied, while the limit moment in pure bending is held 
constant. the is infinitely 
strong against twisting, which would imply an infinite amount of 


In particular, for a = 0, section 


material in the cross section. For this reason the material con- 


sumption will be expressed as 


W = 1/45 (M, + TOI 


= I/TMA (34 


where Vv, and T 
under pure bending and torque, respectively. It 


, are the breakdown moments for the ith beam 
should be 
pointed out that Equation [34] is completely arbitrary, and is one 


} 


of several that reduce to Equation [6] when a@ = 1. However, it 
will serve to give some idea of the material consumption involved 

Exact solutions for both three and four-beam grids were ob- 
tained for a range of values of a, and the dimensionless material 
compared with the material consumption for a 


1) are plotted in Fig. 11 


consumptions 
three-beam grid with a = The general 
shape of the curves indicates that a section weak in torsion will 


give the most economical material consumption. However, it has 


been pointed out that Equation [34] is arbitrary, and a detailed 


investigation of the amount of material necessary to produce a 
given torsional limit must be carried out before any definite con- 
clusions may be drawn 


SUMMARY 


This paper describes and analyzes the type of square grid that 
might be employed to support transverse loads, and a method is 
proposed which leads to arbitrarily close bounds on the collapse 
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Material consumption 


% beam grid 
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MATERIAL 
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ConsumptTions ror Dirrerent Breakpown 


CONDITIONS 
load of such a grid. Although discussion has been made only of 
grids formed by the intersection of regularly spaced beams, rec- 
tangular, or indeed grids of any shape, may be investigated using 
identical methods of solution. 

In effect, the determination of an upper bound on the collapse 
load is an extension of the kinematic principle of Greenberg and 
Prager,* but applied to states which are not necessarily in equilib- 
rium with the applied loads. Equating the internal and external 
work still gives an upper bound on the collapse load, however 
For the grid problems, hinges are assumed (generally in a greater 
number than that required for collapse), and a collapse configura- 
tion is arbitrarily specified, from which the 


internal and external 
work can be calculated 
Since more hinges may be assumed than are required to form a 


mechanism of one degree of freedom, the resulting bending- 


moment distribution will not be an equilibrium state; if the 
kinematic solution is reasonably good, slight adjustments may be 
made in order to obtain a statically admissible equilibrium solu- 
tion, from which a lower bound may be determined 

The conclusions that may be drawn from the minimum ma 
terial consumption investigation must be regarded as tentative 
It was shown that a grid consisting of a large number of beams is 
the most economical in the consumption of raw material, but the 
than this 


In addition, no accurate estimate was made of the ma 


increased costs of fabrication may more overcome 


saving 
terial required in a cross section in order that given bending and 
In fact, it is likely 
that with the advent of limit design methods, new structural sex 


torsional limit moments might be sustained 


tions will be designed in order to achieve the greatest economy 
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Heat-Exchanger Tube-Sheet Design—2 


Fixed Tube Sheets 
By K. A. GARDNER,' MASSILLON, OHIO 


It is shown that ‘“‘fixed’’ tube sheets may be designed VM, radial bending moment, in-lb per in 
in exactly the same manner as “‘floating’’ tube sheets with Vf, = peripheral bending moment, in-lb per in 
the same boundary restraint, provided that a fictitious ng = number of diaphragms per expansion joint 


” 


uniform “‘equivalent design pressure”’ is used in the cal- N = number of tubes in tube sheet 

culations instead of the actual hydrostatic pressure. p = hydrostatic pressure, psi (without subscript represents 
This equivalent pressure is evaluated in terms of tube-side equivalent design pressure, see Equations [9] and 
pressure, shell-side pressure, differential thermal expan- 16 

sion, and the condition of boundary restraint. The de- ? = tube pitch (spacing between centers), in 

sign factors for all tube sheets presented in an earlier = local effective pressure at radius r, psi 

paper? are shown to be well represented by very simple = radius, in 

expressions when the fundamental design parameter = thickness, in. (without subscript represents tube-wall 
xr, becomes large. thickness 

= temperature, deg F 

= tube-sheet deflection, in 

= tube load, lb per tube 


NOMENCLATURE 


7 t y hi: t s mee | = : 1! 
Phe lowing nomenclature is used in thi paper. = Br, dimensionless (see Equation 10 


a = radius of tube sheet over which hydrostatic pressure = linear coefficient of thermal expansion, (deg F 
acts, in = (2k/D)4, in, 
ber 2s = real part of Bessel function, Jo(2y/i , = flexural efficiency, dimensionless 
bei s = imaginary part of Bessel function, Jo(ry/i Ne = ligament efficiency, dimensionless 
ber’ x = d(ber r)/dr @ = tube-sheet slope at radius r, radians 
bei’ z = d(bei z)/dr v = Poisson's ratio 
width of annular expansion joint diaphragms, o = stress, psi 
arbitrary constant Subscripts 
outside diameter of tube, in 
flexural rigidity of tube sheet = Kh*n,,/12(1 
Ib 
base of Napierian logarithms = 2.718 


a = at periphery of tube sheet 
¢ = at critical radius of tube sheet 
d = diaphragm of expansion joint 
) , = t al 
modulus of elasticity, psi radia 
: s = shell 
fraction of total tube-sheet area on which pressure ; 
{ = tube 
directly ib 


INTRODUCTION 
dimensionless factors for obtaining effective pressures, 


deflection, slope, and radial bending moment of tube In an earlier paper* the author presented an analysis of the 
sheets deflection and bending moments in the circular tube sheets of 
heat exchangers of the internal floating head’ type For ex- 


changers in which two tube sheets are connected by the tubes 

tube sheet thickness, in . . 

dhiauel nore expanded or welded into them (as distinguished from U-tube 
) iry constan 

or bayonet-tube exchangers), it was shown that the tube-sheet 

tube-bundle modulus - 


Vi(d t)E,/a?L, psi per in 
2E,t.a/ E.Nt(d t) (Equation [17 dimensionless 


deflection equations are those of plates on an elastic foundation, 
this foundation being the “tube bundle”’ itself. It was stated 
tube length, ir in the conclusions of the paper that the same method of analysis 
 lengtl " 
could b pplied to tube sheets tegral with the s ting 
distance tube-sheet edges move with respect to ear h yuld be pple | » tube et integral ith 1¢ supporting 
, member . and it is the purpose of the present paper to 
other under load, in 
derive design equations for such conditions. It will also be 
Assistant Chief Engineer, The Griscom-Russell Company. Mem shown that the results are applicable, as a special case, to “ex 
ASMI 
Heat-Exchanger Tube-Sheet Design,’ by K. A. Gardner, Trans 
ASME, Jovurnat or Apptiep Mecnanies, vol. 70, 1948, p. 377 
Contributed by the Applied Mechanics Division and presented at diagrammatically in Fig. 1 Phis is the simplest and cheapest 
the Annual Meeting, Atlantic City, N. J., November 25-30, 1951, of type of shell-and-tube exchanger, and while it suffers from the 
THe AMERICAN Society OF MECHANICAL ENGINEERS 


ternally packed floating head” exchangers 
A fixed or ‘integral’ tube-sheet type of exchanger is shown 


disadvantages of inaccessability of the outer tube surface for 


Discussion of this paper should be addressed to the Secretary . . 
ASME, 29 West 39th Street, New York, N. Y., and will be accepted leaning and from differential thermal-expansion stresses, it is 


¢ 


for many performance conditions where these factors 


until July 10, 1952, for publication at a later date. Discussion re idequate 

ceived after the closing date will be returned are of littl consequence The addition of an expansion joint 
J atements jinions advance yapers are be , 
Note Statements and opinion ady anced in papers are to be in the shell overcomes the thermal-stress problem, although 

understood as individual expressions of their authors and not those 

of the Society Manuscript ved by ASME Applied Mechanies the floating-head type offers the same advantage and others as 

Division, June 11, 1951. Paper No. 51—-A-38 well. Fig. 1 also shows one type of externally packed floating- 
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and the other symbols are as used before. Thus deflections u 


| 


are considered positive when the two tube sheets bend toward 
each other and the value of W,, if positive, represents a tensile 


force in the tube and, if negative, a compressive tore The 
dimension m is the distance the outer edges of the two tube 


> 
NT 


> 


sheets move away from each other under the combined effects 


rn 


of pressure and thermal expansion. The load W,, considered 
as uniformly distributed over its associated tube-sheet element 
and combined with the hydrostatic pressures on both sides of the 


tube sheet, results in a net effective pressure on that element of 


a,] L t psi 


HHH EY 


where 


A Tube sheet > Tubes P; = hydrostatic pressure on tube side, psi 
B Stationary head (channel Baffles and tube supports . 
C Floating head Shell 1 


D Head partition f;: = fraction of tube side face of tube sheet wluding tub 


= hydrostatic pressure on shell side, psi 


1 Inrecrat Tuse-Surer Heat Excrancer, ann Br 
EXTERNALLY Packep FLoatinc-Heap EXcHANGER 


ends, on which tube side pressure acts direct] 


\ (« =) 
vead exchanger which, for purposes of tube-sheet computation, . e 
is identical with an integral tube-sheet exchanger with a per- f, = fraction of shell-side face of tul 
fectly flexible expansion joint in the shell. As in the previous side pressure acts directly 


paper, only circular tube sheets with uniform tube distribution, v / y 


supported only at the periphery, will be considered ; ( 


Metuop or ATrTack Nid —bE 


a 


Che nomenclature and fundamental assumptions of the a’L 


wriginal paper? are used. These assumptions are as follows: 
iy: _— . . ; Equation [2] is identical in form with the expression derived 
1 The mean neutral deflection surface of a uniformly per- aa 
in the original paper for internal-floating-head exchange 

orated plate under load may be considered congruent with that 
of a homogeneous plate of the same radius, identically loaded, 


shose flexural rigidity D is a fraction n, of that of the original 


Therefore the fundamental! differenti equation and its gener 
solutions are the same, namely 


undrilled pl ite. 


a“q 
2 The tube-sheet element associated with a single tube is so i 


: ~ dx* 
small compared to the total tube-sheet area that it may be con- 


sidered of differential size and hence the external forces exerted 
on the tube sheet by the tubes may be treated as continuously 
varving 

} The entire surface of the tube sheet over which the pres- 


sure acts may be treated as uniformly perforated, 


rhe first two assumptions make it possible to treat a drilled 
tube sheet as though it were a solid plate of some lesser thick- 
ness on a uniform elastic foundation. The last ignores the band 
of undrilled plate at the periphery which, in any event, is prac- 
tically nonexistent in fixed-tube-sheet exchangers. Malkin’ has 
in effect, proved that assumption (1) is valid for tube sheets 
drilled on a triangular lavout and has, furthermore, derived Hi ber 
expressions which permit evaluation of the deflection efficiency 
Nu however, he also shows that this assumption is, at best, : . , 
in which, for simplicity in nomen 
an empirical approximation for square tube layouts. (See the : : 
~ : . the original paper, H has replaced B 
author’s discussion following Malkin’s pape 

The force exerted by a single tube on the tube-sheet element , 8 
associated with it is determined by the change in its length 
corrected for thermal expansion), its original length, its cross- There are three boundary condi 
sectional metal area, and its modulus of elasticity; thus tion of the three arbitrary const 

] The balance of forces on ei 


given by 


coefficient of thermal expansion of tubes leg I 


average metal temperature of tubes, deg I 


3**Notes on a Theoretical Basis for the Desig Tube Sheets of 
Triangular Layout,” by 1. Malkin, Paper No. ! -120, presented at 
the Annual Meeting, New York, N. ¥ vernber 26-December 1! 
1950, of Tue American Society ofr Mecuanicat ENGINEERS 
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in which p is the average effective pressure acting on the tube 


sheet fac e, 1.€ 


. 2rtfspsde 


Fic. 2 Batance 01 «ces ON INTEGRAL TUBE SHEET 


2 The deflection w at the periphery is equal to zero 


} aT L 


a,T Lb 


the ratio of maximum to the average effective 
pressure, as used in the 


in which Fy is 
original paper, ie 


ber 7, 
H 
bei 74 


ber Tq 


bei or 


condition of peripheral restraint can be expressed 


3 The } 
} 


generally as the ratio of the 
at that point; thus 


slope to the radial bending moment 


BD 


i[ bei x f - bet r| H [ es r. oa : bei’ x | 
: | 


ber’ z, + H bei’ z 


a a 


Simultaneous solution 


in which 


HEAT-EXCHANGER TUBE-SHEET DESIGN—2 


FIXED TUBE SHEETS 
2E J, 


akl 


Bbé, 1 H 
| ber 7, tT bel J, + 
M, Ta 
(F*) é v 
bel J, ber 7, 
VM, Ir. 


Errect OF AN Expansion JOIN 


The parameter K is actually the ratio of the effective elastic 


shell to that of the tube bundk it is the 
ratio of the force required to produce a given strai: the shell 
to the force necessary to produce the same strain in 
bundle. It is thus the of the ability of the she 


tive to the tube bundle) to 


modulus of the 


the tube 


measure rela- 


movement of the vo tube- 


by the introduction of an expansion joint into the she {l- 


resist 


sheet edges with respect to ea h other ability ws reduced 


though there is little general agreement on the ce sigh ol 
pansion joints Feely 


approximate the effect on K for the type of joint which « 


a paper by and Goryl* may be used to 
nsists 
essentially of a series of annular diaphragms welded together 
oring 


ob- 


alternately at their inner and outer circumferences. Ign 
deflection due to pressure, the effective value of K may be 
tained by multiplying the value from Equation {17} t 


l 
, t, c y E, 
i [= E, 
where c is the width of the annuli, t, and BE, are the thickness 
and modulus, respectively, of the diaphragms, and n, is the total 


number of diaphragms in the expansion joint (normally a mini- 
mum of 2 The additional deflection due to internal pressure 
and the and the 


diaphragms still further reduces the effective value of K 


rotation of the Junction between the shell 


An externally packed floating head, Fig. 1, is for purposes 
a perfect expansion jouit and for 
There is, of 


a resistance to sliding friction through the packing gland 


of tube-sheet computation, 
exchangers so constructed the value of K is zero 
course, 
but this effect is negligible compared to the elastic resistance of 


metals 
EQuIvALent Design Pressurti 


The value of A for the internal floating-head type { ex- 
changer covered by the earlier paper? is obviously zero since 
there is no shell resistance to the relative movement of the edges 
of the two tube sheets. The 


Equation [16]) to 


average effective 
the 


pressure p is, 


therefore, equal (by wtual tube-side 


hydrostatic pressure p;, and, with this substitution, all the equa- 
tions of the previous paper are reproduced 

The important conclusion may now be drawn, and the main 
point of thix paper made, that the effective pressure, ction 


and so forth, of integra fixed 


all identical to the corresponding qu 


slope, bending moment 
tube sheets are 
floating tube sheets of the same diameter, identically supported 
M./@,) for the 


is the value given by Equation |16) rather 


(as characterized by which equivalent tube- 


side design pressure oa 


than p,. This is true because the value of C in all the design 


equations is simply scaled down in the ratio p/p 


the conditions of operation of most | ex 


Since hangers do 


not pret lude the possibility that the full hydrostatic pressure on 

‘Stress Studies on Piping Expansion Bellows,"’ by F. J. Feely 
Jr.. and W. M. Goryl, Jowrnat or Appitiep Mecuanics Trans 
ASME, vol. 72, 1950, p. 135 
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one side may be exerted in the absence uf any pressure on the 
other side, it will not often be possible to design on the basis 
of Equation [16] as it stands, i.e., on the basis of the differential 
pressure. Rather, the larger of the two following values usually 


must be used 
Kkl(a,T, a7 


(1 + KF, 


+ Kf, 


‘ 


+ klia,7 
(1 + KF, 


Up, 


which, in the absence of thermal effects, become 


( 1 + Af 
~ PNT + KI ‘) 
Kf, : 
Pa KF, [16d 
Ignoring thermal effects for the moment, it is apparent that 
the equivalent design pressure for a fixed-tube-sheet exchanger 
must always be less than either shell- or tube-side hydrostatic 
pressure, since f, and f, by their nature are always less than unity, 
F, is always greater than unity, and K is always positive. It 
is also obvious that with K = 0, a condition closely approximated 
by fire-tube boilers and heat exchangers with shell expansion 
joints or outside packed floating heads, the design pressure is 
entirely independent of the shell-side pressure, no matter how 
high it may be, and is, in fact, equal to the tube-side pressure 

exactly as in an internal-floating-head exchanger. 

On the other hand, if A is large the equivalent design pressure 
this has a 


approaches p,f,/F, or p,f,/F,, whichever is greater; 


maximum value of approximately one half the actual pressure and 
Be- 
tween these two extremes perhaps the most remarkable features 
of Equations [16c] and [16d] is the fact that the tube-side pres- 
sure is usually the fixed 
design, even though the shell-side pressure may be 


may be considerably less, depending on the value of Fy, 


controlling factor in tube-sheet 


many times 
4S high 

The effect of differential 
impressive when considered as an equivalent hydrostatic pres- 


thermal expansion becomes quite 


thus 


sure 


be 


1 
p= Kil (“* 


If no expansion joint is used, Akl which normally 


has a value of about 10° 
10~* and 1, respectively, the equivalent pressures are of the order 
Thus a difference between shel! 


since a and .F,) are of the orders 


of magnitude, +10(7;, ff 
and tube metal temperatures of only 50 deg F could produce 
effects equivalent to 500 psi pressure. It should be noted that 
the equivalent pressure found in this way is less than that ob- 
tained on the common assumption that all tubes are stretched 


or compressed exactly the same amount; on this basis 


. aT, aT, 
d= KL ( — ) 19 
1+A 


However, it is the “maximum local” effective pressure, not the 
average, which determines the possibility of damaged tubes 
and this quantity is higher in the ratio F,(1 + 


than would be obtained on the simpler basis 


or tubs 


K)/(1 4+ 


joints, 
KF, 
Bounpary ResTRAINT 


By definition, the slope at the periphery is negative if the tube- 
side face of the tube sheet is concave at that point, and the ex- 
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ternal bending moment is also negative if it tends to produce 
Since the moment exerted by the restraining 


the stationary head 
dimensionless parameter 


such concavity. 
the shell (and sometimes 
always opposes the deflection, the 
BD0@,/M,) occurring in Equation [18], which is a measure of 
the degree of boundary restraint, is always negative in fixed- 
other external bending 


influence of 


tube-sheet exchangers unless some 


moment is present in addition to that produced by the shell 


(See Fig. 3, for definition sketch Thus, for simply supported 


o A 
Mg a 





DEFINITION SKETCH FOR SiGns oF Deriections, Epoe 


SLOPES 


Fic. 3 
AND MoMENTS 
only and of shell and flange bolting, respe 


B and C show effects of shell 


tively 


M, = 0), 


for pertec tly clamped sheets it equals zero and, in the absence 


plates, this parameter equals minus infinity (since 
of any other restraint than the resistance to deflection and bend- 
ing of the shell, other values lie between these extremes 
However, in cases where the tube sheet is extended to form a 
flange for the tube fluid heads, a negative external bending 
moment is applied which may be sufficient to make (8D6,/M, 
positive. The effect of externally applied bending mo- 
ments is rather complex, since it depends, among other things, 


these 


upon whether the equivalent design pressure is positive or nega- 
tive and whether the value of 2, is greater or less than 2.65 
the tube-sheet stress may be alleviated rather than intensified 
even in their absence, the complication of the combined effects 
of discontinuity stresses due to shel! and head pressures and dif- 
ferential thermal expansion between shell, tube sheet, and head 
makes the determination of boundary restraint an uncertain 
factor in the estimation of tube 
In view of the foregoing, therefore, 


his judgment and experience in deter- 


-sheet stresses 

it is suggested that the 
designer be guided by 
mining the relative fixity of the tube-sheet periphery as betwee 


clamped and simply supported 


Asymptotic EXPRESSIONS FOR Desicn Factors 
In the previous paper,? the maximum effective pressure, de- 
flection, slope, radial bending moment and stress were given in 


the 


following terms 





GARDNER—HEAT-EXCHANGER TUBE-SHEET DESIGN—2 FIXED TUBE SHEETS 


TABLE 1 DESIGN FACTORS 
20 Dérw 
Fg Fe £ = Fe . — Pu 
Pp pa‘ pa? 
Simple Clamped Simple Clamped Simple Clamped Simple Clamped 
1.0000 1.0000 0.0636 0.01563 4) OO57 0 02405 206. 0.125 
1.0046 0.0628 0.01554 4).0047 4) 0239 § 0.1245 
0 060 
1.396 74 0.0531 01445 4) 0825 0218 ‘ 0.1200 
0.043 
2.253 17 ). 0316 01112 0546 0175 0.1056 
0.0150 00691 0321 O1ll 0 08560 
0.00719 00378 0210 00719 0 06945 
0.00390 00206 0144 00487 0405 0 0581 
0.00240 00354 0345 0.0500 
0 00159 00799 0 04387 
0.00110 02665 0.03913 


0 000803 00510 0 02384 0 0352 


= pa’ y 


op 


Ne 


Values of these design factors were tabulated and graphed in 
most cases over the range, z, = Oto 10. For convenience’ sake 
they are repeated here in Table 1 

With the considerable reductions of equivalent design pres- 
sures which can be achieved for fixed-tube-sheet exchangers by 





the use of Equations [16], [16a], and so forth, the necessary 
tube-sheet thicknesses decrease and the values of rq frequently 
will be found to exceed 10. Fortunately, for larger values of 
r,, the special Bessel functions employed in this analysis are 
well represented by single terms of asymptotic expansions as 
shown in the Appendix, and the expressions for the various de- 
sign factors become quite simple 
For clamped plates, as 


Fic.4 Variation or Errective Pressure at Rapivs or Maximum 
Dertection With 24; Urper Curve, Crampep; Lower Curve 
SIMPLE 


Tuse-Compression Design Factor 


The load on any tube, from Equations [1] and [2] is given by 


in which q is the local effective pressure at the tube-sheet radius 
at which the tube is located. The maximum load, as can be 
seen by examination of Fig. 8 of the previous paper,? occurs 
at the periphery and its value is 


4 


In ordinary floating-he ad-ty pe exe hangers, this is a compressive 
load if the higher pressure is O1 the shell side and a tensile load 
if on the tube side. In the type of exchanger under consideration 


Fy ‘ here, however, the outermost tubes are always in tension regard- 
less of which pressure is higher, unless the thermal-expansion 
effect is sufficiently great to overcome the pressure effect. This 


The error in these approximations at z, = 10 ranges from 0 


can be appreciated better if the expression for the equivalent 
to 1 per cent for clamped plates and from 1 to 4 per cent for 


: pressure, p, is substituted in Equation [34] thus 
simply supported plates, the maximum in the latter case occur- 
ring in Fy with an approximate value of —0.0228 as compared Ww ma?) pl, f.) + fp, — F KkI(a,T, — a,T,) 
to the calculated value of —0.0238. - \ 1+ KF, 
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This can become negative on! 


FKkL(o1 ip, 


n expansion joint or an outside packed floating he 


ill fo other 
In tension 


or the left-hand member is sm r any reas 


outer! 
Insuffi 


nost tubes must always be 


ient attention was paid i 


n the original 
he 


paper*® to t 
riterion governing t ibility of tube bun 


»p the support attributed to it in this method o analysis 


I 


failure of some ring of tubes under a « 
} h the | 


il 
tubes may 


itl 
numeri ol 


} 
Value 
eed the 


bundle, 


the t 


nsile load 


load 


that 


outern compressive 
pou terior it is Pp 


Therefor 


negative 


yssible th 


V Lite vdaditior 
h Maximum ef re 


sure 


tabulated in Table 1 and shown graphically in 


ibove 10, F. is ipproximated ) 


clamped 


simply supported 


FLEXURAL ErricteNcy 


to utilize any of 


the foregoing 
which is defined 


it is necessary to be 
inte 7, 


As 


Ni(d 


leflection efficiency” or “flexural 


rm not defined quantitatively. 


that th llowing values be | 


efficiency,”’ n,, 
It is suggested 
used pending experimental evalua- 
tion o »accurate Vv ilues 


for triangular pitch 


{0a 


for square pitch 10/ 


made 


to Malkin’s 


wr for the justifieatior 


Aper und the 


f these 


l discussion 


expr 


EXAMPLI 


tube-sheet 


ur thi 
ID, with 260 
triangular 
The shell 
ill parts is 30 x 


kness for an integrally 


IS BWG 


welded steel] 


heat exchanger, 20 in 


12 | 1 


Ss centers 


ein 
steel tubes ft long on 1-ir pitch, and 5 semi- 
support 


the 


design 


plates on 27-in 


thickness is */, in., 

psi, the allowable 
ill parts is 12,000 psi, and the coefficient of 
mansion is 6.7 K 107% per deg | 


elastic modulus for 


‘ 
stress tor 


thermal ex 


Sec 


al bucklis 


a 


F, 
Using the ordinary TEMA formul 


tact 


quivaient-pressure ¢ 


APPLIED MECHANICS 


The design conditions are as follows 


Tube-side pressure, 150 psi; metal temperature 


100 FS 
Shell-side pressure, 400 psi; metal temperature 


290 F® 
ud 


mn The fundamental constants 


( 
( 


O.040 


are 


) 07 


or 
he 260 


ile 


0.652 

” 
t 10 
2060 


1 


0.75 


) 


0.701 


10 
in 


t 260 
1e 100 
t 


O69 
0.906 ( ) 
l 


10 


x0 x . 


L.&86 10* 


2. = 10 ) 


The equivalent design pressure may be 
effects N-sidk 


t Xpansion 


12 x 0.91 


x LS6 


0.469h! 10 


ritten terms oT \ 


and 


| 
of tube-side pressure, she 


pressure 


thus 


150(1 + 0.84 72 


0 
10° 


> 


1.86 & 
0.84F, 


Even though the value of F 
be 


@ is not 
helpful to tabulate the values 


of all possibilities of pressure failure 


this point, it wil 
0.84F,)p in terms 


\V/ ) indicates which 
effects are considered to act: 


Tube 


pressure 


Shell 
pressure 


Vv 


Obviously, 
the shell 


tube pressure. 


the worst condition be 
pressure and thermal exp 
In 
iffect the equivalent thermal-expansior 
poses of illustration, this will be assur 


Some the 


will ous action 
f 
ol 


absence of 
many cases, loss ssure would 


but, for pur- 
i no » occur here 


assumption sheet t 


ind checkec 


unity 


the 


would 


tule kness 


must 
} 


now 
ve made aximum Vv would occur if 
were and ™ equ 198 psi 


10V/ 198 12000 


* These per 
ons would not 
taken 


ter atures are 


normal 


together with the 
ct of ther 
‘Standards of Tubular Exch 


ond edition, New York, N.Y 


thermal 
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which suggests a thickness of 1.25 in. as a starting point. On 
this basis 
, ox \*/s 
te $/(1.25)"/* = 4.49 
Fi= 


Fy = 


1.8 (clamped) or 3.5 (simply supported 


0.077 (clamped) or 0.055 (simply supported 


The equivalent design pressure thus becomes 
364 : sod 
- 145 (clamped) or 93 


2 3.94 


2.51 


simply supported 


ligu- 
used on the assumption 
P = 0.25 and(P —d 4 


the maximum radial bending stress is therefore 


uncertainty as to the value of the 
of 0.3 will be 
that it lies somewhere between (2’ d 


21)/P = 0.35 


Is mOrne 


Although there 
ment efficiency, a value 


6 X 0.077 
0.3 


) = 14,300 psi (clamped 


0.055 . ° : 
= 6450 psi (simply supported) 
0.3 . 

At this point the designer's judgment is called into play, for 
the thickness is more than ample if the tube sheet may be con- 
sidered simply supported, but not enough if it must be considered 


clamped. The following facts may be taken into consideration 


1) It 
}'/¢-in plate 
(b) The 
of tube 


is certain that a */,-in. shell cannot rigidly clamp a 
there will be some rotation of the junction 
condition (assuming no failure 


normal operating 


pressure) has an equivalent design pressure of 67 
psi and a bending stress of 6600 psi even if the tube sheet is 
fully cl impe d 

(c) Even if yielding of the tube sheet should occur (which 
12,000 psi 


pressure 


is very unlike owing to the safety factor in the 


stress), the equivalent thermal-expansion 


allowable 


would be partially relieved, and the condition of boundary 


restraint would approach simple support as a result. 


Assuming that the exchanger is not to operate in a cycling 


service, i.e., that it would operate under fairly steady conditions 
for months at a time, the author would be persuaded by these 
that the 1! 


considerations «in. tube-sheet thickness is ample 


and might even be further decreased, provided that no other 
part would thereby become overstressed 

With a */-ir 
0.037 


becomes 


tube sheet, I, = 6.58, . = 5.0. and Fy = 


both simply supported); the equivalent design pressure 


346/5.2 = 


x 0.037 x 67 f 10 \" 
8800 psi 
0.3 0.75 


Inasmuch as good heat-exchanger tube-expanding practice does 


67 and the stress 


not permit much thinner tube sheets for 4/,in. tubes, the re- 
maining calculations will be based upon the in. tube-sheet 
thickness 

The ability of the tubes to contribute the support attributed 
investigated. The load per tube in the 


to them must next be 


outermost ring of tubes is given by Equation [34], and in 
the tubes at deflection, by 
36]; thus, with F, = 5.0 (simply supported) or 2.6 (clamped), 


0.23 (clamped 


maximum tube-sheet Equation 


and F, = —0.52 (simply supported) or 


r X 100 
0.723 & 150 + 0.635 * 400] 


260 


1.21(5p 109 + 254) (simply supported) 


HEAT-EXCHANGER TUBE-SHEET DESIGN-2 


FIXED TUBE SHEETS 


or 


21(2.6p 109 + 254) (clamped) 


21(—0.52p 109 + 254) (simply supported) 


21(-—-0.23p 109 + 254) (clamped 


Three separate terms have been left within the parentheses to 


indicate the separate contributions of the equivalent design 


pressure, tube-side pressure, and shell-side pressure, respectively. 
With the aid of these equations and the table of equivalent 
of tube loads 


design pressures previously made, a similar table 


can be made up as follows 


Pherma 
expan- u ~ 


Clamped 


Tul 
pressure 
y + 108 157 153 
u7 + 333 + S28 
149 + 18 + 13 
+ 56 +193 +189 
+ 204 +175 +175 
41 139 139 
52 +351 341 


Shell 


pressure sion Clampeti Simple 


By inspection it can be seen that the maximum tensile load would 
occur at the point of maximum tube-sheet deflection if the tube 
pressure failed, but the other effects were still felt, and that the 
maximum compressive load would occur in the outermost ring 
of tubes if both pressures failed but the fluid streams remained 


it temperature. The same conclusions could be reached more 


were used and Equation [36 
but there 


| were 


} 


rapidly if Equation [34a 


thrown into the same form is some instructive value 


in the foregoing tabulation 
stress would be 


The maximum direct tensile 


S51 


250 psi 
0.0490.75 


x x 0.049 


Phere tube 


tube-sheet slope is zero at the 


is no bending stress in the ut this point since the 


radius of maximum deflection 


Considered as columns, the tubes have a maximum unsup- 


ported span of 45 in. from the tube sheet to the second support 
plate and could be considered as built-in at the tube-sheet end 
and “pinned” at the support plate as a conservative assumption, 


The critica! load for such a column, by Euler's formula, is 


7 ‘ ‘4 
10° & 0.375 0.326 
4 
2190 Ib 


tube length were used, considering both ends 


ritical 


Leven if the entire 


built-in but ignoring the support plates entirely, the 


load would be 


(-" ) x 30 XK 106° X — 37 26‘) = Ib 
14 / 


is compared to the maximum possible There- 


load of 351 Ib 


fore there is no danger of failure of tube support either by 


buckling or by direct. tension 
In problems of this kind there is always the possibility that 


the shell may not be able to take its share of the load without 


failure, so this point requires investigation also, From the 


equation preceding Equation [9], it will be observed that the 
tensile force acting axially through the shell is 
2rut 
a7? 
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and that this is equal L. E. Sebald, late vice-president (Engineering), The Griscom- 
Russel] Company, in the work which is summarized in this and 
previous papers, 


Cherefore the direct shell is Chief 


He also wishes to thank O. W 
Engineer, The Griscom-Russell 


I] Company f 
mission to publish this paper 


ongitudinal stress in the 


Heimberger 


or per 


Appendix 


“he Bessel ctionse ) sD 
value of p is 364/11 4 he ei function empk ved in this } ape 
115 but that this oc 


e table of equivalent design pressures it will be seen 
largest possible negative epresented 
urs in the absence of an) by the asymptotic expansions given in the 
pressure the large negative value with tube 


pressure 
acting Is 25 t 0.84 X 2.6 


9. The 


in the shell is aecordingly 


v the al 


SUMMARY AND CONCLUSIONS 


to the same assumptions and restrictions as postulate: 
previous paper? fixed tube 


sheets may be designed in 
the same Wa us the tube sl 


eets of floating-head 
\changers Dy merely 


substituting an 


equivalent design pres 
sure (lquation [16]) for the 


actual hydrostatic pressure. This 

a function of the shell-side and 
tube-side hydrostatic pressures and of the differential free thermal 
expansion between the shell and the tubes 
ol this paper, 


equivalent design pressure is 


For the purposes 
i fixed-tube-sheet exchanger need only be con 
ed is one in which the shell pressure can be exerted only 
me side of the floating head; 


thus externally pac ked floating 
head exchangers 


ome within this classification as a special case 
as do fire-tube boilers 


“ir 


ind the use of these expressions and the 
e the equivalent design pressure is normally considerabl 
than either hydrostatic 


r derivatives 

the results given in Equations [25] through [32 7 
pressure, thinner tube sheets are Inasmuch as many of the expressions el ountered InVolVe 
with floating-head exchangers this means, in ratios of these Bessel functions, the following will also be helpfu 
that higher values of parameter ? 


ountered in design. The expressions based 
uptotic expansion of the Bessel functions 


siderable simplification of the basic desigt 


the characteristic 


upon the 

result in cor 

factors as shown by 
inclusive. 

It is recognized that all the results of this paper 


»COSSOT ire 


ind its pred 
somewhat conservative owing to the neglect of the 


g bending moment of the tubes This effect 


interact Wits 
nized in the original paper in so far as its effect on the tubes 


selves was concerned but 


the reaction on the tube sheet 
gnored as negligible This effec 


jue nt pape 


t will be investigated in a 
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Bending and Buckling of an Elastically 
Restrained Circular Plate 


By H 


This paper considers the effect of an elastic boundary re- 
straint upon the deflection, moments, and critical (buck- 
ling) loads of a circular plate. The solutions given are 
based on the classical theory of plates and are exact 
They 
include, as particular limiting cases, the known solutions 
for the simply supported and rigidly clamped edge. The 
physical significance of the results obtained is discussed 


within the assumptions underlying this theory. 


in detail with particular emphasis upon the degree of re- 
straint along the clamped edge. 


NOMENCLATURE 


ing nomenclature is used in the | 


plate radius 

radius of load circle 

radial co-ordinat 

ingular co-ordinats 

a dimensionless radial co-ordinate 

b/a load position parameter 
deflection of median surface 
magnitude of concentrated load 

V., Me, Mee = radial, tangential, and twisting moment 

Y., Ve = 


O,, 08, Tr9 


radial and tangential shearing forces 
tangential, and shearing stresses 
Eh? 


12(1 


radial 
D flexural rigidity of plate = 


= Young’s modulus 

= plate thickness 
Poisson’s ratio 
restraining constant 

load intensity 

real part of 
maginary part of’ 

o* 1 Oo o? 

or? r Or r? of? 

AND STATEMENT OF PROBLEM 


INTRODUCTION 


The problem of the circular plate has been discussed by 


numerous authors for various loading and boundary conditions 
The solution corresponding to a rigidly clamped circular plate 
under a concentrated load at the center was first obtained by 
Engineer, Consolidated Vultee Aircraft Cor- 
Research Engineer, Armour Research 


' Senior Structures 


poration; formerly, Founda 
tion 

For presentation at the National Conference of the Applied Me- 
chanics Division, State College, Pa., June 19-21, 1952, of Tuer 
AMERICAN Socretry oF MECHANICAL ENGINEERS 

Discussion of this paper should be addressed to the Secretary 
ASME, 29 West 39th Street, New York, N. Y., and will be accepted 
until July 10, 1952, for publication at a later date 
ceived after the closing date will be returned. 

Nore: Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those 
of the Society Manuscript received by ASME Applied Mechanics 
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Clebsch (2) gave a solution in form of an infinite 
series for the case where the concentrated load is applied at an 


the 


Poisson (1).* 


arbitrary point, while Michell (3) obtained the solution of 


same problem in closed form. A series solution for the simply 
supported circular plate under an arbitrarily located concentrated 
load Foeppl (4 
metrical buckling of a circular plate 


was given by The problem of radially sym- 


radial 
the plate 


under compressive 


forces uniformly distributed around the edge of was first 
H. Bryan (5 


supports 


treated by G or both simple and rigidly clamped 
boundary 

It is the object of this investigation (a) to determine the de 
a circular plate under a concentrated 
the 


to find 


flection and moments of 
load acting perpendicular to the median plane of the plate, 
point of application of the load being arbitrary, and (b 
the critical loads necessary to cause radially symmetru al bu kling 
of a circular plate under a uniform distribution of radial compres- 
sive forces acting at the circular boundary in the median plane of 


the plate. In both cases (a) and (4) we shall assume the boundary 


of the plate to be elastically restrained against rotation according 
Le., the slope of a line normal to the boundary and 


is taken 


to a linear law 


tangent to the median plane at the boundary propor 


tional to the radial moment 

The results obtained are of practical significance since neither 
rigidly clamped nor simply supported boundary conditions can be 
realized fully under actual physical conditions, and thus any case 
must lie somewhere between the two 


encountered in practice 


limiting cases of the problems under consideration 
FORMULATION AND SOLUTION 


Transverse Bending Wi Boundary Restraint Accord- 
ing to the classical small-deflection theory, the transverse deflec- 


lasti 


tions of the median plane of the plate are characterized by 


q ’ 
Viw = 1) 
D 
where g and D are the load intensity and flexural rigidity, re- 
7 K . 


spectively. The bending moments, shears, and stresses are re 
lated to the deflections w by means of the following equations 


which, for our purposes, we refer to polar co-ordinates (6) 


O*u 
MV, = D ) 
oF? /_ 


° 


O*u 1 Ow 1 
Lp 
or* 


+ 
r or r? 


O*u Ou 
or 


o¢ 


? Numbers in parentheses refer to the Bibliograph 
the paper ; 
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— ( 1 Ow ) 
—2G2 
or \r 6 


Without loss of generality we may assume the point of applica- 
tion of the concentrated load F to lie on the z-axis of the system 
of co-ordinates with origin at O as shown in Fig. 1. Thus a point 
P on the plate is determined uniquely by its polar co-ordinates r 


and @ where 0 < r < aand 0 < 6 < 2m. 








Pig, 1) Eeccentricatty Loapep Cirevutar Plate 

In the particular problem under consideration the deflections 
along the circular boundary r = a must vanish, and the inclina- 
tion of the median surface with respect to a horizontal plane 1s 


proportional to the radial moment, or 


w(a, 8) = 0 


ow 
io| | = aM (a, 0 
Or ja. é 


vhere & is a dimensionless restraining constant varying from zero 


to infinity. In particular, k = 0 corresponds to a simple support, 


while k — © represents the case of a rigidly clamped boundary 
The problem then consists of determining a function w(r, @ 
a with the 


which is biharmonic throughout the region 0 < r < 


exception of the point of application of the concentrated load 
This function 
must, moreover, conform to the boundary conditions of Equations 

To solve the problem 


where it must possess an appropriate singularity. 
formulated we follow the approach 


Symonds (7) and assume a solution in the form 


} 


where Is & particular solution of Eequati m [1 


general solution of the homogeneous equation 


Viw = 0 


is biharmonic with the exception of 
The solution 


w: contains an infinity of initially arbitrary constants and is bi- 


The particular solution u 
the load point where it has the required singularity. 
harmonic in the circular region of the plate 0 <r <a. Upon 
superposition of w; and w: it is possible to adjust the constants so 
that w(r, 9) satisfies Equations [1] and [5]. 


A particular solution of Equation [1] is 


8xD 


—w, = rr logr 
F 


where w, is the deflection at P due to a load F at a distance 
from P as shown in Fig. 1. It is the singular part of the deflection 
function of a plate of infinite extent under a concentrated load 

The function w, may be expanded into an infinite trigonomet- 
ric series. With the notation indicated in Fig. 1 


r,? = pr? 4+ 2 Dr 


cos 6 
With reference to Fig. 1, | 


OP 
OF 
PF 


where z bers character 


and 2 are the complex 1 
tors OP and OF, respectively 


um 


Hence for r 


? log r; 


Similarly, for r - 


2? log r, = (r? + Bb? 2br cos 0) Re [log 


The 


co-ordinates is 


general solution of 


where 
r? log r 
log r 
R, = ile at lili a ** + Dy 


It represents the deflection of a subjected to loads 
at the boundary only. 
Since the deflections, slopes, and 


= 0 


moments must remain finite at 


dD 


ipplving to the solu- 
, and 


12] as well as the series expansions 


The remaining constants are determined 
tion, Equation [6], the boundary conditions Equations [5 
making use of Equations 
Equations [11] of the particular solution. We obtain identities 
in @ giving rise to a system of six independent equations in six 
unknowns. 
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Introducing p = r/a 


ordinate and load position parameter respectively, and solving, 


we obtain the following solution 
8xD 
F 
Srl 
i 


4y a*}? log a a’? 


log a 
8rD 
1 
a 


8rD 
L 
F 
8rD 
F 
8rD . 
( 
F 


where for convenience 


+ J) 


n(n 


The foregoing constants together with Equations [6], | 
9], and [12] yield the equation of the elastic surface 


8rD 
w(p, 6 
a’F ‘ 


? 


=',(p? +6 


25p cos 8) log (p*? 
4. §? 


W, 


« 
+ § > prlV.(p? l cos nO 

n=2 
Substituting Equation [15 
the moments 


8r 
, Mép,@) =2+2(1 + )T + 2(3 + vip 


(io +1 U) cos 6 


2 1)5? sin 8 
- 25p cos @) 


+ (1 


v 
(p? + 8? 
+ v) log [e(p? 6? 


i (n\n 1)p"W,, cos n8 


n=2 


2 
+ (1 v p 


- 25p cos 6 


into Equations [2] we obtain 


- 26p cos 8)] 


AND BUCKLING OF AN ELASTICALLY RESTRAINED CIRCULAR PLATI 


6 = b/a as dimensionless radial co- 


+8? >> nX1 ) 3 + 2Av+i l6a 


n= cont 


pV, cos n@) 


a*T 


v)d* sin 6 
25p cos @) 


21 


log le(o* 
4+ 8 log [e(p 


2dp cos @)) 


ww 
) D> (nxn l)p"W,, cosn@ 


@) sin 6 


COs 


+ 6? 25p cos 6 


psd? 2(1 


1)p"W,,, sin n 6 


6: 3 tie 
— 


p? > 


62 > ni(n + 1) 


n=2 


l)p *ip"V, sin n 0 


15 
ribed in the Appen- 


The infinite trigonometric series in Equations and [16] in- 
volving W, may be summed by methods des 
dix 

Buckling With Elastic Boundar 
under consideration is assumed to be subjected to uniform pres- 
o, at its edge 


y Restraints. The circular plate 


sure The plane-stress system therefore will be 


one of uniform compression throughout the plate, and the possible 
radially symmetric equilibrium positions of the median plane 
characterized by 


17 
must satisfy the differential equation 
DV ‘w = ho.V 


where, because of the assumption of radial symmetry 


ele 


assumes the form 


(- 


d 


v2 
dr 


18 


[; 


Thus Equation 


(43 


(dr 


d 
dr 


du 
dr 


l1 ad 
dr 


r 


where 


Integrating Equation [20] once, 


du 
r 


dr 
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where c is an arbitrary constant. However, for the case of a full 


constant c must vanish since 


(“) 
= () 
dr 0 


It is now convenient to introduce the slope 


plate the integration 


du 


into Equation [21 Carrying out the indi- 


cated differentiation we obtain 


a do 


d dr 


shall require a solution of 


wundar conditions 


the median surface with respect to a hori- 


is proportional to the radial moment 


” 


is recognized as Bessel’s differential equation, 


ich is given by 


Y, represent the Bessel 
respec 


‘ 


function of the first 


tively, both of the first order, and ¢ 

ants, 

conditions at 7 0 require that Cc 0, an 
solution, Equation [24] into the boundary 


le 
results i 
v)J (aa 
on represents a nontrivial solution of Equation 


i which are 
(buckling 


or certain discrete values of aa roots ol 


transcendental Equation [25], and the critical 
| 


s obtained trom 


Du? 


ha? 


26 


where i.e., the critical stress may be considered as a 
function of the degree of constraint along the boundary 


Discussion oF Resutts 


1s shown in Figs. 2, 3, 4, ar 


effect of 
of constraint on the deflections and moments of an 


ind v 0.3 


id 5 exhibit the 


1 cireular plate for which 6 = 0.5 


profiles along a diameter passing through the 


ATTAEM 


load 
vn in Fig. 2. Deflections corresponding to different 
f constraint are illustrated, namely, those where / 0 
l ,and k — rigid clamp 


support), A 1 (elastic clamp 


identical loads the inge ! Nection pro 


lamped plates is bounded by the profiles of 
rted and rigidly clamped plate 
the variation of the radial bending moment along a 
ng through the load point. Here 
Values of 
d plates lie between those of a rigidly clamped 
Because of 
the 
load 


This was to be expec ted because the assump- 


iwwain we consider 
midis 1, and kK— the radial moment 

iTrip» 
and simp supported plate for corre sponding points 
applicat 
tly under the 


ion of 


the singularity of the solution at the point ol 
6 = 0.5, 


concentrated load a = moments dire 


tend to infinity 
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Moments ALON DD 
Tut OED 


tions of the theory of ben 
point of application of 
Fig. 4 


circular boundary 


indicates the distributior 
As before, the 
onsidered, ar 
it the 


losest to the 


med ur d one may conclud hat the 


maximum 


radial moment ge 


Moreover 


boundary will occur at the point of the ed 
which is « 
this lov 


point of apy f sacl 
ition Is ind pe ndent 


again, values of moments 


between those of a simplh 


support 
Fig. 5 permits study of the effect of 
the maximum radial bending 


moment 


moment ts largest lor k — 


plate, and decreases rapidly with 
ported edge is approached asymptotical 
k= M, =0. 
Fig. 6 shov 


restralr 


Finally, 


stress (7, on the iimetr 


This graph is obtained from Equation 
dental I 
a continuous spectrum of 

The known limiting ¢ 
case of simple 


ind 


rigidly 
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Moment at B 
{ESTRAINING CONST, 


0.3 “ ( 


the critical stress (¢, is highe 


i decreases rapidly with a decrease in the 
restraining constant As k — 0 the buckling stress approac hes 
asymptotically the known value for the case of a simply supported 


le of buckling was considered 


first mode 


tively. It is seen that 


rigidly clamped plate ar 


plate. ¢ nly the 


ACKNOWLEDGMENT 
The author wishes to thank Dr. E. Sternberg for many valuable 


suggestions in the preparation of this paper 


Die 


St 
H 


1891 


To facilitate 
infinite series involving W,, appearing in Equations 


were summet 


STRESS ASA FuNcTEN 


First ickling mode 


BIBLIOGRAPHY 
e Mouve { 
Academie Royale 


de 


1829, p. 545 


Theor 


der Elasticitaet Fester Koerper 
1862 
ular Plat by J 


itical Society, v« 


Leipzig, Germany 
Flexure of a Cire 
lon Mather 

Biegt 
hte Bayeris 
fa Plane 


Proceedings of 


Kreisfoermigen Platte 
he Akademie W ias ol. 42 
Plate Under Thrusts in Its Own Pla 


the London Mathematica 


ability o 
Bryat 

p. 54 
Plates S. Ti 
1940 

Plane Str 


and Shells 

New York 
Proble 
Plates,”’ | 

cs, Trans. ASME, vol 


N.Y 


Appendix 
the following 
15!) and [16] 


t {e"" where 


numerical evaluation of the result 


(7 In the following discussion le 


Then 
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denotes imaginary part of. Furthermore 1 


n—1)p"W, cosn@? = 1 5 


where / 


ie 2 : 
> = (2 cos nd = Rel 1 
ad, (n + 1 2(1 


(p? + 1) cos 0 : cos 20 
p 6? 


2 
log (1 z | 
. (1 + 8%? 2d5p cos 8 


l l 
cos 8 log (1 + 6%? 25p cos 0) 
2(1 bp 


— log ( . ‘ dp sin 0 
log (1 —2 | sin 8 are sin - io 
(1 + 6%? — 25p cos 


Combining Equations [27 28], and [29] and letting 2 = dp . ' 
we obtain 1 pv nsinnO = . p* 
n =? op 


l 
i 


ee Pla 4 
ities 0 j : = p? + 1) sin #0 
0p 


6%p? 25p cos 6) 


25p cos 6 sin 8 log (1 + 6%p? 
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Gas Cooling of a Porous Heat Source 


By LEON GREEN, JR 


A limiting case of solid-fluid heat transfer is examined, 
in which a gas passes through a porous wall of high spe- 
cific surface with heat generation within the solid mate- 
rial. Dimensionless temperature profiles in the wall are 
presented in terms of the rate of heat generation, rate of 
flow, and thermal properties of the gas and solid. The 
pressure drop across the wall is approximated by using an 
average wall temperature and assuming isothermal condi- 
Temperature profiles, pressure drops, and pump- 
ing-power/power-output ratios are calculated the 
hypothetical case of a heated graphite wall cooled by 


tions. 
for 
helium. It is found that the thermal dependence of the 
gas viscosity produces a minimum in the pressure-drop 
versus flow-rate curve, and it appears that favorable 
pumping-power/power-output ratios can be obtained by 
the use of high pressures. The problera of temperature 
stability in a gas-cooled porous solid is pointed out and the 
need for experimental work emphasized. 
cooling technique for high-pressure, high-temperature 


Use of the sweat- 
ducts is suggested. 


NOMENCLATURE 


The following nomenclature is used in the paper 


= R/M,(in./deg R 
Btu /lb deg F 


gravitational acceleration, (in. /sec? 


b = gas constant 


specific heat of gas 
I 


weight flow rate of gas, (Ib/in.? sec) 


thermal conductivity of gas, (Btu/in. sec deg I 


apparent (bulk) thermal conductivity of solid, (Btu/in 
sec deg F 
thickness of wall, in 


molecular weight of gas, (Ib/lb mol 


fluid pressure, psi 


rate of heat generation in wall, ( Btu/in.* sec 


rate of heat removal from wall (Btu/in.? sec) 


— ft? I 
universal gas constant 10.73 x x ) 
in.? Ib mol deg R 
temperature, deg F or deg R 


= power required for pumping of fluid (Btu/in.* se« 
length variable 
viscous resistance coefficient of porous material (in 
= inertial resistance coefficient of porous material (in 
= dynamic viscosity of gas (Ib sec /in.* 


p = specific weight of gas, pef 
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INTRODUCTION 


The possibility of using porous construction to achieve efficient 
heat removal from a power-producing solid, although not a new 
Such 


idea, has recently been emphasized by Tsien (1).* con 


struction would reconcile high power density with small tem- 


solid The 


pose of the present paper is to study a limiting case of heat re- 


perature differences in the ind the coolant fluid. pur- 
moval from a porous wall by a gaseous coolant when the surface 
per unit volume of the porous structure is high enough that the 
temperature of the fluid at any point in its path may be considered 
Stud- 


indicate that this assump- 


to be substantially the same as that of the adjacent solid 
ies of sweat-cooled porous metals (2 
tion is not unreasonable 


TEMPERATURE DisTRIBUTION 


he assumptions involved in the temperature-distribution 


problem are summarized as follows 
1 Im the 
ire approximately equal 


2 The gas flow and heat flow are steady and one-dimensional 


porous wall, gas temperature and solid temperature 


Convective effects neglected by this assumption will alter the 
temperature distribution in the approaching gas at low flow rates 
but should not influence greatly the conditions prevailing within 
the wall 
$6 All heat 
t The ther 


! fic heat of the 


duction in the wall takes place in the solid 


onductivities of the solid and gas and the 


gas are constant 


4 gas at weight flow rate G through 
If the rate of 
i heat b 


flow ot 


The o 
a porous wall of thickness 


e-dimensional 


is shown in Fig. 1 


1 


heat generation per unit volume in the wall is q ilance on 


in element of volume in the wall (see Fig. 2) leads to the equation 


LT Ge, dl 


ty A da 


Wall 


A similar balance on a volume element in the approaching gas 


gives 


subject to the conditions that at 


solutions of ke 


Numbers i 


the paper 
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Fig. 3 presents curves of the dimensionless temperature dis- 
tribution in the wall as given by Equation (7) for different values 


of the parameter (G c,l k, As expected, it is seen that uni- 


formity of the wall temperature (desirable from the thermal 
stress viewpoint) is promoted by a high thermal conductivity in 








the solid. Manipulation of Equation [7] gives a relative tem- 

perature-uniformity index independent of the rate of power pro- 
duction in the form of the ratio 

DIMENSIONLESS TeMPEt 

Watt ror Dirrerent Vacut 


shown in 


Fig. 4. The absolute tempera 
however, increases with decre 


mum value of 


lim 
G-—-0O 


For calculation of the pressure drop a 
n the following) it is convenient to use 


erature 


Ne-DIMENSIK 


ERATING HEAT 


temperature 
1 Fig. 5. Sinee 
ELEMENT OF UNIT AREA ' 


1.0, the average wall temperature 


> Heat Balance on I EMENT IN WaLtI portior the flow rate, as would 
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tAGE TeMPERATURE INDEX VeRsus G C; 


FLow RESISTANCE 


Experiments with p etals (3 shown that 


the pres- 
isothermal! flow through 


d bya 


have 


rous n 
sure gradient in 
a moderately fine-grained px medium ca ex presse 


quadratic equation of 


R 


viscosity, and the 


where 0 gas constant, yp is nam 


ients a@ and 6 are length parameters char- 
+} f 


wteristic of the e pore material 


US 


illed the viscous a nertial 


terial, respective oefhicier with din 


L-?), character 


of “‘creepit o ‘ inertia for regligible The 


ensions 


US 
esistance terial in the re- 
gimme , 
vertial coeffi 
the addition 

1 within the 
A commonly 


da 


olds number 


flui 
correlation drop versus 
a Re 


equiva- 


pressure 


flow-ra media is formulated in yn- 


factor based upon the so- 


ler pi ulated by invok ilogy |t 


en the laws of Poiseuille, whict valid only for 


flows 


twe ire 


low-velocity would explain deviation 
y the 


are visualized as straight, parallel capil- 
showed that 


inalogy 


his 


from Darcy's law at higher flow rates b onset of turbulence 


within the which 
The 


length parameter is not 


pores 


laries porous-metal flow studies (3), however 


a single alwavs sufficient to correlate 


a large i f flow 


parallel] 


pressure-drop versus flow-rate curves over 


rates, which in general are not even 


curves 
suff 


cross one another in some cases of rently 


porous 
structures, 


The necessity of having two independent length parameters for 


of arbitrary 


medium 
of the di 


characterization of porous complexity 


the a 


ids number presentation 


considerably redu ivan mensionless friction- 
Although a 
i porous structure of 


by 


onvenmience In 


factor versus Reyno teyn- 
vumber and fr 
desc 


in dimensionless gro } 


ley 4 geom- 


13 
is achieved by s« 


fying ription may defined 


writing Equation 
no « ilculations 
With the introduction of heat genera 


the 


» doing 


tion within the solid material, as in present problem, the 


complications involved in dimensionless are so 


representation 
1ot justified, and greater simplicity is achieved 
13 


great that its use is 
by 


in the foregoing fort 


direct calculation of the pressure drops from Equation 
pressure-drop proble m would require 
the heated 
temperature trom 


d 


An exact solutior 
Eq i i 


wall, substituting 


integration of ross the thickness of 


point the pr 


per 


Equation [7 orrespo! ilue of the namic vis- 


ling v 


4 POROUS HEAT Sol 


RCI 


the « 


coefficients a and B are 


cosity of rhis procedure, of 
that the 
For thin 
ileulation using the 


inde pen lent « 
walls of high thermal conductivity, ho t 


sim} 


average wall temperature Equation 


is sufficiently ac 


urate 


Heuium Frow Turoven A Grapurre Wat 


As 


rous sol 


a specific example of gas cooling of a heat producing po- 


id, calculations of the anticipated performance of a porous 


gr aphite heater for helium gas have been made Cir iphite “a 


heater material for use with inert gases is of interest because it 


has been found to exhibit remarkable mechanical 
1,5 
ind the 


Grade 


properties 


at 


t 
« 


The hypothetical heater 


nperatures of ine scence 6 
material is assumed to 


0 


wall thickness is taken as 0.5 in 


National (¢ 


graphite, pertinent properties of whicl 


be arbon Company commercial porous 


Table 1 


sn intrepid 


are presented in 


rhe high-temperature viscosity of he 
the data of 
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| Akin | 
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the flow rate is a double-valued function o 


owing to the increase of the gas viscosit 


ind operation of the heater at the minimum 

y is indicated. The most extreme case 
! is investigated further in Fig. 9 which 
the pressure drop across the wall for inlet pressures of 5 


vivantage of 


tm, respectively, thus emphasizing the 
Because 
Table 1), 


r at high pressures of the relative 
I 


the porous structure (se¢ a small pressure 

ross the wall is imperative 

of merit commonly applied to heat exchangers is the 
the fluid to the 


pumping power 


the power required for the pumping o 


power removed as thermal energy Since the 


used in forcing a perfect gas through the wall is 


-removal rate 
Y 
« 


merit becomes 


pressure 


This ratio is ; function of flow 


APPLIED MECHANICS 


WEIGHT w RATE 
Porous Wat. Versus 
Heat GENERATION 


Fu S Presscre-Square GRADIENT IN 


Weicut Frow Rate ror Dirrerent Kates oF 
i portion indicates region where gas was calculated 


Dashe t 
extray 


4000 300¢ 

' ’ 
RADE 30 POROUS GRAPHITE 

120.5 in 

q* 10 Btu/in3 sec 

T, * 500 °R 


pee 


Pressure Drop Across W 4 
10 Brv /In.* Sec anv 


same case as in Fig. 9. It appears that in a high-pressure system 
the W/Q ratio may be reduced to a value low enough to permit 
useful power production from the heated gas. Gilliland (9) has 
shown analytically that the W/Q ratio for a gaseous heat-transfer 
system varies inversely as the square ol the system pressure, and 
this prediction is obeved by the data in Fig. 10, as is shown it 
i pressure ol 40 atm would be 


the W @ ratio tor the present case to 0.1 per 


Fig. 11, which indicates ths 


quired to reduce 
cent 


be emphasized t we f 1 results do not sh 
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low W’/Q ratios with 


’ coolar 


that 


with certaint) 


conclusively gaseous 


achieved merely by using porous construction 
In preliminary experiments by the author,’ however 
nitrogen-cooled NiLF, 


heat exchanger with a power density of 0.088 Btu /in 


using a 
tubular 
sec, W/Q 


were ob 


electrically heated, porous-wall 


depending upon the flow rate 
which «¢ 


averaged about 2 atm, it may be presumed that 


ratios from 0.2 to 2.0 


tained. Since the mean the svstem \hausted 


to the atmosphere 
a system pressure of about 20 atm could effect a W/Q ratio oi the 


pressure in 


order of 1 per cent It mav certainly be concluded that the idea 


® Unpublished work conducted at the Jet Propulsion Laborator 


of the California Institute of Technology, 1948 


GAS COOLING OF 


A POROUS HEAT SOURCE 


appears attractive enough to justify an experimental investiga 
tion of its possibilities. 


TEMPERATURE STABILITY 


The curves in Figs. 8, 9, 10, and 11 immediately suggest opera- 
tion of the porous heater at the peint of] minimum pressure 


drop. If for some reason nonoptimum operation is desired, a 


choice of two alternative positions (high-temperature, low-flow 
for each value of the 
The order of 


commencing the coolant circulation and the power generation 


or low-temperature, high-flow) is available 


pressure difference established across the wall 


then determines which position is attained. 
Instability of the wall temperature, however, n 
the 


hot spot in the wall exist owing to some inhomogeneity of the 


arise from 


temperature dependence of the gas viscosity Should « 


porous material, for instance, the increased viscosity would tend 
to cut down the flow at the point, thus increasing the local 
The 


conceivably could increase in this fashion until failure 


tem 
perature and aggravating the situation wall temperature 
occurred 
In the foregoing experiments by the author, no such instability 
was noted, but in this case the power density was low, operation 


was well into the high-flow, low-temperature region, and a con- 


siderable amount of heat was conducted out through the elec- 


trodes. Because of the difficulty in analyzing the temperature- 
appears that an 


stability problem in the porous heat source, it 


experimental approach would be more fruitful. It may be pre- 


sumed, however, that stability will be promoted by high thermal 
conductivity in the porous solid, high heat capacity in the gas, 


and high flow rate. Stability considerations indeed may require 
operation on the high-flow side of the optimum point, where the 
flow 


of Equation [13 


resistance is dominated by the quadratic nonviscous term 


Experimental work will be required to resolve 


this problem 
CONCLUSION 


The foregoing discussion has explored the problem of heat re 


moval from a porous heat source by a gaseous coolant. Con 


sideration of the specific case of a helium-cooled graphite wall 
showed that efficient operation from the point of view of pump 
ing power would require the use of a system under a pressure of 
at least several atmosphe res The use of such high pressures a“ 
the extreme t« mperature desired poses & serious pr »blem of duct 
ing or containment of the system, since the creep of most mats 
In this 
Duwez ( 10 


rials at such temperatures prohibits their use 
the 
for use in the cooling of rocket motors ippears to o Ter interesting 


conner 
thon sweat-cooling”’ technique proposed by 


section of a high-pressure, high-temperature 


Analysis (11, 12 


have shown that effective 


possibilities, A 
sweat-cooled duct is idealized in Fig. 12 
13, 14 
whieved by a relatively 


and 


experiment ooling of the wall 


ean be small flow of coolant, as is indi- 


cated in Fig. 13 


Since the effectiveness of the porous-wall cooling method (in 


cases with a temperature gradient along the wall) has been found 


to be greatest for wall materials of low therma mductivity 


11), the use of ungraphitized carbon is indicated for the present 
1000 F, for 
100 | rate of onl 


0.2 per cent of the main stream rate would reduce the inner 


proposed applic ition For a stream temperatur rf 


instance, a coolant flow (from a reservoir at 
w“ all 
surface temperature to about 2450 F The pres¢ e of the po- 
rous shield probably would eliminate the need for | gh-tempera- 
In this scheme the 


duct would be subjec t to only a small « OMpressive stress re quired 


ture materials in the outer casing. hot inner 


by the pressure difference driving the coolant, thus minimizing 
the creep problem, and the high pressure would be contained by 
the cool outer casing. Sweat-cooled ducts as applied to a gas tur 


been 


bine (the coolant air being taken from the « ompressor ) have 
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high enough to exploit the compactness and high permissible 


power density peculiar to a porous heat source 
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The Marcus Method Applied to Solution of 


Uniformly Loaded and Clamped Rectangular 


Plate Subjected to 


By C. C. CHANG! AND H 


This paper presents a solution to the small-deflection 
problem of the uniformly loaded rectangular plate with 
all edges clamped, when the latter are subjected to uni- 
form tensile or compressive forces in the plane of the plate. 
Since an exact solution to the problem is hardly possible, 
the Marcus method was used for the approximate solution 
of the biharmonic equation, subject to the conditions of 
zero slope and deflection at the boundary. Deflections and 
moments are presented for a square plate with four par- 
ticular values of tensile end load, and for the rectangular 
plate with seven length-breadth ratios and with one value 


of compressive end load. 


INTRODUCTION 
HE problem of the uniformly loaded rectangular plate sub- 
I I 


jee ted to forces in its plane is frequently encountered in the 


If the 


edges of the plate may be assumed free from bending moment, an 


design of such structures as a ship’s bottom plating 


exact solution to the small-deflection problem is easily obtained in 
3 The re- 


sults of such calculations and when a pair of opposite edges are 


the form of Fourier series using the Levy method (1 


subjected to uniform tensile or compressive forces have already 
been presented 2 fesults for the case in which al 


edges are 
subjected to uniform tensile forces have also been given (3). 
We now consider the problem when all edges are « lamped and 


subjected to uniform tensile or compressive forces. With the 


O*w J*w 


usual nomenclature, the differential equation to be solved is 
O*w q 
oy* D 


P (& “ ) 
D \ oz? oy? 


a positive value for the load P per unit length denotes ten- 


Vw ——- 
or‘ Ir*Oy* 


wher - 
sion and vice versa. Taking the co-ordinate origin at the center of 


the plate. the boundary conditions are 


Ow 

=(Qatr= 
or 
Ow 


Oy 


=Qaty = 
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where a and 6 denote the lengths of the sides of the plate. Having 


determined the function w, the bending and twisting moments per 


unit length of plate are 


Ow 


Dol 
Oroy 


Since an exact solution to the problem is hardly possible, use is 
made of the Marcus method 


This method has recently been used by E. 


4) for obtaining an approximate solu- 
tion to Equation [1] 
G. Odley (5) in his analysis of clamped rectangular plates sub- 
It is apparent from Odley’s 
paper that the Marcus method gives results which are sufficiently 
accurate for most practical purposes. The main of the 
method is that the meeting of specified boundary conditions pre- 
sents little difficulty. 

We assume that a solution to Equation [1] 
the form 


jected to hydrostatic pressure 


value 


may be written in 


w= KXY {4} 
where K is a parameter and X and ¥ are, respectively, functions 
of x and y only (Adopting the nondimensional co-ordinates 


= 2 . 
$ = 2r/a; 7 


and substituting Equation [4] in Equation [1], the intensity q, is 


found to be 

exp |< 2% 2 d*X d*¥ 
) , t 
, a‘ dé a*h* dg? 
Y d*x 


—4PK + 


a’ dg? 


X dy 


=” b dn‘ 


dn? 
X d*} 
b? dn? 
If the function X of z were known, a determination of the func- 
tion Y of y could be made as follows: 

From Equ ation [5]. the load Pp, on a strip of the pl ate of unit 
width is 


/ ve 1GKD | 
q ir = ( 
a’ b4 


Pat , \ ay 
4Db? “*) dn? 


1 
d‘X 
[ ~ dt, B, = [ 
o «dé /0 


However, the actual load on this strip is ga and, after writing the 


d*‘} 
* dnt 


obs 


aX 
dé? 


where 


A, = : € 


following expression for K 





ga‘ 
4(4DA, Pa®B, 


Pa* | d*?Y b‘ | A, 
+ 

8D | dn? a‘ LC, 

Q) 1) = 0 [7 


K = 
we obtain 
a‘y 
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There are two forms of solution to this equation which are re- 
quired for our problems. For 


B, Pa? |* A; Pa? B, 
C. 8D "8G, 4D C, 


(A sin $,7 sinh ¥,9 + B sin $,n cosh 9,7 


the solution is 


+ C cos $7 sinh ¥,n + D cos d, n cosh ¥,7 


“= ) 
8D 


Pa* B, 
iD C 
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Pa? |? A 
8D as 
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If the function Y of y was known instead of X, similar calculations 


while for 


B, 
C, 


e solution is 


Pas 


where 
4D 


4D 


would yield the load p, on a strip of the plate of unit width and 
perpendicular to the former strip 


16KD | 


where 


Writing 


gb! 
144DA, — Pb*B, 


Pa? 2: dg? -_ 


Pa*® B, 
; the @ and y are reduced to the forms 


we obtain 
d‘*X 49 a* B, = 
dg b? | ¢ 9 


As before, there are two forms of solution of this equation which 
are required for our problems. For 

Pb* B, 

4D C, 
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C, 8D | ~ C, 


the solution is 
(A” sin $,£ sinh Y,& + B” sin yf cosh y,£ 


+ C" cos ¢,£ sinh ¥,& + D” cos ¢,& cosh ¥,£) 


y* . 
Pb? B, (3s Pb? ) 
iD C, Cc, 8D 

Pb? )| 


Pb? | d*X a‘*| A, 
8D | dé bc, 


¥-—-1) =0 11] 


X=1 


where 


Pb? B, (> 
iD C, 


Ww hi e for 
v 
Cc 
. 


—(A'" cos ¢,'E + B’ 


+ 


the solution is 
X =1 
where 


Ph? 


iD Cc, 


Pb? B 
iD Cc, 


and 


We- tlt 
=\i-q> 
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maximum moment o rs at the center of = | + 0.020244 cos 2.21337 1.00587 cos 0.189427 
ind a/l 5 or 2, the maximum 
A eee: PERRET ee lo enahteiealines ipproximation to the deflection surface is then 
Kil m mome oves t« rd the shorter sides qd 
num moment moves toward r sid 1 q 
24 16D 
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Measurement of Recovery Factors and 


Friction Coefficients 
Flow of Air ina Tube 


for Supersonic 


2—Results Based on a Two-Dimensional Flow 
Model for Entrance Region 


By J. KAYE,' T. Y. TOONG, 

The first part of a program to obtain reliable data on the 
rate of heat transfer to air moving at supersonic speeds in 
a tube has been devoted to measurements made on adia- 
batic supersonic flow of airinatube. The details of these 
measurements have been described in a previous paper. 
The calculated quantities such as the local apparent fric- 
tion coefficient, recovery factor, Mach number, and so 
forth, were obtained from the simple one-dimensional! flow 
model for which the properties of the stream are uniform 
at any section, and boundary-layer effects are ignored. 
The analysis of some of the same data given in the previous 
paper is undertaken here with the aid of a simplified two- 
dimensional flow model. The supersonic flow in the tube 
is divided into a supersonic core of variable mass with the 
fluid remaining in the core undergoing a reversible adia- 
batic change of state, and a laminar boundary layer of 
variable mass. The compressible laminar boundary layer 
increases in thickness in the direction of flow, and then 
undergoes a transition to a turbulent boundary layer. 
The two-dimensiona! flow model is limited here to the 
region where a laminar boundary layer appears to be pres- 
ent in the entrance region of the tube. The results of the 
analysis based on the two-dimensional flow model indicate 
that where the flow in the tube boundary layer appears to 
be laminar, the measured pressures and temperatures in 
the tube for adiabatic supersonic flow of air could have 
been predicted, with sufficient accuracy for engineering 
problems, from measured data for supersonic flow of air 
over a flat plate with a laminar boundary layer, and with 


zero pressure gradient 


ano R. H. SHOULBERG, 
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without heat transfer, have been made for many years 


at each cross section has been found 


developed flow is attained 


neering problems where fully 


au adiabatic wall conditions 


i boundary layer 
r isentropic core 
station numbers 
upstream stagnation conditions 
atmospheric conditions 


sentropie conditior 


INTRODUCTION 


and 
Most of 


Measurements con with fluid flow in ducts, with 


erned 


these measurements have been made for subsonic flows in pipes 


\ simple one-dimensional flow model based on uniform properties 


udequate for calculation and 


design purposes, provided the entrance region is neglected or fully 


Fortunately, there are many engi- 


leve lope d flow is obtained in a pipe 
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or where the entran ! the total length of 


flow 

Several types of engineering designs exist where the entrance- 
region effects have a large influence on the performance of the 
equipment, such as in ducts in wind tunnels, heat exchangers 
ducts in gas-turbine power plants, and other applications where 
pipes having a relatively small length-diameter ratio are used 
For short pipes, the simple one-dimensional flow model has beet 
replaced, in most mathematical treatments of the flow in the er 
trance region, by two-dimensional flow models or axially sym 
metrical flow models. These more compli ated models have not 
come into general use for engineering designs and calculations, 
but are used mainly as a guide in estimation of deviations from the 
results based on the one-dimensional flow model 

Phe simple one-dimensional flow model has been used with some 
degree of success in computation of mean and local friction co- 
efficients for supersonic flow in a tube It has been shown in a 
previous paper (1),* however, that this simple flow mode 
fully adequate in the interpretation of the calculated quantities 
ind that a more detailed analysis is needed for full comprehension 
of the phenomena which occur in the entrance length of the tube 

s result is not surprising if it is realized that it is usually not 
possible to maintain supersonic flow in a length of pipe corre- 
sponding to that required for fully developed flow In fact, in Fic 
many measurements given in the previous paper (1), it has been 
$2 The following assumptions were made 


ound that the “entrance” region extends over the entire 


diameters of length of supersonic flow simplified two-dimensional flow model 


4 simplified two-dimensional analysis of supersonic flow in a 
i ia ‘ I , — = 1 A laminar boundary layer exists over the entrance re 


tube has been examined in detail. The objectives in this analysis , 
the tube 


have been to analyze quantitatively the experimental data in 2 The air remaining in the re underex n 
4 , 1 core 1dergoes a 
terms of th nin boundary laver at tube entrance, to « are . . 
ern he lami art indar iver at tube entrance, to ompare change of state from the stagnation pressure to the pre 
quantitatively the calculated friction coefficients and recovery fac- : 
given section, 


tors with those for supersonic flow over a flat pl ite, and to predict 3 Thestati pr ure is uniforn ' } I 
, » static pressure is uniform across each section of 


ritit ; tear ; ' r ] . » 
quantitative he behavior of supersonic flow sca Mog re- 4 Air is a perfect gas with a cons 
ri t = pis wi i I dary ver c 
gion of variou iminar bounda iver and an specific heats (k = 1.40) over 
isentropic core are present s also evident that the simplified ‘ 

. onsideration 

wo-dimensio flow model n be used t \ » superse ( - oo. 
two-dimensiona w mode ised to analyze the supersoni ag ey Sar er ee eee 


flow ith | transfer in of the laminar 


velocity which is further ipproximate 
the boundar iver equa } 
rOR Two-DIMENSIONAL Flow Mont in Fig. 1 
’ ' 6 The 
piihed two dimensional flow model discussed here 
the round tube is arbitrarily divided into two parts a- 
shown in F 1. The first part corresponds to the core which al- It should be noted that assumption 6 
the Prandtl number is unity Sinee the 


most fills the tube at the entrance section and decreases in size 
the direction of f The second part corresponds to the less than unity a temperature discon 


= usually of small size compare the temperature gradient (d7 


laminar boundary layer which 
with the core at the entrance section and grows in thickness in the tion 6 is not used, however, 
lirection of V the e diminishes in size, mass is trans merely ret an approximat« 
laver. The arbitrary thick vel n boundar 
rresponds approximate nalysis is similar 
in the boundary la 


» understandi 
paragraph are consist » be use 
laver. His basic assump 

ould be divided into 
ire assumed to 
gion where viscous forces 
in the tube it is assumed here that 
rees on the fluid remaining in the core 


negligible and that the viscous forces are con- tinuits 


centrated mainly in the boundary layer. Hence it is possible to 
of the fluid remaining i Continuity: 


determine simply an 


the core by using isentropic : 
Geometry: 


Equation of 





KAYE, TOONG, SHOULBERG 


p-RT 


Equation of state: p 


Definition: M.? V.2/gkRT 


Assumption: V, 


The following relations hold between the upstream stagnation 


condition and any tube section shown in Fig. 1 
Isentropix 
Isentropic 


knergy: we 


Assumption: 


From Equations follows that 


9 


a 
«J 


Note that 7’, is, in reality, defined by Equation [11] 
Application of the momentum equation for steady flow to a 
control volume which encloses the fluid between two adjacent 
stations of the tube yields 
+ w,V 
tr D(L ° 


g 


+ wV, ] 
q 


= 


and the wal! shearing stress is given by 


<9 


friction coefficient f;.:/, is evaluated in 


terms of the fluid properties midway between the two adjacen 


where the true local 


stations 


The additional relations used in this analysis are as follows 


Definition: c, 


Definition: ( - T.) ° 17 


Vethod of 
straightforward 


The method of 
using the measured values of pressure and 
adiabatic wall temperature along the tube. 


Computation 


computation is 
Since a sample calcu- 
lation is given in the Appendix, it will suffice here merely to indi- 
cate the method of computation 

Starting with the measured stagnation pressure and tempera- 
naining in the core can be 


ture, the states of fluid re isentropic 


evaluated at each station from the measured station 


The recovery factor 


pressure.,® 
Equation [17] can be evaluated imme- 


ulations of the details of the boundary 


liately without further 
layer. This method of calculation of recovery factors 
1949 } - 


tation of the first report or 


was used 
previously (June some of the present authors in presen- 


this work (3 


, and has hee n used later 
The details of the laminar boundary layer are 


n Eq 


is obtains 


by other authors 
then obtained fror tions [1] to [13 The true local fri 


tion coefficient 1 from Equations [14] and {15}. 


RESULTS 


In a previous paper the measurements for adiabatic super- 
sonic flow of air made in forty runs with three test combinations 
are presented in detail. It is evident that the use of the simplified 
two-dimensional flow model to convert ail the data in these runs 


would involve a considerable amount of calculation. Therefore 
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only six runs made with test combination B have been selected 
for computations made with the various two-dimensional flow 
models. These six runs correspond to those selected from test 
comparisons made with the earlier test 
This 


comparison of calculations 


combination B for the 
combinations given in the previous paper (1 election of 
the same six runs permits a direct 
based on the one-dimensional flow model with those based on the 
simplified two-dimensional! flow mode 

One of the objectives of the two-dimensional mode! is to com- 
pare the results quantitatively with the well-known behavior of 


the development of a laminar boundary layer on a flat plate with 


zero pressure gradient in the direction of flow. For this reason the 


characteristics of the boundary layer on a flat plate are used as a 
basis of comparison on the various charts which will be discussed 
It should be noted that this basis of comparison 1s a temporary 
expedient to be used unti] an exact solution for adiabatic 
supersonic flow in the entrance region of a round tube is availa- 
ble 

The six charts 2, 3, 4, 5, 6, and 7, 


4 cover the range from 
laminar flow to turbulent flow 


Figs. 
n the boundary layer of the tube 
None of the data obtained in the original measurements is shown 
on these charts since all] such data are given in the previous paper 

1). Thecaleulated quantities based on one flow model are plotted 
on each chart without reference to the quantities calculated from 
Thus, in the curve of Mach number versus 


the 


the other flow model 


length Reynolds number, the values of Mach number and 


Reynolds number are taken from either the one-dimensional flow 
model or from the two-dimensional flow model 
Each of the 


Mach number, 


six charts, Figs. 2, 3, 4, 5, 6, and 7, presents the 


yer thick- 
and local friction coefficient For the 


mean stream temperature, boundary 


ness, recovery tactor two- 


the Mach number and mean stream tem- 
and the 


dimensional flow model 


perature reter to the isentrop ore recovery tactor is 


given by [Equation 


Cod Since the phrase s two-dimensional flow model” and 


one-dimensional flow mode occur frequently in the discussion 


forms 2 


of results they will be referred to in the abbreviated 
DFM and 1-DFM, respectively 

Boundary Laye Figs. 2, 3, and 4 compare the calcu- 
lations based on the 2-DFM with those based on the 1-DFM for 
those values of the 


evidence indicates a laminar boundary layer exists for some dis- 


Lamina 


the diameter Reynolds number where all 
tance downstream from the tube ent 
In Fig. 2, for the smallest \ the 
Mach number based on the 1-DFM decreases from about 2.2 to 
1.4 from the tube entrance to exit Mach number 
based on the 2-DFM decreases from about 2.6 to 2.2 
pattern is also found in Figs. 3 and 4 for larger va! 
Thus the 2-DFM 


velocity in the supersonic core in the 


rance, 
ilue of the Reynolds number 

whereas the 
This same 
the 


a larger 


teynolds number. results in both 
nearly con 
1-DFM 
mean stream 
for the 2 


Lube and a more 


stant value of this velocity than that based on the 


A similar result is found in the calculation of the 
for the 1-DFM 


temperature 
DFM. The 


temperature base 


ind the core temperature 


ratio of mean stream 


don the 1-DFM increases 


exit by about 45 per cent in Figs. 2 


« , ana ere 


temperature to 


ignation 
from tube entrance to 
is the ratio of 
the 2-DFM 
Thus the 


constant 


ore temperature to stagnation temperature i " 
ent in Figs. 2,3 vi 4 
2-DFM is more nearl: 
based on the 1-DFM 


ing the 2-DF™M is the existence of 


increases | mily about 20 per 


core temperature based on the 


than the mean stream temperature 


The basic hy pothesis under), 
a laminar boundary layer in the 


This hypothesis can be tested by 


entrance region of the tube 


omputation of the thickness 
The ratio of this 
and 4 This 


ratio increases from about 10 to about 45 per cent from the first 


of the laminar boundary layer in the tube 


thickness to the tube radius is shown in Figs. 2, 3, 
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Howarth's results is very good in 
measurement L, which should 
zero boundary-layer thickness, | 
at the end of the curved cor 
boundary-laver thickness Is sm 
zero boundary-layer thickness 
used in Howarth 


result 


Apparently, in Figs, 2, 3, and 4 


tion corresponds to an entrance 
condition | is ty f supersonic flow since the laminar end ¢ 
boundar r does not fill the tube at exit and a fully develope probab iused by the transitior 
tlo rn is not attained houndary layer, corresponds t 
In addition to computation o it i ol he lamina tities at about the same ler 


boundar iver in the tube on 1 ) he 2-DFM, th ! rig measureme 


thickness 1 be computed from t theoretical predictions for previous pap 


straight line wher 


it plate with zero pressure gradient in the directior 


incompressible flow ve results of Blasius (4) pre average deviation o 
the those at the first ar 


» equation 
however, that 


lO4L DD) vy (Re " we | DFM v 


kness Dy 


Howarth (5) for a compressible boundat I ‘ 
mnhiirmed 
2-DFM 


the equation 


1O.4(L Dil + 0.19870) 1M Vv (Re 19 The recovery factors based o1 

could be represented by 
Che thickness of the laminar boundary layer calculated from iverage deviation of less than 0.4 per 
Howarth’s work is shown in Figs. 2, 3, and 4. The agreement the first and last stations, and by a 


between the thickness calculated from the 2-DFM and from irger value of the average deviation 
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lon the 2-DFM exhibit tly the 
that shown by the original measured values of 
This result implies that the re- 
a tube based on the 2-DFM 


with those on a flat plate with a 


covery tactors Dast simost Xa 


same variation as 


the adiabatic wall temperatur 


overy factors for supersonic flow in 


should be direct] mparabk 
free stream condition identical to the core condition in the 


The few measurements available on recovery factors 


speed with a laminar boundar' 


Rubesin, and Tendeland (¢ 


around bodies at supersoni 
This 
about 0.86 for a Mact 
body of 


gradient. Ir 


are summarized 
summary indicates the recovery factor is 
revolution 


the r 


number of about 2 I ibout a paraboli 


ressur 


even in the 


prese 
measured recovery factor for a Mach number o 
laver is OSS] 


fle late with 
10 
2-DFM, 


ind OLSS3, as 


length Reynol 
covery tactors 
M ach number ‘ 75 
Figs. 2, 3, and 4, 1 tively The 
is 0.883 + 0.005 ment should be 


OSU] 
average of these three 
onsidered a coin 


dence until the ressure gradient on both the tube 


plat 
number o1 t overy 


flow and the igated in greater detail. The 

effect of M i h 

DFM appears to be sma 
The analysis of I 

over a flat plate with a compressible 

dicts that the recov y factor should vary as the 

Prandt! number or a 


the 


factors based on the 


mons and Brainerd (7, 8) for superso 


laminar boundary lay 
square root of the 
ompressible boundary layer, some loubt 


exists as to mperature to be used in evaluation of the 


Prandtl number 2, 3, and 4 show a curve corresponding to 
the Prand+t! number evaluated at the core 
perature based on the 2-DFM 
the 2-DFM are in excellent agreement with the theoretical pre- 


the square root of tem- 


The recovery factors based on 


SUPERSONK 


FLOW OF AIR IN A TUBI 189 


» 


n Figs. 2 and 4 


diction und are about 2 per cent larger than the 


predic ted values in Fig. 3 
ution for 


The original measurements of static pressure distr 


flow of air in the tube as well as t local 
coefficients based on the 1-DFM are contained 
l It 


for the 


idiabatic supersoniv 


apparent [riction 
in the 
DFM 


use of the 


there " l- 
The 


Iric- 


previous paper was indicated 


Was Inadequate entrance region in the tube 


2-DFM, ho 
coefficient since } 


wever, should result in a 1¢ local 


tio snges in velocity profile are taken into 


sccount approximate! Henee the true local friction coefficients 


should be ir eement with those for flow ov 


ilues b 
In Figs. 2, 3 the true 
2-DFM are compared 
1-DFM 


ver a flat plate with zero pre 


witl arent tion 
with those incompressible 


Phe va 


curve marked “‘flat plate 


ients based on the and 


ri ior 


sure gradient 
incompressible flow are given by the 

Phe effect of 
i M acl 


about 5 


sminar ompressibilit laminary boundary 


number ¢ 5 ‘ true friction co 


Phe contrast between the 


two flow models is immediately apparent in gs. 2 and 3 where 


ient curve for the 2-DFM is greatly 


the slope of the friction-coeff 


slope of the curve based on the 1 
of the 2-DF M is now 


} as that for the flat plate 


reduced compared with the 
DEM. In addition, the slope 


urve for the 


en sume in Figs. 2 and 
‘ the 
flat- 


the 


wreement shown in Fig. 4 betwee 
on the 2-DFM and the 
could be considered accidental if it 
Fig 4 There 


yveryv factors based oa the 


The extraordinar 
true | friction coefhicient rane 
plate Values were not lor 
pertect 


2-DFM 
und ver 


evidence shown it 


wgreement between t we Ted and 


those for a flat plate with a laminar boundary layer v good 


values of the corresponding boundary- 


chart ¢« 


between the 
his 


e that co nparison with the flat plate value t 


igreement 

thicknesses ontains the 1 t ubstantial 
temporary 
interpretation ol the 


phe- 


it means of 


na which occur in adiabatic supersonic flow of a na tube 


entran egion. Furthermore 


this cl indicates 
2_DFM so as to 


mg the 


the 


of supersonic 


flow of 


or rina tube us 


a flat plate 


coethi 


supersoni air o1 


iImetion rent 


the true 


because of the 


boundary pom 


i true ir 


al 
0 


siderable error can be 


extrapo 


ikculations Dase 


he 1-DFM for the 


tub 
youndar 
ir boundar 
base n the 2-DFM wou 
with th 


values in tl 


rION OF Surersonic Frow Bevavior 


good agreement between the values of the atninar- 


boundary-layer thickness, recovery factor, and local true fricton 


coefficient, based on the 2-DFM, and the corresponding values for 
1 flat plate indicates that this 


a laminar boundary layer or 


simplified and approximate 2 DI M ma be Use* 


1 in reverse to 
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rhis method of prediction of distributior pressure and tem- 
perature for adiabatic supersonic flo yul sufficient], 
vccurate for Many engineering problems s up numerous 
ther possibilities, such as the prediction of pressure and 


prea 


tem- 
perature for supersonic flow in a tube witt at transier Work is 


in progress on this phase of the program 
lata on supersonic flow in a tube 


CONCLUSIONS 


Since exact solutions are not available for the adiabatic super 
sonic flow of air in a round tube with a laminar boundary layer 
leveloping from tube entrance, an approximate and simplified 
two-dimensional flow model is proposed for engineering pur- 
poses The results based on this two-dimensional flow model 
have been compared with the corresponding quantities for 
adiabatic supersonic flow of air over a flat plate with zero pressure 
gradient. The agreement between the corresponding quantities 
is good for all cases wher« r boundar iver appears to 


exist in the tube 





predict hav ‘ personi lo ot in a tube, and 


possibly for other fluids and for other pass . Hence it appears 





possible to predict lor engineering designs the distribution of pres- 

sure and adiabs lt erature along the length of tube from 

data based on a laminar be lary laver for supersoni flow on a 
flat pl 

The an ) ised ry t of the distribution of 

pressure and i t is identical with that 

stagnation st > and 

to use two new rela- 

pressures and temperatures 

ind either of Equations [19 


ilong a flat plate could be 


0.664. \ Re 


There would then be a sufficient numb independent relations 
to specify completely the distribution of pressure and temperature 
for the supersonic flow in the tube where a laminar boundary 
layer and an isentropic core were assumed to exist. The prediction 
method should be used, however, only up to a length Reynolds 
number of about 10*, corresponding to the start of the transition 
to a turbulent boundary layer 
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Fie. 7 


The results, based I he two-dimensional flow model, com- 


pared in detail with those based on the one-dimensional flow 
model customar ed ia much better understanding of the 
phenomena oO rrin ns rsor oO ir in a round tube 


Some of the 1 nsional flow model for 
the entrance reg e¢ nated by t the two-dimensional 
flow model 

The Mach 1 rs b | 


are slightly lar 


on the two-dimensional flow model 


1 more constant along the tube length 


than those » one-dimensional flow model 


Phe thickness « 


nsional flow 1 | is in goo eement with that based on 


base 


boundary la 


er based on the two 
aime 
Howarth’s results t 
zero pressure gradi 

Where 


laminar, the re 


the flow in the t boundary la appears to be 
used on the two-dimensional flow 
and with 


model agree withir er cent both with the measured 


the theoret 


supersonk flow oy 


a lar al inda ‘ Where the flow 


the r ove! 


er a flat 
plate wit! 
boundar ver vent t urbulent 


based on the . “i Vv mn ] r pr »babl 


cidence, with supersonic flow 
flat plate wit! 


Where 


laminar, the 


sundary iayer 


boundary layer appears 


ients based on the two-dimen- 
agreement with the experimental 
, 


layer on a flat plate with zero 
layer 


sional flow model are 
alues for a compressible be lar) 
Where 
. ae 


appears to be turbulent, the loca 


pressure gradient the flow in the tube boundary 


friction coefficients based on the 
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two-dimensional flow model are, roughly, half of the experimental 
values for a turbulent boundary layer on a flat plate 

The success of the two-dimensional flow model indicates that it 
may be used in reverse to predict the pressures and temperatures 
flow of air in the entrance region of a 


laminar boundary 


for adiabati 
tube 


supers« mk 


with a laver In other words, the 


data for adiabatic flow of air over : plate with a 


supersonk 


iminar boundary layer may be used to predict the behavior of 
adiabatic supersonic flow in a tube to a good approximation for 
engineering problems. In addition, it is probable that a similar 
method of prediction may be used advar tageously lor supersonk 


flow of air in a tube with heat transfer 
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Appendix 


' 
ts based on 
mensional fl th the adi 

section rete “i tol prev paper (1) as test 
combination B s pole ulation » included rhe 
original measurements t x runs, t ensions of the 
apparatus, and the other 


l A 


dimensional flow model is given { 


summary 


Sample Calculation. The original data 


taken fron A ppendix A of the previous pape method of 


computation of the nozzle discharge coefficier al with 
that used in the sample calculation of run No in the previous 


paper (1). The values are summarized ¢ 


Gt," 28.45 Ibm /( 
Re,,).” 5.684 « 10° 
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= 0.960 
w,* = 0.01119 Ibm/se 
w = 0.01074 lbm/sec 
A = 0.1978 in.? = 0.001373 ft? 


30 of the “Gas Tables’’ (2), corresponding to the 


of 0.04731 for the first station, it is found that 


From Tabk 
Ps Po 
M = 2.64 
T./T., = 0.4182 

T, = (7,/T.JT« = 


value of 


0.4182( 565.3 236.4 F abs 


rhe recovery factor is found next by 


0.9394 0.4182 
0.806 


1 0.4182 


8.834 < 0.04731 & 144 


53.342 236.4 


= 0.004772 Ibm, ft 


gkRT, 2.637 V 32.174 1.400 X 53.342 * 236.4 


= 1988 fps 
4 486 lbm/(see ft 


OF 0.004772 LOSS 


the boundary layer 
2 = 993.9 


V,2/(2gc, 565.3 
993.9 2 
483.1 F abs 
2 X 32.174 & 778.16 & 0.2399 


0.04731 & 144 
483.1 


8.834 xX : ; 
- = ().002336 Ibm /ft 


53.342 X 


0.002336 993.9 


From Equations [1], 
0.01074 
2.321 < 0.001373 
2 X 0.004772 


0.7676 

0.002336 
0.7676 0.001373 
0.001373 0.001054 


0.001054 = 0.010000 Ibm/sex 
= 2.321 « 0.000319 


0.001054 sq ft 
0.000319 sq ft 
YASH & 
0.000740 Ibm /se« 
24/0.5018) VY 0.001054/r = 
0.1239 


r= 1 


The shear stress at the wall and local true friction coefficient 


are obtained from the momentum theorem 


0.1978(0.04730 
0.000740 


0.05993) + l 


32.174 
32.174)(0.009235 


wr X 0.5018 x 2 
( 12 X 12 


= {8.834 x 


(0.010000 1988 + 93.9 ] 


x 1937 + 0.001505 » 


= 0.7889 Ibf/ft? 


x 0.7889 


0.004772( 1988)? 


0.00253 
+ 0.005650( 1937 )? 


G.D(L/D) x 0.5018 
Le ; 12(0.619 > 


9.486 1.614 


(Re,), = 10-8 








Lateral Vibrations as Related to Structural 


Stability 


By HAROLD LURIE, 


The apparently different physical problems of lateral 
vibration and elastic stability are limiting cases of a single 
phenomenon, the more general expression being the mode 
of vibration with end thrust. It is shown that the square 
of the frequency of lateral vibration is approximately lin- 
early related to the end load. The linear relationship is 
exact if the mode of free vibrations is identical to the buck- 
ling mode. In all cases, the load corresponding to zero fre- 
quency is the critical buckling load. Experimental tests 
were conducted on elastically restrained columns in the 
form of rigid rectangular frames. It was found that the 
relationship bet ween the square of the frequency and the 
load is practically linear, and that the extrapolated load 
corresponding to zero frequency coincides with the buck- 
ling load. Rigid-joint trusses also were tested. In some 
cases the relationship between the square of the frequency 
and the load deviates considerably from linearity, but the 
approximate buckling load can be predicted. In the case 
of thin flat plates, tests showed that the buckling load can- 
not be predicted from vibration tests in practice. It is 
indicated that this is probably because the linear plate 
equations are not valid owing to initial curvatures in the 


plate. 
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= transverse load per unit area 


) = mass per unit area 


INTRODUCE TION 
In addition to the mathematical similarity between the simple 
some striking 
A 


onsidering forced vibrations and 


problem of free vibrations and of elastic stability, 


when certain modifications are introduced. 


led by 


of an initially curved member 


inalogies occur 


wd example is provi 
stability 


leflection function of a uniform beam is expanded in a series of 


Suppose the static- 


tior Phen it is 


modes of vib 


the 


unctions corresponding to the i 


deflec- 
hb 


vell known that, under a harmonic load shape of the 


obtained 


on curve for a forced vibration with fre quency @ is 


multiplying every term of the statie-deflection function by the 


i/[r—(2) 


Lu 


resonance factor’ 


d 


deflection 


in be l 
imn is expanded in a Fourier series, the shap 
under P is obtained by 


the tion b 


ilso shown * that if the shape of at irve 


an end thrust i m ing every 


efle 


curve 


iheation 


tion func 
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inloaded wn 


1 the 


It is this similarity between the onance mag 


might lead one to suspect some 
The 

a beam which is pinned at each end 

the 


nification factor that 
I 


readily in the case of 


n ¢ be rroborate quite 
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out any external forces being applied to stem equation 


governing the free harmonic vibrations 
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lhe members of a structure are otter 
1 determination of their natural 
in order to prevent resonance f suct w of variation o 


natural frequency with end thrust car ablished, it would 


lead to useful practical applications i form of nondestructive 


tests for the determination of the st 
GENERAL Eg 


pression P is to decrease the Before examining any specif 


vibration. In addition, the variation of the sible to study some general expressions 


frequency with P is linear; and finally, the fre tions of continuous systems such as be 
» when P approaches P. which is the buckling ntial equation for the vibratio t 
the general form 
irises as to t g ] f this varia- 
1 behavior spected in structures 
supported unifort cal intuitior vhere the subscript ¢ tal 
ne to believe that if 1a | compression does de- variables. Liu] is a lu 
ny structure mly the unknown defle« 
respect to the space variables 
the beam 
tion of 
rge number 
s, rings, and 


ious boundar 


In the foregoing general equa 


{) is proportional to the densi 


o find the In the problem under consicde 


us made t is of primary interest Mathem 


nd stability tion is given by a 


lifferential equation 


together with the initia 
Separating variables 
th the eigenv 

urit by 
ind proposed a method for determining 
requency ™ wurement His 


Is propose d vibratiol 


probler 
is particular \ 
relation- L* ¥ ' D 
some further or plates subjected to edg 
hu’s work needed per unit length Equation 
technique more fully, 
further theory 


ormation is twofold 





LATERAL VIBRATIONS AS RELATED TO STRUCTURAL STABILITY 
Hence from Eq ns >| are it can be concluded that 
10 


(s L* is a function of the end load, the right side of this equation 
will be a function of the end load. Hence Equation [10] will 
give the general rel stionship between the square ol the frequency 
and the end load d mi ire both functions of z, | quation 
10) indicates that a1 ondition for (w,*/w to be pro 
portional to the end load i.e., the mode of free 
vibration must be ’ d by the application of an end load 
The value of d load at which the natural frequency is 
reduced to zero can be deduced from Equation [9 For w,*? to 
he zero, L* must vanis Putting this condition in Equation 


‘ the natur 


frequen Oo beam will become zero when the 


, An examination 

end load assumes ilue giver . 
two modes can be ¢ 

nometric functions vanist 


ship to hold quatio 


But the values of P given bv t equation correspond to the s 


hue kling loads (7 ) natural lrequency is a 
duced to zero , d load approaches the critical buckling 


' Phese two requirements 
load 


: hyperboh functions or sine 
Vanishes in the case of } } 
functions ) ) ‘ conditions will yield purely 
hyperboli lu ws. h t sine soluti 1e only possible 


conditions f \ i correspond to supported 


v t E - ‘ ends. It 





Again the values of \ in as given by this equation 
correspond to the static buckling loads of the plate (1 

Hence the apparently different physical problems of free where 
lateral vibration and elastic stability are limiting cases of a single 


phenomenon, the more general expression being the mode of buckling 
vibration with end thrus The free modes of vibration and the 
buckling modes correspond to the conditions that P = O and : 

orresponding to ntl 


This illustrates that the sq f the frequency varies linearly 


APPLICATION TO STRAIGHT BEAMS vith the end thrust ‘ frequency becomes the 
‘ ; itural frequeney for fre ibt 1; at zero frequency, the end 
Having derived an expression relating the square of the fre- ; 
: thrust is equal to the buckling 
yuency to the end load in general terms, it is now advantageous to 


} 


itions of this result are quite consider 
examine in greater detail the form that this equation assumes in . saiciere 


itural frequen r vibration for 


the case of straight beams 
Age deen dase a » end load and then with an end load 


The harmonic vibrations of a beam with end thrust P are gov . 
i} , he buckling lou two points can 
erned by the eqt 
to intersect 
vill correspond 
buckling load 
on cauatio dicates that t . 
Therefore, from Equatior Equatior indicates that the line 
ind P will no longer be valid. The que 


where 


ruin : CONSTANT 


Tt is Proceeding in a sin r manner to that for straight he 
1c solution o 
possible to derive the corresponding equations for a flat 
= - A a , If the plate is subjected to compressive end thrusts NV, and N, per 
unit length, and end shears N,, per unit length, the dynam 


equ ition becomes (1) 
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the 


ibration modes of su 
18 ’ 
ar? sion per unit length 
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vi 


mM 
wration of a membrane « 
in addition to being ider 

This was proved by Schae 


the following general result 


where D is the flexural rig and yw is the formls b 
rand Haver 

lead to the equation Hence 
thin plate ot poly gonal shape and unifort 

the edges and 


irequency @* 


mass per unit area 
Again harmonic oscillations 
10 p supported along al 
thrust N per unit length, the 


* . 
w 
It is interesting to note th: 


ra polygonal plate with an 


not valid in the limit when the 
The boundary would then be 


tion only, Equa- ‘ - 
conditions no longer would be 

20) redu ; : 
the deviation from linearity i 
irbitrary planform would be « 
seen readily in the ist 


bration and of buck 


lasticl 
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supported 


is evident from Equation 
angular plate uniformly compressed in one direction, the 
ion can be written as 

t boundary 
hdd appre 
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} " ‘ 1 ’ vecurate represe 

e concluded that Equation [21] would yield a linear ep 
the fundamental mode 
| equen 


Howe ver 
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. Phe amplitudes assumed in Ray 
-direction, sim : 
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by energy calcu 
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s} ip between w,? and Ma 
conditions under which this linear law is still exact, viz 


ite uniformly compressed in the 


yerpendi to the c . 
on. he determined easily 


opposite sides 

} ne re ¢ 

d having various edge conditions : 
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The equation of motion is obtained * 
ind so the method gives ap upper Dou 

principle, the approximation to the fund 


much closer than the 


to the true amp! 
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The energy met 
a This is 
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forces, each of which if ting separate 

value problen rl in general 
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cluded from Rayleigh’s principle that 


V>PW 


P w 


Pu w,? 
Equation [28)} gives a lower limit to the frequency of the axiall 
1 ~ 1 


loaded beam. The equality sign will hold if the modes are identi- 


cal in all three cases. This checks the previous result obtained 
for uniform beams which are simply supported. 

Now Rayleigh’s principle can be made to yield an upper limit 
to the fundamental frequency by the use of an assumed mode of 
ind total potential 


Equating the kinetix energies of a 


vibration 
beam 


s a 
EI 
tl 


Phis relationship is exa 


wtual deflection function of 


the vibrating strut, i.e., that which satisfies the differential equa- 


tion. However, it results in a good approximation if the differ- 


ential equation is disregarded entirely and v is anv function of z 


which satisfies the boundary conditions and which is a moderately 


accurate representation of the vibrating mode. By using an 
approximate functi the relationship for uniform beams 


hecomes 


Equation [30] gives an upper limit to the frequency of the axi 
i } 1 


ally loaded beam. Again the equality sign will hold only if the 


modes are identics n which case it chee 


ised to find the 


ks with Equation [28 


This equation can be maximum deviation of 


ind w? for particular cases It is obvious that 


line rity betwee! P 
if the buckling mode is used as the approximate deflection func 


tion o(x), the coefficient of the PP, term in Equation [30] will 


be unity. Similarly, the free vibration mode can be used for 


v(r), in h case the coefficient of the (w*/w)* term will be 
unity 
As an example of this upper limit in a particular case, consider a 


The fundamental buckling mode is 


(1 on 


Substituting this into Equation [30] the expression becomes 


beam fixed at each end 


3pw**l* Pl 


l6nrtkl inkl 
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STABILITY 


But for a clamped bar 


There 


fore Equation [31)can be written as 


5(“) P 
1 > 0.9635 + 
, @ P 


Now let the deflection function be represented by the mode of free 
it the middle instead of at the 
takes the 


Vibrations If the origin be taken 
ends of a beam clamped at both ends, the value of v(z 


form (1 


cosh mr 
| 
cosh 
» 
4.730. Substituting this into | 


where juation 4) 


pression becomes 


0.9704 


If Equations [28], [32], and [33] are plotted as in Fig. 1, the true 


curve of (a@,*/a, 


bounded by the 


against P/P, will lie within the triangle 


three straight lines Hence the shape of the 


probable curve may be estimated from the figure t is noticed 


that the deviation from linearity is extremely small 


For any other boundary conditions, similarly smal! deviations 


from linearity would be expected. In all cases Rayleigh’s 


method gives a very close approximation both to the natural fre- 
quency if the buckling mode is used, and also to the critical load if 
the free I 


Vibratior 


mode is used This method likewise can be 








Frequency Wirn Enp Loap sy Eneror 


MertHops 


VARIATION OF 





200 JOURNAL OF 
applied to flat plates with arbitrary boundary conditions, with 


similar conclusions 
EXPERIMENTAL Tests ON Beams 


In order to test elastically restrained columns, Chu (6) had 


investigated the case of a rectangular frame subjected to end 
thrusts. Although he obtained a linear relationship between the 
load and the square ol the Irequency the results were not con 
Theoretically the should buckle in the 


energy state, which corresponds to a symmetrical mode. 


lowest 
How 
ever, the buckling load corresponding to the extrapolated value of 


clusive frame 


Chu’s curve was approximately equal to the critical load for un- 


symmetrical buckling, Fig. 2. Owing to local vielding at the 


bic. 2 anp UnsyMMeTRICcAL Mopes 


SYMMETRIC AI 


Knife-edges under application ol load, the d imping had increased 
in extent that the test could be conducted only up to 

half the load. Chu that 
measurements that the was 


in the unsymmetrical mode of vibration 


to such 
remarked 


buckling 
indicated 


ipproximatels 
although the structure 
stable one could not 


however, conclude that the structure would buckle in that mode 
first 

In view of the foregoing tests, it was decided to carry out fur 
ther tests on a similar rectangular frame However, the frame 
supports would be improved, and a more refined technique of ex- 


citation would be developed. It was considered desirable to con- 


duet the tests right up to the physical buckling load in order to 


correlate the vibration mode with the actual buckling mode 


5 


The frame, made of 5/i-in. X 1-in. cold-rolled-steel strip brazed 


together at the four corners, measured approximately 12 in. > 


20 i The 


hardened-steel seats 


load was applied through knife-edges resting in 


Instead of measuring the free vibrations of 
the frame, an oscillator was built which would impose an exciting 


force on the frame By varving the frequency of the oscillator so 


as to develop resonance the natural frequency could be deter- 


This was arranged by con 
The ar- 


strain gages were at- 


mined at each loading increment 


structing a magnetic circuit with an alternating flux 


rangement is shown in Fig. 3. Dynamic 


tached to both vertical members of the frame in order to deter- 


mine whether the vibrations were occurring in the symmetrical or 


unsymmetrical modes. Measurements were recorded auto- 


itically on a two-channel recording osc illogr iph. 


The 
Chu’s results could be repeated 


first test was conducted in order to determine whether 
The first question to be settled 
was whether the frame, with a stable unsymmetrical vibration 
mode, would remain stable above the symmetrical buckling load 
Load was applied in increments, and at eac h load the electro 
The 
varied until the occurrence of resonance, which was evidenced by 
It was found that the 


magnet was energized frequency of the oscillator was 


maximum amplitude on the oscillograph 


frame vibrated in the unsymmetrical mode, corroborating Chu’s 


results. The load was increased until it exceeded the theoretical 


symmetrical buckling load of 7330 Ib, and the frame remained 
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stable Frequency measurements could be recorded above this 
mode of 
ibility 


unsymmetrical mode 


load In order to determine the buckling, the load was 


It was found 
The 
Fig. 4. kx- 
j 


gives a buckling load of 


increased continuously until inst resulted 


that the frame buckled in the varia- 


tion of frequency with load is shown as curve 1, 
trapolating the curve to zero frequen 
11,150 Ib, 


oad of 11,360 Ib for this frame 


theoretical unsymmetrical critical 
The 


curve tor 


is compared to the 
test also confirms that the 


deviation from linearity of the such an elastically re 


strained strut is so small as to be negligible in comparison with 
experimental error 
Having established the fact that the mode of buckling is of the 


same form as the mode of vibration, a second test was conducted 
Similar results 
The actual 


as compared with the 


or t second frame ot the same dimensions 


occurred; the points are plotted as curve 2 in Fig. 4 


buckling load was measured as 10,250 lb, 


extrapolated value of 11,450 lb and the theoretical value of 11,360 
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Ib In this test an endeavor was made to obtain some points 


representing the symmetrical mode of vibration, but at most fre- 


quencies insufficient power was available However, two points 


were obtained, and these are plotted as curve 3 in Fig. 4 The 

fact that this curve lies above that for unsymmetrical vibrations 

seems to imply that, for some reason, it is a higher-energy stat 
The third frame was constructed of similar dimensions, except 


that the knife-edge seats were arranged eccentrically by moving 


them 0.2 in. toward the center line of the frame This would 


insure symmetrical buckling, and it was desired to determine the 
vibrations were 


corresponding vibration modes, However, the 


still unsymmetrical, and the points are plotted as curve 4 in Fig 
4. The 


In order to discover this apparent discrepancy with theory, it 


measured buckling load was 5500 Ib 


seemed logical to examine the experimental setup in the light of 


the theoretical assumptions. The boundary conditions would 


first come under suspicion All four knife-edges were firmly sup 
ported in V-blocks which would prevent any lateral movement of 
the four corners. Although this appeared to be a minor point 
the support system was modified. The two left-hand knife-edge 
90-deg V-blocks 


hand supports were replaced by hardened-steel knife-edge seats 


supports were as previously, but the right- 


Hence lateral movement of the right-hand supports was now per 


mitted. The modified system is shown in Fig. 5 


Mopiriep ARRANGEMENT Pic. 6 
or Knire-Epces 


The frame tested with the modified boundary conditions was 


similar to those in the first two tests. This time, however, the 


became the stable configuration for all load- 


symmetrical mod: 
The variation of frequency with load is shown as curve 5, Fig. 4 
Linearity is verified to within experimental limits, the extrapo- 
lated value of 7330 Ib corresponding exactly to the theoretica 
buckling load rhe frame actually failed at 6950 Ib in the sym 
metrical mode 

In view of the results of this test, the reason for the behavior of 
The effect 
of maintaining the horizontal distance between the knife-edges is 
sufficient to cause the unsymmetrical mode to be the stable mods 
These 


the very important result that the buckling load of such a frame 


the frames tested previously now becomes apparent 


both for vibrations and buckling tests thus showed up 


can be predicted by simply measuring the frequency of vibration 


under various loads. The correct mode of buckling will be satis- 
fied automatically owing to the fact that the vibrations will also 
occur in that mode It is seen now that Chu’s extrapolated value 


have corresponded to the actual critical load for his 


would in fa 
boundary conditions had he continued loading up to buckling 
Two important points must be realized in using this experi- 
mental method. The first is that the members of the frame must 
be long enough to fail by elastic instability, i.e., no short column 
effects Then the 
will not vary with loading 
portant in static buckling tests, because the equivalent end fixity 
tests, 


boundary conditions must be such that they 
This second condition is not im- 
at buckling can be computed. However, in the dynamic 
it 1s important that the boundary conditions for each vibration 
test be identical to the boundary conditions which would prevail 


if buckling were to occur 


LATERAL VIBRATIONS AS RELATED 


PO STRUCTURAL STABILITY 


ApPLicaTIoN TO Ricip-Joint Truss 


In the 


mately linear relationship between the load and the square 


previous sections it Was pointed out that the a proxi 


ot the 


frequency could be expected to apply only to elastic bodies whose 


boundary conditions did not vary with load The rectangular 


omplied with this condition as the hori 


ind hence the buckling ver 


frame which was tested 


zontal members carned no load 
} 


rs were restrained by the con nt elastic re 


vided by the horizontal members the case of a joint 


however, the axi member increases with in 


loes not remait 


truss 


creasing load izidity at the end of an 


particular member constant Consequently the 


deviation from the linear relationship can 


main small in such a structure 


In order to ascertain the effect of the 


vith load, it was decided to conduct tests on simpk 


The truss configuration which was chosen is sho 


trusses 


grammatically in Fig. 6 The length of all members was 12 in. and 


section os mild steel, and al 
Hardened-stec 
knife-edges at A and ( 


were allowed to pivot, allowing 


« in, square The material w 


rigid connections 
knife-edges and seats were used Phe 


it B and | 


lateral motion at these points 


joints were brazed to make 


were fixed, while those 


ittaching a vibrator to the center of, 
This 
By varying 


Phe frame was excited by 
member Al 


vibrator was a modihed permanent-magnet spe aker 


and normal to, the in the plane of the truss 


the frequency of energization of the voice coil, the diaphragm of 
the speaker could be made to vibrate over a considerable range o 
tion was measured by d 

AB, BC, and BD, 
The plete test 


a 


frequencies. The frequency of vibr 


namic strain gages attached to the members 


and recording the output on an osc illogr apt 


setup is shown in Fig. 7 





Jorn’ Tauss Test Serve 


he frame was then loaded in increments, and the lowest reso 


nant frequency recorded at each load It was found that each of 
the instrumented members had a different resonant f equency for 
uvidition there was one fre- 


This 


frequency for the complete truss was the one which decreased svs 


1 particular load on the truss, but in 


quency at which all the members resonated resonant 


tematically with load, while the individual frequencies seemed to 
little Th IR8O Ih. the 


suddenly failed by torsional instability out of its plane 


vary load was increased until, at frame 


This un 
expected failure prevented the test from being continued to the 


} 


desired failure in the plane of the truss. The test results are 


plotted in Fig. 8 as curve 1 
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cond frame was then constructed identi 


al to the first, with 


eption that the members were twice as wide, viz., '/: in. x 
In addition, lateral restraint was provided by denying the 
mualing bar These 


seen in t ind again the lowest 


mm 
‘ 


op | possible angular rotation restraints 


ean be The frame was loaded 


th 


failed at 


ig. 7 
resonant frequency ! 


buckling in the plane of the truss 


or the complete frame was recorded at each 
$945 lb when the AE began 
Phe plotted in 
As would be expected in analyzing the data 


all the 


load frame member 


test results are 


Fig. 8 as curve 2 


it was iound that at resonance members vibrated either in 


phase or 180 deg out of phase 
Che actual buckling load was lower than the extrapolated value 
ly 5600 Ib, and this is apparently due to slight im- 
the The theoretical buckl 
computed by a method de veloped by Winter (10) and 
Bijlaard (11), and yielded a value of 4150 Ib. That 
the carrying capacity of this truss would exceed this value was 
igain to be expected in view of the fact that the two points A and 
C 


of approximate 


] l 
pertectc ir truss 1g load tor this 
truss 


Wis 


extended by 


were prevented from moving, whereas the theory assumed no 
restraint rhe restraints were necessary in order to prevent ur 
due complication in the test setup 

An examination of Fig. 8 vields some interesting results As 
predicted, the relationship is not linear 


the fact that 


but more significant 

the deviation from linearity is considerably larger in 
curve 2. It might seem reasonable to deduce from these tests 
that the deviation from linearity increases as the section modulus 
of the members inere 


ines 


This could be interpreted physical! 
as follows: For a’given moment of inertia, as the depth of the 
the ch 


ixial loads in adjoining members is decreased 


members is increased, 1ange in the end restraints due to 


EXPERIMENTAL TrsTs ON PLATES 


ng literature failed to vield any experimenta 


rectangular plates 


ons It appears that the 
have been carried out are on small circular plates 

rophone diaphragm Therefore it was decided to set up 
both unstiffened and stiffened flat 


compression while frequency measure- 


a testing | 


wrogrual 


plates could be 
ments were mace 
The plates used were 


0.040 in 


made of 248-T duralumin sheet, 12 in. 

12 in thick A special ring gage, placed between the 
I } 

the testing machine 


rements of load 


specimen and the platform of was used to 


measure the small in This gage, mounted above 


the test specimen, is shown in Fig. 9 which also shows the vibrator 
used the 
d by 


in truss tests rhe frequencies of the plate were 


measure mounting a strain gage on each side of the plate 
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noninne 


circumlerence 
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The first specimens to be tested were simply supported flat 


plates. The four supports were V-grooves, and the edges of 


the plates were fitted with steel knife-edges The natural fre- 
quency decreased only slightly with increase in load, and a typical 
set of points is plotted in Fig. 10, curve A 

Tests similar 


were then conducted on plate $ loaded 
long clamped edges, had the unloaded edges supported on knife- 
edges and 
Fig. 10, curve B for 


lated by 


which, 


Similar results were obtained, these are shown in 


The natural frequency zero load was caleu 


the method described in reference (12 


In subseque nt 


tests, all four edges were clamped, with similar results 


In order to check the reliability of the experimental 


the 
, the 


setup 
two unloaded edges were allowed to deflect with no restraint 
loaded edges only being clamped. 
Fig. 10 
result to be expected because of 


In this case the 


approximate! 


inear theoretical curve »wed closely, a 


curve ¢ A 
the fa 


late acts 


is u 
wide simple column 


A series of tests were then condu The 


plate with the 


was of the same dimensions as t 
rddition of a single channel stiffener riveted longitudinally 


‘71 P 
ests 


s were 


along 


the center of the plate th various edge 


onditions, but the resu ffened 


As 


under compression do not 


i result of the tests 


the 


A possible explanatior 


foregoing ippeal I thin plates 


lollow rela- 
thin 
ire in the plate 


joad-Trequency 


might 


that tk rvat 


tion in practice iM 


that 
plate is never perfectly flat, so 
effect 
the 


Was also borne out u 


has some The tests indicates that the buckles 
beginning of loading, and grad inere This 


ind Batdort 


grow with ist 
tact 


ion equations should be used in stu 
1 


who showed theoretical 


irge-detle 
f small devia 
, "i 
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tween the fre quency 


recent 


Massor 


ind end thru 


paper DS net reiathonship 


for the « 
reular plates with 


Initial curvature 


Masso 


ular 


Kar 


ir plate equations, 


ipproxim 
pl ifort maded 1 the 
load-frequen itionship was computed 


his 


various values of 


clamped cir 
The 


initial deflections of 


yur 
v re 


the plate, and results aré 


otted in Fig. 11. The curves are plotted 


rameter proportional to the initial ce flectior iz 
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caused the considerab rom linearity 


Therefore it is in 


tions in such 


using 
The 


to estimate th 


justification for 


equ d in practice nethod of 


using frequency m I ts in ler buckling 


load cannot be factory results in the 
plates in 
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thickness o 


compres 
where 
the 


eliminate to test the 


the difficulty might be 
asuring the frequency with each increment of 
, 


curvature will be decreased wit! 


sng ir ¢ juations may tx 


rapolation zero frequetr along the negativ 


thus give the theoret 


For linear sy irequency of lat 
ions is practically line ited to the end thrust 
tionship is exa 

to the bucklin 

determine tl 

vided the 


our 


g the 


straight-line 
to the 


t Sal re 
tin jua 


extrapolation t t ency gir good appre iatior 


buckling loac h etermination il load frequen 


load 


corresponding ) ‘ t oundary nditions which preval 


measuremet t iy of predict 


whereas a theoret unjustifiab| 


tain conditions which are not realized in practice 


In cases where onditions do change with load 


is in rigid-joint trusses thod seems to be applicable under 


tlized reh in th on might 


these li ons I 


restricted conditions. is direst 
be valuable in 


curve were est 


iblished as 


the curvature « irequency* 
1 


symmetrical, 


then the buckling load could be predicted from 


tests up to approximately half that load 


On the other hand, the method does not seem to have a practi- 
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Both free vibrations and forced motions due to cross- 
winds may create important problems in the design of pipe 
lines supported above ground. An analytic investigation, 
based on simple beam theory, shows that the flow of fluid 
in such a pipe line has no beneficial effect upon the vibra- 
tions. The fluid velocity causes a dynamic coupling of 
the simple modes of vibration so that the normal modes 
of vibration are of complex shape with 90-deg out-of-phase 
components. The solution is presented for free vibra- 
tions and for steady-state forced vibrations, and it is shown 
that large amplitudes may be developed if the amount of 
damping is too small. It is shown that at low fluid veloci- 
ties there is negligible effect upon the vibration of the pipe 
line, and at a certain high critical velocity the fluid flow 
causes a dynamic instability. The present analysis re- 
vises the conclusions which appeared in an earlier publica- 
tion. 
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whose total mas ioment of 


s, pipe plus fluid, per unit length is 


inert id modulus of elasticity F Let p be the mass 


fluid 


ia of pipe J, a1 


of the nd let it have a constant city t 


pipe, tl 


in lengt! 


per unit lengt) " 


relative to the ting the effect the 


the 


ordinate me 


Is egies upon 


change of the pipe during vibration 


} 


long 
It isd 


tion describing the m« 


' and the verti 
esired to derive the 
the As the 


center the pipe 


complicated motion owing to the fact that each point or 


assured © pipe be ‘ displ ue 


ment of the pipe differential equa- 


tion of pipe fluid flows along 


the curve described by the line of has a rather 


the center 


vertical acceleration, an angular acceleration, and a 


With 
Hami 
of the motion is overlooked 

Ax 


angent 


ine has a 


hanging curvature a problem of this complexit it is 


iple 


convenient to utilize *to insure that no part 


ton’s prir 


an element of fluid flows along the velocity 
to the 
Assumir 


f the 


pipe 


t ind a velocity vertical 


the 


pit 4X18 4 


small displacements of 


and the 


direction g Ze 


component o velocity is “component 


The 


plus that of the fluid, which is 


kinetic energy in a lengt! 


jl 
12 


The rm the right side represents the kinetic energy of 
the pipe shell 
fluid. The 


theo 


ind the second represents the kinetic energy of the 
in the 


strain energy length dz is, according to beam 


anics- Dynamic 


and New 


st 


Mect 
Nostr 


See, for instance 
Hudson 


D Co 


Housner and D 
1950, p. 260 
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each term may be integrated by parts so as to eliminate the vari- 


ous derivatives of 6z. When this is done, there is obtained 


Jus r, (mea 2pr 2 pv*z 


where all the integrated terms have disappeared because of 


Elz") 6z dx dt = 0 


the 


boundary conditions. Since 6z is arbitrary, the only way for the 


integral to be zero is for the « xpression in square brackets to be 


identically zero Equating this to zero gives the differential 


>* 22 . 0*z 

El ; . . {1] 
ax * a Z ol* 

that the beam 


his equation states is acted upon by three dif- 


forces The first term on the right represents the 
d with the 


the beam; 


ferent inertia 


inertia force associa change in direction of v, enforced 


by the 
acceleration because it travels along a curved path. 


curvature of that is, the fluid experiences an 
The second 
term is the inertia force associated with the Coriolis acceleration 
the fluid is flowing with velocity v relative to 
(Otozr 
The last term represents the 


which arises because 
the pipe, while the pipe itself has an angular velocity (0? 
along its length. 


at any point 


sociated with the vertical acceleration of the pipe 


inertia force 
Free VIBRATIONS 


Equation {1] differs from the more usual 


contains a mixed derivative 


Vibration equations 
) Orol Although this 
erm contains a first derivative with re spect to time, its influence 
which } 


that it 


s not the same as a viscous damping term 


) t rhis can be seen by writing the equation in the forn 


072 
>t? 


The 


the solutions of 


considered to be an applied force free 


normal modes which are 


t) set equal to zero For a simply sup 


onstrained so as to move in the nth 


dt then the applied force is propor- 


force is proportion al to the displ wen 
g out of phase with the velocity so that the vibrational 
of the beam is dissipated in doing work against the force 


f F(r-t force is propor- 


lisplac ement but 


is case is not proportional to the displacement. 


ris syt 


il about the mid-point 


intisymmetrical and therefore 
ot associated with damp- 


072 Orol 


and it is impossible for a s} 


is that it represents dy- 


nami ou mmetrical mode 
of Equation [2] t xist by i f free vibration, for corre- 
sponding wit ymmetrical force, 90 


deg out of pha ‘ ll the antisymmetrical modes, 
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and an antisymmetrical mode produces an F(z, t) which excites 
all of the symmetrical modes 
1} must be 


Thus a normal mode of Equation 


les of the 


a linear combination of all the 


homogeneous form of Equation [2] witl and anti- 


symmetrical modes 90 deg out of phase supported 


span thus has normal modes z;, satisfying I 


=> A 


n 


tn—1 SIN (2n 


rhe normal mode which interact with one 


Zz; Is composed of terms 


inother through the mixed derivative tern 


An odd term, say 
’ 


must be in dynamic balance with the rees exerted by 


all the even terms, and an even term must h the 


coupling forces exerted by all the odd t Jach coefficient 
A, is coupled with every other coefficient 
The coefficients and the natural! freque letermined 


When Expression [¢ 2 ibstituted in I 


1] the mixed derivative g 


as follows juation 


With these substitutions 
collected in groups according 


oeff 
satisfy |] 


equations, for n 
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If the 
frequency equation whose roots determine the natural frequen- 


i's are eliminated from these equations, there results the 
cles Ww 1e normal modes; then = | these equations 
all the 
modes are suppress¢ d except the first two, the equations ire solved 


The 


setting A 


determine the remaining A for each normal mode. If 


readily. frequency equation is 


32prw 


and the coefficient 


the trans-Arabian pipe line* these give 


for the fundamentai mode 


For the parameters of 


4; = 0.06, w = 15.6 


These values are based on a fluid velocity of 15 fps which is 
the figure used in the previous paper* and which is somewhat 


larger than the actual! velocity in the line However, it is seen 


that the 15-fps vel = not large enough to excite the unsym- 


metrical A; compone ippreciably, and it has a negligible effect 


upon the natural frequency which with zero fluid velocity is 


@ 15.6 to three significant figures.‘ It may be verified that 


the strength of the coupling decreases rapidly, so that if A;, A 


and A, are 1, the magnitude of A, is small compared to 
that of A, 

namic coupling is 
efficient A, is the largest and the 
Ay 


script goes AWA 


mode of vibration. In general, the dy- 


it in the ith mode of vibration the co- 


coefficients A i and 


\y i gnitude as the 


sul 


Forcep VIBRATIONS 


I juation |6) may be used to represent 


the solution state forced vibrations. If the aerody- 


namic forces on the pipe are assumed to be distributed uniformly 


along the span and vary sinusoidally with time, the differential 


equation ot m 
a! 

El 
art 


has been expanded in a Fourier series 
leads to the 


where the load p s wt 
Substituting the normal modes of Equation 16 


equations 


Bending Vibrations of a Pipe Line Containing Flowing Fluid 

by H. Ashley and G. Haviland, Jovrnat or Apriiep Mecnanics 

Trans. ASME, vol. 72, 1950, pp. 229-232 
‘ This is the frequency using J] = 0.125 ft* 


A PIPE LINE CONTAINING FLOWING 


‘ 
Aw| El(2k)* (*) per 


FLULD 


2prw » 
l 
Asx, thus 


These equations determine the coefficients A, and 


steady-state vibrations 


it means that 


giving the complete solution for the 
When w is such that these equations give A = © 
the aerodynamic force is in resonance with a mode and the ampli- 
tude becomes increasingly large 

There is a singular case when the amplitudes may become in- 


finite, namely, when the fluid velocity has the critical value 


El x* 
pl? 


This may be seen readily from the differential equation for static 


loading which is 


O*z , o% 
El] > pr p 
or" or? 


This is the same as the equation for the pipe carrying a lateral! 
load p and an axial load ? = pr? It follows that the pipe will 


buckle when P reaches Eulers load, that is 


rE! 


j2 


pv* 


which gives rise to the buckling, represents the inertia forces pro 


duced by the curvature The critical velocity is, of course, 


large 


very 


and is not of importance for pipe lines 


The foregoing analysis shows that if the pipe line is undamped, 


resonance will occur when the frequency of the aerodynamic 


requency of vibration, and ampli- 


In the 


force coincides with a natura 


tudes of unlimited magnitude will develop actual pipe 
juently the 
will de 


amount of damping present and the magnitude of 


line there will be, of course, some 


damping and conse: 


amplitude will remain finite Phe amplitude attained 


pend on the 


the exciting for 


Phus it is seen that the vibration of the pipe 


line is a problem of the same sort as the vibration of the Tacoma 


Narrows Bridge, and if the pipe line is constructed so that it has 


a very small amount of umping, then it is possible that rela 


tively large amplitudes may be built up 


the parameters for the trans 


Arabi 


has only a very 


The values of pipe line 


being such that the fluid v it small effect, it 
is of some interest to examine the explicit solution for zero fluid 
velocity 


Assuming an exciting force ¢ sin wt acting upon a simply sup- 


ported span and setting in Equation [1] results in the 


steady-state forced vibration 


osh 2 A 
cosh A 


where 


‘ | ma? (: ‘ 
Vir (3) 
12,1 = 66 

So the reported data im- 


ll, w/(29) = 
13 Ib per ft 


Using the values m = O,z = '/, 


in., there is obtained q = 
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13 sin 2.4 tif the differential Equa 


plies an exciting force F 
tion [1] describes the motion 

A force of 13 Ib per [t is too large to be reasonably expec ted 
However, the foregoing analysis does not take into account an 
axial force which likely is in the pipe As is well known, the 
effect of 
vibration 


an axial force is to decrease the natural frequency of 
The fluid pressure in the pipe is not sufficient to pro- 
duce an appreciable change in ireque ne but if the pipe line Is 
arge temperature-expansion stresses exist, the 


fluid 


pressure in the trans-Arabian pipe line is 500 psi and the net axial 


fabricated so that 
frequency can be changed greatly. For example, if the 
compressive stresses are 20,000 psi, then the natural frequency is 
reduced from 2.5 eps to 1.2 eps which coincides with the frequency 
of the aerodynamic force. It is not unlikely that the axial force 
in the trans-Arabian pipe line is much larger than it is thought 
to be, and a relatively small aerodynamic force is approximately 
If this is the case, 
vibration can be controlled by the installation of vibration damp- 


in resonance with the pipe the amplitude of 


ers along the line. However, the essential question yet remain- 
What magnitude of aerodynamic force is possible when air 
This 
must be investigated if the problem of pipe-line vibrations is to be 


ing is 
is flowing past a cylinder which is near a plane boundary? 


solved completely ° 


Appendix 


An investigation of this problem also appeared in an earlier 
paper®. The differential equation studied there was 


072 o? 


) 
El 
0. 


Oral dt? 


his equation he glects the inertia forces produce d by the curva- 


ture of the pipe and includes only one half of the inertia forces 


due to the Coriolis acceleration, so that it underestimates the in- 


fluence of the fluid velocity. The solution of this equation is of the 
same form as that of Equation [1], and for a simply supported 


beam it. is the sum of normal modes 
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rhe coefficients and the natural frequencies w, are determined by 


substituting the ith mode in Equation z the Expansions 


and collecting terms according to w 


Equating the coefficients to zero gives the « 


u(7) 


El(2n 


Phese equations determine the frequet uw normal 


mode, and setting A; = | they determine the remaining A It 


follows that the solution of Equation [9] is the sum of the normal 
Since the differential Equation [9] 


modes z, of Equation [10 


is that of a conservative system, the solution shows, as it should 


that there is no energy loss, and the free vibrations once estab- 
lished are permanent. In the pre 
that the flow of fluid in such a pipe line produces marked damp 


The 


in the previous paper was based on a 


vious paper it was cone luded 


ing tendencies and thus may reduce theseverity of loading 


solution of Equation [9 
approximation and in with the 


orrespondence 


power-series 
authors it has been agreed that if the power series contained a 
to Equa- 


the cor. 


sufficient number of terms, the solution w 


tion [10 It also has been agreed that Equat 


rect description of the vibrations 





Heat Conduction in a Compressible Fluid 


By J. D. COLE! ann T. Y 

The conduction of heat in a compressible fluid medium 
is studied in this paper. The compressibility of the fluid 
has the important effect that an excess of heat is not only 
carried away by conduction but also causes a pressure 
wave to move in the fluid. A closely related problem is the 
effect of heat conduction on the propagation of sound 
waves in a gas. In both cases the total energy is divided 


between thermal and pressure waves. 
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effects ol 


nduction in & compress 


omitted. This, of course, is not a 


realistic any kinetic theory of matter indi- 


‘ nption since 


ates a close relationship between the coefficients of heat conduc- 


tion and it The main justification for the omission is 


viscos 
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tion of mass, momentum, and energy applied toacontinuum. It is 


issumed that the medium is a perfect gas, so that the perfect gas 


| iw 
Rol 


holds 


be written 


and viscosit ted. The basic equations car 
is follows 
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solu 
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20 


d [32] w 


36), and using Equations 


*) Vie 
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The desired fund solution T, 


complex inversion Equatior 
it 1 


imental > is given by applying the 
to Equations [42] and [43] and 


with some rather compli ated 


2) 
s seen that 
integr ils 


we have to deal 


The success of the method depends very much on 
applying conformal transformations of the complex plane to 
simplify these integrals. If we consider I and successively 
the transformations 


apply 
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The fa 


sent the fund 


toring of the oy 


wrator L-quation 33} permits us to rey 
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we obtain the following result 


where the contours C,,. are shown in Fig. 1. @; is a contour start 


ing at stan ingle of 15 deg with the real axis ending up at 


at an angle of +45 deg, passing the real axis to the right of the 


singularity at z l @, starts 


essential 
place but passes to the left of z = 1. Thus we have the important 
result that there is a discontinuity in the representation of the 
fundamental solutions across the characteristic of the systems in- 


t For 2 c,t, © 
1 the 


dicated Equation 6 = 2 can be de 


formed into the imaginary axis anc fundamental solutior 


represented as a real integral 


APPLIED 


ind ends in the same 


MECHANICS 


For |z 


ind it can be shown that most of the contribution come 


> ¢,t the contour @; can be made to pass through 


neighborhood of this point on a contour where 


48] and [49 wher 


and density 


Equations are good 


sure velocity are given by 


 yr(*)! 


sign 


2 qy7(7,. ) at 


uations 


fundamenta 
ind 


Of the 


iS 19) illustr 1e remarks 


roduction unit amount of heat introducs 


ne mere 18], 


ly diffuses about the origin, Equation 
conduction but some diffuses about the 


19), 


seen by integrating Equations 


the usual lines of propagation of sound waves 
18} and [49] from z 


, and using 
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that 1/7 of total heat diffuses about the origin, with the co- 
efficient of diffusion «/y; 1 1/¥ of total heat diffuses about 
r = +ct, with the coefficient of diffusion [(7 1)/(2y)\« As 
«k becomes smaller all of these disturbances concentrate about 
their centers and approach 6-functions as « — 0. By examining 
the formulas for pressure, density, and velocity, Equations [50 

51!, and [52], it can be seen that there are appreciable pressure 
und velocity changes associated with the waves near 2 +ct, but 
much smaller changes near 2 0. Near z = O the changes are 
mainly in temperature and densit Thus we see how the intro- 
duction of heat at the origin produces both a thermal wave and a 
pressure wave It should be mentioned again that the results in 
the foregoing are asymy ftlarge. For very small values of 


it conduction is the dominant process Ihe pressure 


wages take some short time to form. From our solution we also 


can see that heat conduction causes the attenuation of sound 





waves as 1) ¥ «i Thus heat conduction has an effect similar to 
that of viscosity which causes attenuations like 1, Y vt. Thus the 
picture when both heat conduction and viscosity are considered in 
this problem should be about the same as here 

The solution for higher dimensional waves also can be worked 
out, and the results are very similar with the usual geometrical 
distortions. These cases together with applications to various 
boundary-value problems Wi be given in a future paper 
For boundary-value problems, where solutions are given by 
superposition of fundamental solutions, the results are similar to 


those given here 
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It is important that the data contained in technical papers be made readily available to designing engi- 


neers. 


In order to satisfy these needs of industry, this section of the Journal includes a concise presentation 


of data and information drawn chiefly from papers previously published by the Applied Mechanics Division 


of The American Society of Mechanical Engineers. 





The Computation of Flexural—Torsional 


Buckling Loads 


By H. I 


Solution by ordinary means of the cubic equation for 
the torsional buckling load, first derived by Kappus, often 
A method is 
loss of accuracy involved. 


leads to insufficiently accurate results. 
presented which avoids the 
The procedure is straightforward and no graphical means 
are involved. However, one formula is to be solved by 


trial and error. 
INTRODt CTION 


FT NHE problem of primary instability of columns, including 
flexural and torsional buckling, may be considered as defi- 


nitely solved with the comprehensive treatise by Kappus.? 


The considerations underlying Kappus’ theory, which is an im- 
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Hoff 


same pre 


the gener il case 


ross sections has given a concise strain 


tion of the issing the several 


aerivia 


ling modes that in occur, a rameters 


the same cross-sectional properties t most suita 


n the pre sent paper 
of the authors, however, gave full attention to the nu 


ng load of an unsupported 


given by Kappus? to that 
Research. 
Profil im 

vol. 14 


tructural 


MICHIELSEN,' DUIVENDRECHT 


NETHERLANDS 


topic. Formulas and a nomog permitting 


an easy determination of the criti this nomo 
gram often leads to an unsati ses ever 
insufficient accuracy of the sult his car n from the 
following considerations 

From the nomogram three roots w 
ilways read, Only 


It is obv 


when the se 


ritical stress) are 
tical importance lous, howeve 
teriorates rapidly l 
lowest one 


Now it 


basically 


is well known that a short 


torsional mode, whereas ir flexura 
form (Fuler-l 


symmetrical 


uuckling) will preva If the cross section is almost 


with a single axis of svmmet of tran 


sition between the two modes is ther 


this range 
the lowest two roots of the ec 
stress, are near 
at al If the cros 
two roots coin 
case the cubic 
a linear one 

A direct 
nomogram gives the 
After ¢ 


mentioned to the deter 


difficulty umbet 


roots 


whereafter the 


7/3). It is obvious 
computations no satisia 
The relative accurac 
times unta 
iq, the qu 
determin ‘ ‘~ 
This is so because the w g const 


eas the 


numbers 


to the centroid of the munimun 
(, determining mate buckling stress, may 


, 
nount to a small fract This is the case 
with short 


Hoff as 


In this paper a mett 


or inst 


ince 
vas used by 
satisiactorily accurate 


iphical means 





MICHIELSEN—THE COMPUTATION OF FLEXURAL-TORSIONAL BUCKLING LOADS 
Mernop or Computation Substitution of Equations [5 into Equati 
rhe equation to be solved is given by Kappus* ss 

7, (0,0 p 


Now the quantit 


ix introduced. The physi interpretation of 
minimum value of ¢ This will be discussed further in Append 


equation written as 


p 
Equation [9] is esse tical with Equation [19] of 
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The quantities are detined accordi 


eading to 


dw.nhds 
df .tds 


CENTROID 


47 _— 
x In the same way 





Fig.1 Derix 
First a transformation to the principal axes 
to be made With 


o 


}quation 


Kither 


omput 





JOURNAL OF 


solution of 


Equation [16] or [17 
nd error It will be shown that a rapid approach is nearly al 


vs assured 


u 0, Equation [18] becomes 
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29 
20; 


In the last two of Equations [4] the relative accuracy of 
extremely small values of p,* or p:* can be unsatisfactory, the 
absolute accuracy being of the same order as that of o,* and 
o:*. This, however, is of no importance, as long as (¢,* @;) is 
not of the same order as p;*/o;, or (4,* — o:) of the same order as 
P2*/ G2, 

If, however, (o,* o 


respectively 
— 0, the right-hand 
approaches p,. In that case p; rather 
than p,? should be known with an accuracy comparable with 
that of (a,* + o,) 
\ sufficiently accurate determination of p,? or p:?, when (¢, 
o;) and (4,* o:) are of the same order as p,*/o, and p:"/¢:, re- 
spectively, is made possible through the following transformation 
The last of Equ ations [4] leads to 


member of Equation [28 


The same applies to p,? in Equation [29] 
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Consequently 


Pip: = 


o a 
bats 


With the aid of Equation 
ind p 


30) p, can be calculated from p.: or 


vice versa. Since p uppear only in quadratic form in the 


pertinent equations, the sign is of no importance in Equation 


30) 
CONCLUSIONS 


A sufficiently accurate computation of the primary buckling 
loads of unsupported open-section columns can be performed in 
the following way 

From the properties of the column the parameters @,, 7,, p,,, 
p,, and p, are determined, according to the Equations [2] and 
and 
the 


3], and in agreement with Equations [47], [46], [82], [56] 
[S81] of 


reference 2. Furthermore, ¢,* is calculated from 


Equations [2], [4], and [8 or directly from 
NEC* 
31 
I 


Equation [31] is a modification of Equation [48] of reference 2, 


in so far as the warping constant C with respect to the centroid is 
replaced by the minimum warping constant C* with respect to 
the shear center 

The parameters mentioned are substituted into Equations 
[4], [10], and [13 If necessary, the first Equations ll 
[14] as well as Equations [28], [29], and [30 
prove the accuracy 

A first approximation to the buckling stress a, is then found 


and 
are used to im- 


from Inequality [22]. Solution of Equations [16] or [17] by 


trial and error leads to the actual value of gp. Equation [17] 
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usually is the most efficient one. If it fails to give a satisfactory 


Equation [16] always can serve as a reliable substitute 


Equation [27] should be used instead of Equation 


accuracy, 
lf 3,° 
16 


> Os, 


Appendix | 


Substituting Equations [2] of the present paper into equation 


23) of reference 2 one finds 
rpn\* Ip Ip 
‘ 7 + ( ) 0, 2( \ ) bes 
2(? )p + 


Setting the derivatives with respect to 2, and yp equal to zero, 


one has 


The values of zp and y for which ¢, is a minimum, are then 


found as 


Substitution of (rp/t times the first of Equations [35 


/p/t,) times the second of Equations [33] into Equation 


leads to 


From Equations 


op," + 2p,,P2P, Tp, 


Setting p,, = O and replacing accordingly the subseripts 2 


y by 1 and 2, respectively, one can transform Equation [35 


OG, minF\O = 60,0 o-p o.p 


shows 


Comparison with Equation [7 


o.* =d@ 


If x, and y, are introduced to denote the co-ordinates of the 
shear center with respect t« the prin ipal axes through the cen 
troid of the cross section, Equation [34] leads to 

p 


o 


Through Equations [38| a physical interpretation of p; and p 
is found, in addition to that of o, 


Appendix 2 


Example (a 
Let it be assumed that 








According to Equation [10 
6, = (325.9 
328.25 


a,* 328.25 


From Equation [11 
14.8 X 
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14.8 x 
14.8 
641.75 


14.69 = 328.25 3 


13 


313.56 64 


With the assumption 


Equation [13] yields 


. 38.1 
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50.8 


o.* = 52.795 
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With the aid of Equation [14] one obtains 


14.7 14.69 


103.6 


ss 
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= 2.085 


@, = 52.795 50.71 = 103.5 


Application of Equation [17] results in 


2.085 
627.0 < 76.2 


(641.8 o 103.5 


Substitution of leads to 
2.085 
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No reliable result obtainable. Equation 
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Therefore 0p = 13. 


9430 Exan ple (b 


Let it be assumed that 


é, 
56 
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sumption 
LS 


] quation | 15) yields 


14.695 

14.735 
Application of Equation 
(14.695 0.04 
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Substitutior 


2.085 
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0.850 
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No reliable result obtainable. 


is tried now 
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Therefore, ¢) = 
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Let it be assumed that 
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and Equation 
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15.58 
28.8 X 0.031 
7p) (17.00 


Substitution of leads to 


Sp = 15.7 15.77 
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15.58 Pas 
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From Equation 
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(16.23 Cp 
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16.32 em 
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0.0288 
568 &X 
16.32 16.32 
27.9 & 0.081 0.0355 
0.43 


56.7 


Therefore, o = 15.78 


If Equation [27] 


is used 


Substitution of 
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28.0 0.031 0.0246 


0.62 


1 + 
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Torsional Vibration Modes of Tapered Bars 


By H. E. FETTIS,' DAYTON, OHIO 


The author offers a solution for the differential equation 
which governs the torsional vibration modes of a shaft 
or bar, these modes being simple harmonic functions of 
time. 


4 HE torsional vibration modes of a shaft or bar which are 
simple harmonic functions of time are governed by the 


differential equa 


(o-ORDINAT! 


l solution of | 
spanwise co-ordinate, ft 
Trequenc radians pet sen 
ingular displacement, radians 
ross-sectional mass moment of inertia, in-lb sec? 
torsional rigidity, Ib-ft 
paper is concerned with the solution of Equation [1] for 
ntilever bars in the case where J and GJ are proportional to the 
urth power o the spanwise co-ordinat« It is seen that this 
ist assumption holds exactly f lid linear tapered bars of 
rope material,? a l m be used with a good 
ipproximation even in ¢ here these conditions do 
| gz n aireraft wi with tor hox-t ype 
the fourth-power Inserting the end « 


lowing equation re sulting equations to 
i ! 


\o 
Aa 


values of J and GJ at some convenient relerence 
sumed to be measured from the projected tip of 
Fig 
ratio,” A, | be defined as the ratio of the dimen- 
i bar to that at the free end. Thus the 
two ends are taken to be, respectively, ean onlv have a non 
the length of the bar being L =a 1 A}. the transcendental ec 


must be solved subject to the cor ditions 


4 Xa 0 
a 0 
j Sines 2 is proportior 


ot zero amplitude at the fixed en : 


‘I, Equation Ss 
mines the natural frequer 
Wright-Patterson Equation [8] there is asso 


orresponding solution of equ 
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mode shapes are shown in Fig. 3 (a) and (¢ In Fig. 2(a) and (b CABLE 1 VALUES OF THE FIRST FIVE CHARACTERIST 
: SUMBERS gl, FOR A LINEARLY TAPERED, UNIFORM SOLID 
the variation of the first five characteristic aumbers BL, with CANTILEVER BAR, FOR VARIOUS TAPER RATIOS 


taper ratio is shown In Fig. 3 (a) and (¢ the origin has been 
shifted from the projected tip to the fixed end by means of the 
transiormation 
DN 8 5 ; 19280 


A l Lo ‘ a : “oe 


6.77280 
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The mode shapes in terms of the new co-ordinate z, are defined by 2 . 3 24 


the equati 


1 t+ ABL?[A — (A — 1a] } sin BL 
A (A : ICA 1)? ABL,*| sin BL, 
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root of the characteristic equation 8). 


iuthor also has derived the characteristic equation lor the ol the 


re a concentrated load is attached to either or both ends in the 





Analvsis of Elastic Two and Three- 


Section Short Columns 


By C. M. TYLER, JR., 
Short columns subjected to repeated loading or in which 
permanent deformation is not permissible should be de- 
signed on the basis of maximum fiber stress. This stress 
may vary considerably in apparently identical columns 
owing to small inaccuracies of shape and application of 
load. The effects of all probable deviations of both column 
and load can be represented by an “‘equivalent eccentricity 
of loading’ or by an “‘equivalent curvature of the column” 
1, 2). These equivalents are so chosen that calculated 
fiber stresses will equal the maximum probable stress en- 
countered in the actual column. In the paper, analytical 
formulas in terms of equivalent eccentricity are developed 
for two and three-section short columns. These may be 
considered as extensions of the secant column formula (3 
to more complicated shapes. Specific numerical examples 
of sectional short columns may be analyzed by a relaxation 
method (4). Sectional long columns have been treated by 
Thompson (5) and Radomski (6). 


NOMENCLATURE 
The following nomenclature is used in the paper 
ixial load on column 
of load 


end moment 


eccentricity 
eccentric equals Pe 
bending moments at a and 6 


leflection at a, normal to axis of column 
length co-ordinate of column or of a section 
noment 


bending 

bending moment at zr due to a unit dummy load at ¢ 
applied normal to axis of column 

length of column 

length of section of column 

modulus of elasticity 

moment of inertia 


equals V El P 


subscripts referring to sections of column 


load stiffness parameter 
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\ sectional column subject to an end moment can be analyzed 


apletely if the bending moments at the ends of the sections are 


wn Each section can be treated as a beam-column with 
different moments at each end 
The moment at the end of any section, say, .,, will equal the 


moment V5 due to eccentricity of lo iding plus the produc t of the 
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nk GENERAL Twe 


Syuuerrican Tares 


Equation this expression L, and J, are the length and 


load stiffne tion, and 7 is measured from the end of the 


section 


A sectional column with end moment Mp will determine an « 


pression like Equation [2] for each section. These expressions 
can be substituted in Equation [1] and integrated to yield simul- 
taneous equations in terms of moments at the ends of the sections 


With these 


section can be found from 


moments known the bending moment on any cross 


Equation [2 
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For the 
moment Vo, hx 


two-section column, Fig. l(a), with load P 


juation 1! can be written as 


" rf Mir) Ls 
—. a eS 
“% M2 
fom ( 
aa El 


be put in the same form as the first by 
transformation 2° = L 


The second integral may 
making the 


El, P 


Then letting J? 


V 


M(x) and M(2z 
and integrated in Equation [3 


When expressions for ire obtained from Equa- 


tion [2 the ratio of M, to Mo is 


found to be 


The « 


« similar 


ase of the three-section column, Fig. 1(6) can be solved in 
manner. Two applications of Equation [1] 


and Mo. 
and M, 


yield two 
Solution of 
to VM 


simultaneous equations in terms of W,, M,, 
these equations vields the ratios of VetoM 


i there 


eset 


Vv Vv 


ion column, Fig. I(¢ 


ind Ly, = La, 


is svVmmetrical about 


Mquations 5! and |6] then 


NUMERICAL APPLICATIONS 


For numerical caleulations the designer must decide upon the 
probable equivalent eccentricity of the load on the column, The 
end moment M isequ il to the produc t of the load P and its eecen- 
tricity e. With MW, known the bending moments at the ends of 
the sections can be calculated by Equations [4], [5], [6], and [7 

The results of these equations may be presented graphically. As 
an example, Fig. 2 is a chart for two-section short columns with 
sections of equal length. As an aid in using the equations Table 1 


contains values of 
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from related functions (9 
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derived 


The maximum bending moment of a section may not be at its 


ends. The location of the maximum bending moment of a se¢ 
of length L and M 


pression (8 


tior 


with end moments MV is given by the ex- 


and its magnitude 


These equations are inconvenient to use. Their results are pre- 
sented graphically in Fig. 3 outside the 
range of this graph, the following approximate equation can be 


used 
L ( us) 
ese + 0.10 cot - 
J J 


For cases which are 


Musx MM, 
Mo My 
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Brief 


On the recommendations of the Executive Committee 
of the ASME Applied Mechanics Division, it has been 
decided to initiate a section devoted to brief notes on 
technical matters in mechanics. These notes must not be 
longer than 750 words (about 2'/, double-spaced type- 
written pages, including figures) and will be subject to the 
usual review procedure prior to publication. After ap- 
proval such notes will be published in the next issue of 
the Journal. The notes should be submitted to 
Secretary of the ASME Applied Mechanics Division. 


Notes 


the 


Bending of Orthotropic Beams 


By H. D. CONWAY,' ITHACA, N. ¥ 


stress 


resolve 


pr IBLEMS 


into the solution of the 


orthotropic plane themselves 


differential equation 


) )(°* 
dy? ) 


iditions, 


\ number of orthotropic problems have been solved in a recent 


paper Use of polynomial stress functions will be made here for 


two practically important problems of orthotropic beams, suc 


is those made from wood Phe corresponding isotropic solutions 


ire well known 
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Note on Energy of Bending of Plates 


By H. L. LANGHAAR,' URBANA, ILL 
"THI 


trop 


differential equation for small deflections of a flat 
Viw = q/D 


expression 21 Viw,u w* 


elastic plate does not contain Poisson's 
r itio explicitly, because the * & 
formula for the total potential energy of the 
of the of 


This circumstance reflects the fact that the surface 


which oecurs in the 


plute satisfies Euler's equation calculus variations 
identically 
it tegral of the fore going expression may be transformed into a 
the 


theorem of differential geometry Cc 


line integral on boundary by means of the Gauss-Bonnet 


uses for which the line inte- 
gral cancels from the natural boundary conditions of the plate 


have been discussed in the literature, but seemingly not from the 


most gene ral yu int of view 
There 
irdless of the magnitudes of the deflections, the expression 


for the strain energy of bending of a flat isotropic elastic pl: 


is no need to limit the discussion to small deflections 
Reg 
ite of 


nt thickness ts 


l ’ aa i 1 2v 
( =: dD / / + + dA 
ae 


are the principal radii of curvature of the de- 


onst 


and 7 


which 7 


formed middle surface of the pl ite. The derivation of Equation 


1] employs the assumptions that plane stress exists in any 


mina that is parallel to the middle plane of the plate, and that 


shearing stresses on planes parallel to the middle plane of the 


Furthermore, the strains must be so small that 


area dA be 


the deflections may be large, if the plate is thin 


ite are 


zero 


the differential may considered to be invariant, but 


be written alternatively 


pf fem 
i 


Equation [1] may 


mean curvature” 


middle 


known as the 


f the 


V 


Gaussian 


K 


eurvature’’ o 


Phe quantities and are 
ind the 
of the plate 

Let the 


let the bounding curve of the deformed middle surface be C. 


deformed surface 


of the plate be simply connected, and 
Let 


middle surface 


the curve C be composed of several smooth ares that form exterior 


4,, Bs, at their vertexes; then 


Sf Kaa = 27 —30,— fx, ds 
i ( 


is the geodesic curvature of the curve C on the deformed 
the known the 


theorem,’”? 


ingles 


where Kk 


middle surface of plate. Equation [4] is as 


Gauss- Bonnet 


Theoretical and Mechanics 
Mem. ASMI 
Geometry by D 
Mass., 1950, p. 155 
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of the Society 
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Since an additive constant in the strain energy is of no interest 


Equations [2] and [4] vield 


2D f f weda+D(—v) fxd 
i 


( 


prov ided that the angles 6. are constants 

The absolute value of the geodesic curvature «, of the curve ¢ 
it the point P is the magnitude of the orthogonal projection 0 
the curvature vector K of the curve C on the plane that is tangent 
TI 
vector K is directed along the principal normal to the curve ( 
of the curve ¢ 


to the deformed middle surface of the plate at the point P 
It is equal in magnitude to the curvature 1/R 
»desic is zero Ti 


to the 


By definition, the geodesic curvature of a ge 
K of 


w else 


curvature vector a geodesic is normal deform 


middle surface, it is zero. Consequently, if the curve ¢ 


reduces to 


ares of geodesics Equation 5 


2p f f meas 
i 


is ¢ omposed of 


U;, 


Iwo special cases of Equation 6] are noteworthy 


(a) If the plate is polygonal, and if the deflection w vanishes 


the edges remain straight ! 


at the boundary, Consequently, t 


* are geodesics, and Equation [6] is valid. This case has 


Nada 
j 


bh) If a continuous plate or slab is supported 


edge 
been mentioned by 
by a gridw 
the same si and 


Ze 


of flexible beams, and if all panels have 

same deformation, the deflection curves of the beams are geod 

they 
the 
ipplies for any panel 

by Prof. L. E 
If the 

curve, and its curvature vector K is independent of the 


Then 


ics, since ire plane curves whe principal normals 


normal to defort Consequently | quation 


This « » WAS nted out to the wr 


pei 
Goodman 


deflection is zero at the edge of the plate, ¢ 


leform 


tior 


where dw/dn is the normal derivativ the flection at tl 


boundary of the plate 
An important case of Equation [7] arises, if the plate is clampe 
Then du 


shows that the line integral of x, is 
additive 


at its boundary dn 0. Consequently, Equatior 


ndependent of the de 
the 


valid 


flection, Since an constant in strain energy i- 


irrelevant, Equation [6] is ac 
of a simply connected plate 


This 


ording|s for any type of ce 


that is clamped at it- 


formation 


boundary. conclusion has been derived previously by 


Stodola.* 
A multiply connected plate may be reduced to a simply con 
This device 


the integration theory of complex functions. If 


nected one by means of cross cuts is well known it 
the widths of 
the gaps between pairs of cuts approach zero, the integrals o 


geodesic curvature along the cuts cancel each other Cons 


quently, Equation [4] remains valid for a multiply connects 
plate, if the curve ¢ 


the deformed middle surface of the plate, and the inner parts o 


is interpreted to be the entire boundary 


the boundary are described in the opposite sense to the outs 


It follows that Equation a multip 


provided that the 


boundary. 6] is valid for 


connected plate with straight edges, edges 
remain straight and the angles between consecutive edges re- 
main constant. Likewise, Equation [6] is valid for any multip! 
connected plate with clamped edges 

3 “*Elastische Platten,’ by A. Nadai, Julius Springer, Berlin, Ger 
many, 1922, p. 275. 

4 “Ueber die Schwingungen von Dampfturbinen-Laufradern,” by A 
Stodola, Schweizerische Bauzeitung, 1914, p. 251 








Discussion 


Stresses in Pipe Bundles 
EH. Lee.t 


writer in the 


The following points were not made clear to the 


prese itation of the paper 
1 The interaction between tubes has been considered as a reac- 


tion force acting through the welds. Is it clear from the sym- 
metry of the problem that no couple is required in addition to 
this point force reaction? In general, one would expect such a 


to arise, but it would be zero if there were no tendency for 
tween the tubes at the welds 


couplk 
relative 
2 Itw 


directions a 


rotation ty 


is stated 


that the equivalent elastic constants in two 


t Since there are 


ngles were almost the same 


three directions of symmetry at 120 deg to each other, it seems 


that this might imply complete symmetry, by analogy with the 


similar case of the moment of inertia of an isosceles triangle, where 


the three axes of symme try determine the ellipse of gyration to 


be a circle lf such general symmetry applies in this case, it 


would provide heck on the numerical analvsis 


AvTHors’ CLOSURE 


The 


raised by him 


suthors wish to thank Dr. Lee for the interesting questions 


In reply to the first possibility 


contact points was not considered in 


juestion, it is quite true that the 
ot 


moment reactions at the 


settled before the solution 


lastic 


the paper, and this point should be 


used in the calculation of the equivalent ¢ constants can be 


considered to be complet rhe truth of the matter is that under 
the action of the force reactions the angular deflection of the con- 


tact points is the same for ulpacent pipes and hence there is no 


need for moment reactions 
To see this, consider the tube assembly, s: 
6 of the 


ect to the 


, for tension p, along 
in Fig of this 
Tube I, subj 
circular shape to some other shape 


1, hence 
slope change at point B 


the z-axis, as shown paper, or in Fig. 1 


closure forces F,, F,, distorts from its 
Clearly there is no rotation 
at point no moment arises there. There is, however 
But the tube II deforms to precisely 
the same shape as tube I, since it is subjected to the same pattern 
of forces; hence no relative rotation occurs at the contact point 


B eithe Note 
that this argument, which is also valid for the p,-loading in Fig 


Thus the moment that might arise at B is zero 


By H. Poritsky and G 
1951, issue of the Jo 
vol pp. 241-250 

Associate Professor of 


Providence, R. I. Met 


Horvay, published in the September 


RNAL OF Apptiep Mecnanics, Trans. ASMI 


Applied Mather 
ASME 


ies, Brown University 


fundamental assumption of slow 


f the paper, is based on ou 

nging stresses in the equivalent solid, which underlies the cal 
ulation of the elastic constants E,, E, Me, My G We 
that the state of l that 


ariy 
these stresses over distances of one tube diameter can be con- 


assum 


stress is nes untiorm 80 Variations in 


sidered as insignificant 


The authors are much indebted to Dr. Lee for raising the second 


question, and for suggesting isotropy of the structure in th 


attribute of 


-plane This is an important with 
the calculation of E,, E,, et 


in the original paper was carried out by 


structures 
equilateral triangle symmetry 
analogous but quite dis- 
similar procedures which did not suggest the relations E E. 
M, uy, G = EF,/Al+yn, 


> 


To prove these relations, consider the structure in Fig. 2 of this 
' 


osure, loaded b § direction 


in applied tension p in the a is 


In this direction the F., uw, constants hold; hence, writing 


wu =icos7r/3 jsin zr 3 


v= —isinr/3 +jcosF 
the strain dyadic becomes 
p M,P 
~ wu vw 
E, Ek 
P - 
| ii 
1E, 
M, 
There result therefore the strains 
p 


1E, 


1E, 
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On the other hand, the stress dyadic is ’ ; Bo upd) up) + K; Jol usp 


puu = p Phe ' ij + ji) + 3 , 3 denote real and imaginary parts. The 
determined from the two conditions that 


and is equivalent to the stresses I ace, z = 0 


t,, = 0, o¢, = prescribed function of 


If the cylinder is rather short then, in place of Equatior 


These stresses give rise to the ions of the form 


3z, sinh Bz) [BrJi(Br) + «KJofB 


sI 
i 

E, ) > [6] must be used. 

= 3 p/iG | 


On the Extrusion of a Very Viscous 


Liquid’ 


6 lesired relations 


My 
P.G. Lunp.? The author’s attention is invited to 
From the proof it is also obvious that the particular value, a= . 7 : 
: = “ : rhom.* The type of problem treated in this paper is one of three- 
m7 /3, of the u-direction in which the properties F,, uw, of the i- . 
- Paetie . dimensional flow with axial symmetry, and the specific proble1 
direction repeat is inconsequential once the existence of the , - ; : 
that of sudden enlargement of a circular pipe by a diametra 
second direction wu is established, isotropy of the structure in the am : ; 
. factor of 4 The results are used to give the end correctior 
plane is insured, , 
; , ; capillary-tube measurements of viscosity 
\ third question raised at the time of the presentation ol the 
paper concerned the S, stress system in an isotropic cylinder for 
. y ° ae 
axially symmetric loading. We plan to treat this problem at Plastic-W ave Propagation Effects 
greater length elsewhere, but a brief outline of it here will be of . . . paola ° 
in High-Speed Testing 


interest, & 


‘he problem is the determination of the stress system arising . one 
be 1c |} to sa ssl — ‘ 5 D. 8. Woop.?’ The authors have presented a most useful 
when a self-equilibrating system, o,{r, z), ol normal tractions acts . , 
sentevice= fete & 2 si : . method for determining whether plastic-wave propagation effects 

on an end face, 0, of a long evlinder | 
may be neglected or not in any given experimental arrangement 
R ; for longitudinal impact testing. In the case of impact experi 

o,rdar 0 tf = 

» = ments in which plastic-wave propagation in the specimen must be 
| | bl ‘ light! considered, it might be well to emph size that it may be neces 
he t tment of this problem differs slightly trom th I ut- . 
e treatmen a — oo sary to consider the waves in the hammer and anvil of the testing 
lined at the end of the paper for a nonisotropic body For the ' , : 
. arrangement as well as those in the specimen I niortunately, ar 
isotropic cylinder one assumes a Love stress function ¢ in the 


infinitely rigid material which could be ‘ eliminate this 

form complication does not exist 

KIA B The authors’ example of how plastic strain-wave effects ir 
specimen can lead to an apparent but false increase in energ) 

his satisfies the biharmonic equation sorption per unit volume illustrates an error which has been made 

by several investigators, including the writer, in their interpre 


Vig =0 10 


tation of the results of impact tests 


leads to stresses ¢0,, 0a, 0 r.. that « OX s third 


: AvrHors’ CLOSURE 
derivatives ot ¢@ Phe conditions 
We wish to thank Dr. Wood for his useful com 


iin cases it will certainly be necessary to consider 


+ ndition n i the anvil and testing arrangement. In the 
Im pose he onditions o 


linal waves, this can be readily handled by the method of ¢} 
cteristics, but the analysis will be more complir 

test arrangement discussed in 

with elastic strains, were unde 

placements of the order 0.1 

the give 

0.001 in., 
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DISCUSSION 


ystem, and 80 It is per 


parts « l 
anvil to hay igid and at rest. This assumption 
is used 


e with t material which 


istic-plastic body For a 
analysis at small n magnitudes it would be 


vdded 


we for the elastix component of 


proximate satisi ac 
necessary to re 


move these assumptions and accept the complexity 


ive a negligible influence under the 


trair rge compared with elasti 


On the Inextensional Theory of 
Deformation of a Right Circular 
Cylindrical Shell’ 


N.d to the plane 


pendicular to the 


lore For loads symmetric which is pe 


the cylinder and divides the cylinder in 


1X18 OF 
two equal halves, the inextensional shell theory discussed in the 
a ring theory Timoshenko’s® for- 


paper reduces essentially to 


mula becomes the well-known expression for the shortening of 


the diameter of a circular ring subjected to equal and opposite 
radial loads applied at the end points of the same diameter if 
2DL is replaced by EJ. As 2L is the length of the shell and D 

Eh*/12(1 — v*), the difference between the two formulas is only 
the factor (1 — v*), that is, about 9 per cent For very short 
the shortening of the diameter of the shell 


calculated from the shell theory may 


cylinders, therefore 


he about 9 per cent too small 


consequence of the assumption that the anticlastic curvature 


s suppressed 


cylinders, however, the inextensional theory 


be seen from the following 


For very long 


must be incorrect, as can irgument: 


If A is the shortening of the diameter of a cylinder whose length 


to diameter ratio is 50, the shortening of a eviinder twice as long 


ind having the same diameter is A/2 from the inextensional shell 


theory. On the other hand, any material added beyond a dis- 


ince of 25 times the diameter cannot change the state of defor- 


mations in view of Saint Venant’stheorem. Theformula mentioned 


vields a shortening of the diameter which decreases to zero as 


the length of the cylinder increases indefinitely while in reality 


the diameter must shorten a finite amount even when the cylinder 


Yuan’'s* 


for concentrated loads, 


is infinitely long. This is borne out by investigation 


This argument naturally holds only 
ind the inextensional theory gives good results when the cylinder 
] 


is loaded uniformly along generators 


F. K.G writer does not agree with author's 


introductory 


OpQvIst The 
general 


statement 


as to the application of the 


theory of bending of cylindrical shells, if this is intended to in- 


> 


clude the type ol loading covered by Figs 1 to 3 of the paper. 
The 
point of view of the general theory for a pipe of infinite length 


For pipes of finite length and closed ends, corresponding and even 


vuthor himself quotes Yuan who treats Fig. 1 from the 


more general problems have been treated by the writer in a series 
of publications 
1, as well as the case of a simply 


1} Covering the case of Fig 


1951, issue of the 


ASME, ve 73, pp. 341 
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Head of Department 
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Ss. W. Yuan, Quarterly of Applied Mathematics, vol. 4, 1946, p. 13 
Professor 1 Kungl, Tekniska Hégskolan, Stockhol: 


of Aeronautical Engineering and Applied 
Institute of Brooklyn, Brookly N. Y. 


Swede 


gen 


support concentrated load 


(theory am I it rene mn * load along a 
generatrix 

2) Giving utlIon “along 
ratrix of the 


3] Giving a table numerical values of the transcenden- 


unctions concerned 


AUTHOR s CLOSURE 


The foregoing comments by Professor Hoff will lead t 
understanding of the interrelationship between th 
phenomenon involved in the system which is under consi 
ind the mathematics whereby we choose to analyze this physi 
It is indeed fortunate that the 


phenomenon discussion of the 


printed papers offers an opportunity for this interrelationship to 
be pointed out since the original paper 


m the 


whose object was solel 


to report result of experimental investigations, did 


some 


not admit the inclus such a discussion 


The additional bibliography supplied by Professor Odqvist is 


welcome emphasize, however, that the 


deal with evlinders with closed ends rather ian 
leal th lind th i i th thar 


paper 
papers ¢l 


those, nan , cylinders with open ends, which constituted the 


il paper 


the origi 


A New Method of Calculation of 
Repeat Factors for Turbines 


opic ol 


and Compressors' 


D. J 


heat factor represents an important item in the flow-path and 


BLooMBERG.? In steam-turbine design practice the re- 
performance calculations. With pressure ratios of 1000 or more 
this reheat factor can account for 20 or more per cent of the 
adiabatic available energy Each design group has its own way 
of deriving these reheat values—mostly by taking step-by-step 
energy values from a “‘state line’ drawn on the entropy chart 

In actual turbine-flow paths, the first stage derives no benefit 
reheat. The 


amounts of expansion and efficiencies 


at all from remaining stages will have varying 
For precise ¢ ilculations 


the step-by-step method is required 
usually are few 


I] pressure ratios are comparatively small 


In turbines utilizing air or other gases, there 
stages Also, the over- 
In these ca 


relative 


ses the reheat factor is usually neglected because its 


magnitude inconsequential Ar gain in energy 


resulting from reheat is usually offset by radiation and other heat 
losses 


The proposed method of calculation (which of necessity is 
limited to the 


stages operating 


situation in which we have an infinite number of 


it constant efficiency) will be useful mainly in 


preparing comparison graphs involving different gases 


While 


among the steam-turbine 


STEPANOF! the use of reheat factor is rather 


A. J 


general designers and further refine- 
*“The Load C: Action of Thin Cylindrical Shells 
The Royal Swedish Acader 
of Engineering Sciences, Proceedings No. 164, 1942, pp. 53 

™“Action of Forces and Moments l 
Along a Generatrix of a Thin Cylindrical Shell,” Jounnat or Appiiep 
MecHANics ASME, vol. 68, 1946, p. A-106 

5 HAllfe . rs Strength of Materials 

“T Kultur Edit Stockholm, Sweden, 1948, pp 

By Joseph Kaye Wadleigh, published in 
1951 ApPLinp Mecuanics s 
vol. 73, pp. 387-392 
Aircraft Gas Turbine Engineering Division, General 
Mem. ASMI 
Department ngers 


ASMI 


written 
in Swedish with sur ary in English 


Syvmmetrically Distributed 


written Swedish 
December 


ASME 


issue of Tue Jour L OF 

? Engineer 
Electric Company, Schenectady, N.Y 
ameron Ex 


Rand 


Engineer, ( 
Company, Phillipsburg, N. J 


gineering 


Men 








232 


ments are introduced in the method of calculation of the reheat 
factor, very well exemplified by the present paper, there is a 
definite trend in thinking, putting a somewhat different aspect 
upon the phvsical meaning and use of the reheat factor. 

In the field of turbocompressors, need for a reheat factor 
never has existed for reasons that follow. It should be realized 
that efficiency of a turbomachine (turbine or compressor) is a 
fixed definite quantity depending entirely upon the degree of 
and mechanical perfection. Like that for 


hydrodynami any 


machinery it is equal to 
Efficiency = output/input 
The output of a turbine can be measured easily but the input 


cannot be ascertained directly This led to the introduction of 

adiabatic” efficiency based on the isentropic enthalpy drop for 
the input. However, the true input of the available energy is 
lower than the enthalpy drop; hence the true output/input effi- 
“adiabatic” conventional efficiency, 


ciency is higher than the 


the difference increasing with the expansion ratio. Therefore a 
correction factor’ became necessary to connect the actual out- 
put with the isentropic input and known adiabatic stage effi- 
ciency 

It has been shown that the true input to a turbine is repre- 
sented by the change oF available energy across the turbine, or 
The latter is 


Keenan’s term for a special case of Gibbs ‘‘zeta”’ function, widely 


what became known as change in “availability.” 
used in thermodynamics of chemical processes under different 
Efficiency of a multistage turbine based on availabil- 
“effec- 
and remains constant for the same individual stage 


names 
ity is a true output/input efficiency, called by Keenan 
tiveness” 
efficiency, thus requiring no reheat-factor correction 

Statements that efficiency of a multistage turbine improves 
due to “recovery” of heat resulting from hydraulic losses is most 
misleading and has all the marks of infringement on the second 
law of thermodynamics 

With turbocompressors the input can be measured directly 
but the output (foot-pounds per pound, or head in feet) is cal- 
culated from the observed pressure and temperature rise by using 
theoretical adiabatic formulas, for uncooled compressors The 
idiabatic conventional efficiency is lower than the true output 
input efficiency, because the true output (bead in feet corres- 
is higher as a result of temperature 


The true output /input effi- 


ponding to a given pressure ) 


caused by hydraulic losses 


rise 
ciency, sometimes referred to as “internal hydraulic” is closely 
upproached by using “polytropic”’ formulas for the output cal- 
The polytropic n for gases obeying the 


= RT is easily calculated from the observed pres- 


culation exponent 


state law p 
sures and temperatures. For new machines the polytropic ex- 
in estimated polytropic efficiency 
using known relationships from thermodynamics. The poly- 
trop (ce ntrifugal 1s 
ower than that of the individual stages but for hydraulic rather 
The multiplier for the adiabatic 


ponent is determined from 


efficiency of a multistage compressor 


than thermodynamic reasons 
efficiency to obtain the true output/input (polytropic) efficiency 
than unity and not less than unity as the “reheat fac- 
This is easily 


is greater 


tor reasoning in reverse’ would indicate cal- 


culated with available formulas in terms of polytropic efficiency 
ind the known gas properties 
Polytropic efficiency as a best approach to the output /input 


efficiency applies with equal accuracy to gas turbines; thus in 


this field use of reheat factor is hardly justified 
trends in theoretical thermody- 


for 


It is not expected that new 


namics will change the design methods now in use steam 


turbines, but it is important to bring the theoretical reasoning 
behind them up to date 
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AuTHORS’ CLOSURE 


Mr. Bloomberg's comments are appreciated by the authors 
It should be noted that the paper was intended to review the 
laborious step-by-step methods of computation of reheat factors 
for steam, as noted by Mr. Bloomberg. It was also intended to 
present a simple method of computation of reheat factors for all 
gases for which the equation of state, py = RT, is valid, without 
the additional requirement of constant specific heat 

Mr. Bloomberg indicates that the proposed method of calcu- 
lation is limited to the case of an infinite number of stages in the 
turbine. In the paper we have shown, however, how to handle 
In fact, the 
proposed method has demonstrated the customary rule-of-thuml 
the 


that the ratio of the reheat gain for a finite number of stages to 


the case of a finite number of stages in the turbine 


used as an approximation by turbine designer, namely 


the reheat gain for an infinite number of stages is approximate! 
equal to one minus the reciprocal of the number of stages, or 


1 
Pins ec n 


Although we agree with Mr. Bloomberg's opinion that the de- 
sign of present-day gas turbines may not justify the refinement of 
the reheat factor, it is probable that future gas-turbine designs 
may find it highly useful. 

Dr. Stepanoff’s remarks on reheat factors applied to turbo- 
compressors are interesting. It is evident that the utility of re- 
heat factors for design of turbocompressors depends on the 


The 


gain in energy 


particular definition of “efficiency which is selected 


authors agree with Dr Stepanoff that the phr iSes 
“recovery of heat due to hydraulic losses’’ are 


completely misleading and should never be used 


due to reheat” and 
He re again we 
feel that terms should be carefully defined to avoid the pitfalls of 


ambiguous reasoning 


The Safety Factor of an Elastic- 
Plastic Body in Plane Strain’ 


° 


F. 8S. Suaw.? 
problems of considerable interest related to the behavior of metal 


To the structural engineer one of the classes of 


structures is the so-called elastoplastic stress distribution At 
present, however, to solve such problems is generally impossible 
best piece of 


Accordingly, the authors have provided the next 


information, namely, bounds on the “‘failing load 

The authors’ work is a further extension of information in the 
field of limit design which is now receiving practical application 
In this instance, however, it is necessary to 


by the 


in structural work. 
initial as- 
First 


ideal pl is- 


keep in mind the severe restrictions imposed 
sumptions concerning the idealized materia] discussed 
(and this point is now well know) the assumption of 
ticity”’ behavior is, at best, a poor approximation to the behavior 
Then too 
importance, is the fact that for a load a little smaller in magnitude 


of most structural metals in use ind of considerable 


than the ultimate load that can be carried by the « omponent, no 


information other than a ratio is available showing how close to 


failure the component may be. Thus imperfections and varia- 
tions in material properties, and so on, cannot readily be taken 
into account. For design purposes, then, the theory should be 
used with considerable care, and supplemented by experimental 
information. 
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AuTHOoRS’ CLOSURE 


As Professor Shaw explains, idealizations always impose re- 
strictions on the practical use of theoretical results. However, 
the engineering limitations of a perfectly plastic solution are no 
more severe than those of a perfectly elastic solution and are often 
less In particular, any imperfections in the material in the 
sense of local] regions of low-vield strength ean be handled rela- 
tively easily because the theorems of the paper still apply 
Also, initial stresses, which may completely nullify an elastic solu- 
tion, have no influence on the plastic solution. The major dif- 
ficulty in limit analysis involves the question of whether or not 
the ductility 


of the material is sufficient for the solution to be 


valid. Much experimental work is needed on this problem but 
the answer will certainly be as clear as is the significance of elas- 


tic solutions in fatigue problems 


A Statistical Distribution Function 
of Wide Applicability' 


T.C. Tsu? 
vised a distribution function of truly wide applicability, as evi- 


The author should be congratulated for having de- 
denced by the seven examples presented in his paper 

Since the writer is currently concerned with the problems of 
particle-size distribution in aerosols, he is interested in the 


pos- 


sible utilization of the author’s method to reduce the necessary 


amount of experimental work. In this connection he would like 


to ask the following questions 

1 In applying the author’s distribution function it is necessary 
If the distribution 
function is a true representation of the observed data, then any 


to determine the parameters z,, ro, and m 
three sets of the values of P and x would be sufficient to evaluate 
In the author’s examples he did not 
Would he care to 


these three parameters 
specify how his parameters were obtained. 
discuss this point briefly? 

2 When the function is applied to an unknown distribution, 
how many observed data are necessary to yield the parameters 


? Considerable practical value would be added to the 


reliably‘ 
author’s function if its application could result in a saving of ex- 
perimental work 


3 The relations shown in Figs. 1, 2, and 3 in the paper, appear 


to represent the equation 


P 


rather than 


Could that be a misprint? The values for log (7 —7.) in Fig. 2 do 


not correspond to the given values of x and zu in the second ex- 
ample (size distribution of fly ash), the discrepancy being appre- 
Could there be 


If so, would the author kindly show a corrected figure 


ciable when x is small some numerical errors? 


R. A. MuGere The author's treatment is definitely 


a con- 
Extended Limit Design Criteria for Continuous Media,”’ by D 
C. Drucker, W. Prager, and H. J. Greenberg Applied 
V athematics, vol. 9, no. 4, January, 1952, pp 
1 By Waloddi Weibull 
the JouRNAL oF APPLIE 
293-297. 
Associate 
State College 
Oakland, Calif 


Quarterly « 
381-389 


published in the September, 1951, issue of 
Mecnanics, Trans. ASME, vol. 73, pp 


Professor of Engineering Research, The Pennsylvania 


Pa. Mem. ASME 


State College 


tribution to the literature on distribution functions. The range of 
fields treated in his examples is also impressive 
However, the reason for introducing the minimum value z 


and ignoring the maximum z,, is not entirely clear. Probably 
it relates to the original applications, which may have been the 
Cystoidea or the yield strengths and fatigue-life data of steels 
Now, for such a case as Fig. 2 of the paper, one would expect 
the maximum parti le to be more tangible , and also more s gnin- 
cant practically , than ther 


Incorporation of both 


nnimum 


a maximum and a minimum value of z 


into the form 
( z-—Zu ) 
= 
€ 


which will again reduce to Equation [5] as 


will bring Equation [5 


becomes infinite, 


ind to the Rosin-Rammiler type of equation‘ as 7, vanishes 


Of course one may start with any distribution function where 


the argument has infinite range, and convert it to one where the 


range is finite. This has been illustrated in the case of the log- 
norma! distribution by Van Uven® and more recently by Mugele 
The 


Rosin-Rammiler 


ind Evans.* atter reference also gives a critical review of the 


ind other distribution functions 


A word of warning also should be added in regard to the ex- 


amples in the paper: They contain some arithmetical and dimen- 


sional errors However, when these are corrected, the examples 


illustrate excellently the general statements of the text 


F. A. MceCuntrock 
the author is 


The distribution function suggested by 


attractive because of its simplicity, the ease with 
which it can be applied to studying the size effect, and its implica 
tion of a lower limit to a distribution In order to apply the cis 


tribution impartially, however, soine systematic means of fitting 


it to experimental data should be used. For a simple distribution 
the following proc edure appears us¢ ful 

The parameters, z,, z,, and m can be chosen so that the dirst 
three moments of the distribution function coincide with those of 
the data The nth 


first calculated from the cumulative distribution 


ink 


Differentiation gives the frequency distribution 


moment of the theoretical distribution is 


r 


F= 


The nth moment about z, i 


exp | 


On changing the variable of integration to 


Rosin and | 
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Rammler, Zeitschrift des Vereines deutscher Ingenieure, vol. 71 
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this becomes 


This integral can be expressed in terms of the Gamma function 


1 and third moments about the 


measure 


Since a; is a function of m only, the value o in be chosen so 


that the values of a; for the theoretical distribution and the ex- 


erimental data coincide. Then since the second moment about 


the mean, that is, the square of the standard deviation, of the 


experimental data is known, the relation 


n/z,2 = TL + 2/n 


be solved for z,. Finally, the relation 


Mi’/z, = (2 ‘ 3 t ‘ 9 


solved for z,, since the mean of the experimental data is 


Plots of the quantities ay, u ind wy’/z, are given 


La 


ind 2 of this discussion 


ric 


rhe writer would like to ask what procedure, preferably system- 
itic, should be followed in the case of a complex” distribution 
An extension of the looks impractical, 


et the 


foregoing procedure and 
writer would like to try applying the distribution in other 


would be 


cases For example it interesting to see whether the 


other data on the ST-37 steel reported by Miiller-Stock would re- 
sult in the same division of the population as found from Figs. 6 
ind 7 of the paper 


APPLIED 


MECHANICS 


AUTHOR'S CLOSURE 


The 


Che proposal of Professor Tsu to take any 


author appreciates the comments made |! discussers 


three sets of the values 
P and z is quite correct but does not use the 


rhis method may be 


data efficiently 
improved by taking the set from a smoothed 
curve. Up to the past year the author’s usual method has been 
to plot the data as shown in the paper and to choose the value 
In this way 


the distribution is simple or complex, but 


u 


O give the best straight line 


to ce 
entirely free of subjectiveness 

About a year ago the author decided 
the 


ris the mean and o the star 


start by standardizing variable 


oa, where 
nating two ol the parameters 
hen takes t 


ior insta 


tribution function t 
exp 


\ curve paper for different values of a 
ardized Gaussian distribution, may be prepared B 
the points (P, z on this paper, it is easy to decide 


whether the dis 


tribution is simple or complex and to estimate, with a good ap 


proximation, the value of @ 

As to the third question, the parentheses are an awkward mis 
print. The values for log 
y= 1.5 > 
the 


been inere 


in Fig. 2d 


20u but to x, =1 


o not correspond t 
It should be 


values and 


the given value 


that 


20 


mentioned r-values are mid-point should 


correctly 


the ¢« 


lave ised by Thus the value 0 u is 


orrect one 


The introduction of a maximum value 


Mu 


not found 


m proposed by Mr 
function. It 
necessary to introduce this new parameter it 

of materials, probably because the theoretical 


the techni 


gele is a valuable extension of the was 
ld of strengt! 
strength n 
perhaps a hundred times higher thar 
iy be quite 
Professor McClintock 
three moments is quite good if the distribution is sir 
The 
f computing the parametersand has mentioned it (witl 
the 


however, he ver ap} 


in other fields conditions m lifferent 


ere 


The method proposed by 


population not too small iuthor has been aware 
possibility 
some different notation 
Actually, 
but admits that it may sometimes have its advant 


As to the 


lor gamma function) in an earlier 
paper.* vied this method 


question of a systematic procedure when the distri 


bution is complex, the author is sorry imit tl so far he has 


found no better method than to cut course 


ichin 


comput i 


not very satisfactory, but a simple el 


recent omplet facilitates the Us 
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Hao? aND K. J. Triccer.*? In applying the new 
analyzing stress and strain distributions during chip 


the higt 


machining operations the metal behaves lil ul illy plastic 


formation iuthors assume tl 
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DISCUSSION 

material that does not work-harden. Experience obtained by the also upon the effective time of heating. It is difficult to analyze 
writers in the actual cutting of three kinds of steel and one typeof this aspect of the built-up nose solution without further experi 
electrolytic copper fully justifies such an assumption. In fact, mentation 
this constitutes one of the major conclusions reached by one of the Quantitative verification of the authors’ new theory needs an 
writers 4 vears ago after conducting research in metal cutting in accurate determination of the tool-chip friction angle A, which is 
the Machine Too! Laboratory of the University of Manchester.‘ not the same as that determined under the condition when shear 
A ductile material undergoing high-speed chip formation acts as occurs over a narrow band. It cannot readily be calculated from 
f no work-hardening occurred during the time that deformation — tool-force dynamometer data and chip thickness measurements 
takes place along the shear surface as indicated by line AC in Fig Examination of the nature of tool-chip contact under a magni 
1 of the paper However, this does not mean that the deformed fication of X30 to 50 after the tool has been in use ft sho 
metal (chip) is not strain-hardened. It was first shown by E. G ne reveals the « ‘ f two distinct regions 
Herbert that the hardness of the chip may be several times that of wctive tool face This may serve as a clue for 
the undeformed chip material. That such a paradox can be ex- mental determination of tl ngle @ by actually measuri: 
plained by the dislocation theory for the plastic deformation of length EC in Fig. 6 6 is known, together with the tool 
crystalline solids has been given in a recent paper force and chip-thickness data, the combined Mohr circle diagra: 

It should be pointed out that the yield stress & referred to by the Fig. 8 of the paper, « drawn without difficulty 
authors is the dynamic yield stress. It differs appreciably in It is hoped that iuthors will continue their work in the 
magnitude from either Tresca’s yield limit or von Mises’ quad- future 
ratic limit. Both these quantities refer to quasi-static loading 
conditions. For structural] steels, it has been found that the dy- R. S. Hann Phis paper presents an interesting analysis 
namic vield stress is in the neighborhood of the static ultimate the metal-cutting process and seems to be the first to give 
tensile strength analysis involving the built-up edg Previously it has been cor 

The authors indicate that the influence of the compressive stress sidered that the type 2 chip was formed without a built-up edge 
on the dynamic vield stress is much less than that assumed by The writer, in tests where the chip appeared to be of type 2, has 
Merchant. Recently this has been brought out in the cutting in found a very small built-up edge in all cases—indeed a few thou 
vestigation of spheroidized SAE 52100 steel using a wide range of — sandth of an inch in size Chips that are produced by an abrupt 
speeds, feeds, and three different tool-rake angles. No corre- ending caused by the tool entering an interruption often carry 
lation has been found to exist between the so-called machining — their end a very small built-up edge which can be seen under the 
constant and the “assumed” slope of yield stress-compressive — microscope Such chips, even though their back side is ver 
stress curve Hlowever, Merchant’s plasticity equation does smooth and highly burnished, have been formed with a built-up 
gree with experimental data to a good degree of approximation — edge The region BC in Fig. 6 of the paper serves to burnish the 


chip so highly that all roughness caused by the flow along FC is 


kK¢ 


friction angle at . rac ted. and general] t has b | “ suc . 
© (hear ansie) 4 ( on a es eradicated, and generally it has been believed that such chip 
{ 


ool chip interface were produced without a built-up edge Consequently it 
re! ; : appears that there is considerable truth in the authors’ conelu 
rhe authors speculation of the existence of a tiny built-up edge sion that 0 built-up edee always existe 
which arises natural] is a consequence ol the mathematical anal- The analysis presented considers flow to take F lace in the 
ysis may serve as an explanation for the foregoing controversy 


] 


: direction only, Fig. 6 That this is not strictly true is evidenced 
and should be explored fully by future experimentation 


; ”y the curvature of the chip rhe authors state that urva- 
k:xcept for the small percentage which is retained as latent heat 
in the deformed chip, nearly all of the energy expended in metal- Remnants | 


ngineer ‘ i Machine Compan 
machining operations transforms into sensible heat, resulting in a Mass. Mem. ASMI 





large temperature rise The effect of such temperature rise on 
the mechanism of chip formation has been reported recently.** 
In distinction to what the authors have stated, the dynamic 
vield stress is practically unaffected, if shearing takes place in a 
narrow zone, since the heating effect does not actually occur until 
ifter the shear energy has been imparted to cause plastic deforma- 
tion. On the other hand, such heating due to deformation and 
that produced due to tool-chip rubbing do have an appreciable 
effect on the mechanism of chip formation through their influence 
on interface friction 

With a slip-line field configuration as shown in Fig. 6 of the 
paper, the deformation process is greatly complicated. Shearing 
takes place over a fan-shaped region AF( Both temperature 
ind strain-hardening may exhibit their influence, depending not 


only on the magnitude of temperature and strain involved but 
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ture of the chip is due to residual stress and thermal] strain con- 


ditions. The fact that curved chips may be heated to incandes- 
without changing curvature tends to disprove the fore- 
Also, they will be found to retain their curvature to the 


dissolved in acid 


cence 
going 
last as they are It appears from experiments 
that the chip curvature is a result of flow along the u-slip lines as 
ertain cases the chip velocity of the 


metal adjacent the tool face has been observed to be as much as 


well as the v-slip lines. In 
oO per cent higher than the velocity of the metal forming the out- 
side face of the chip. We may consider the grains in the uncut 
metal to be essentially equiaxed so that after deformation the 
become elongated showing the directions of principal strains. It 
1 of this discussion that a nonhomo- 


and that the 


vill be evident from Fig 


geneous state of strain exists authors’ analysis wil! 


not apply accurate ly in such cases 


M. KuGene 


nt step forward in the 
ing. However, while the present-day plasticity theory used by 


the authors is sound and well establishe 


MERCHANT This paper marks another impor- 


le velopment of the science of metal cut- 


, it appears that the prac- 
tical assumptions made in applying it to the metal-cutting prob- 
The result is that 


the resulting limiting chip-stress solution 


m are open to some question while the 


jualitative findings of 


igree well with observed facts and in some cases offer a better ex- 


planation of these facts than had previously existed, nevertheless 
the quantitative predictions of that solution do not agree as well 
th experimental data as do the predictions of the more approxi- 


minimum-energy solution 


The first case which the authors consider is that of machining 


without a built-up edge. In applying the theory to this case 


they assume that the yet it is known 


vield stress is a constant, 


iffected by a number of factors Among the 


these 


that this stress Is 


is that of normal] stress (hydrostatic 


mportant of 


the effect of which has been demonstrated by Bridg- 


1 Rotner™ (the writer is indebted to Prof. M. C. Shaw of 
latter 


for his kindness in bringing the reference to his 


Measurements made in our own laboratory by 


ind Whisler 
t Data obtained in these in- 


ittention 


meny ind by Krabacher ? on a Bridgman-type 
varatus confirm this same effe 
tigations, as presented by Ernst'* are shown in Fig. 2 of this 
liseussion. Introduction of this effect of normal stress on yield 


ss into the minimum-energy solution has a profound effect 
upon that solution, making it a much better approximation to 
Therefore is it not likely that the introduction of this 


same effect into the authors’ limiting chip-stress solution would 


reality 


ilso have a profound effect on the quantitative results? 

The authors consider next the case of machining with a built-up 
edge. The fact that the authors’ limiting chip-stress solution pre- 
dicts the presence of a built-up edge marks a real step forward in 
Little 


developing the mechanics of this case previously 


the mechanics of cutting progress has been made in 
However, the 
quantitative predictions of the theory as to the size of the built-up 
edge, the values of friction at which it comes into existence, and so 
Pri- 


on, are questionable because of the assumptions made 


* Assistant Director of Research, The Cincinnati Milling Machine 
Company, Cincinnati, Ohio. Mem. ASME 

*See Authors’ Bibliography, reference (19 

“Change of Mechanical Properties of Metals Under Hydrostatic 
Pressure,”’ by 8. I. Rotner, Journal of Technical Physics, USSR, vol 
19, 3rd edition, March, 1949, p. 408 

't**Torsion and Compression Testing,”’ by J 
of Cincinnati thesis, 1947 

'2*Torsion and Compression Testing’ by E. J. Krabacher and 
K. W. Whisler, University of Cincinnati thesis, 1949. 

*“Fundamental Aspects of Metal Cutting and Cutting Fluid 
Action,”” by Hans Ernst, Annals of the New York Academy of 
Sciences, series II, vol. 53, 1951, pp. 805-823 
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vield stress again ix assumed to be « constant and inde- 


marily, 
pendent of the several factors which are known to influence it 
This assumption alone could have a profound effect on the quanti- 
tative results of the theory as already indicated Then, too, the 
authors assume that the bottom of the built-up edge is in fric- 
tional contact with the work surface, whereas photomicrographs 
of chip formation have shown that the built-up edge actually is 
it this point. Continuous 


plastic flow takes place in this region in the s 


continuous with the work material 


ume manner as it does 


on the shear plane and strain-hardening also occurs. As a result 
the bottom portion of the built-up edge gradually builds down- 
ward, increasing in size, until it becomes so large unstable 


ind is then carried away with the machined surface, after which 


the cycle begins again. (It is this process of the building up and 


sloughing off of fragments of built-up edge on the machined sur- 


face which accounts for the roughness of su faces when pro- 


duced under condition where a built-up edge 
In the composite solution, resulting fro ition oO 
presented a 
limiting cl 
They in 


yond the 


the two cases already ads ussed, the 


comparison of the quantitative predicti Llp 


stress solution with experiment in Fig. 1 
dicate that the reason the 


theoretical boundary det 


imental 
e figure is that 


elastic deformation modifies the l-chip int 


solution he toc 


face If this is so, then it would appear th ita should show 


at least some downward trend towar 
the boundary, in keeping with the do 
ary However, no such downward tret 
nor even in other data obtained at higher 
shown experimentally in Fig. 12 of the pa 
+} 


ugh t 


the minimum-energy solutior 
theoretical point of view, actua upproxir 
the « xperiment al results than miting chi 
stress solution in Fig. 12. The« ation giver 
the minimum energy solut 

propose: 


cussion. (The symbols use 


the writer.“ but the ordina n ’ ire the same 
those of the authors’ Fig 

It can us 
mental data to the line labeled 26 + 


the minimum-energy solut i is 


ation of the exp 


12 
be seen that the degre 


77 deg predicted b 
better than the approxi- 


iting chip-str solution in Fig. 12 


t to this 


mation predicted by the li 
t since it involves 


However, the authors obje« eemen 


the proportionality constant relating vield 


the conclusion that 
stress to normal stress must be of the order of 0.23 (cot 77 deg 
since according to the minimum-energy solution the sum 26 + rf 
a should equal approximately the comp! nt of the slope 
angle of the vield stress versus normal stress line in a plot such as 
that shown in Fig. 2 of this discussior ertheless, Bridgman 
Rotner,” Kemeny," and Whisler'? all 
values of this order for the proportionalit for a variet 
n the values of 26 4 


and values of the con 


Krabacher obtair 


of metals. For example, a con 


parison be 
Tt a, obtained from metal-cutting 


plement of the slope angle of the vie obtained fron 


torsion-compression tests I rade in this laboratory on four different 


steels, is given in Table 1 It « ” seen that the degree of cor 


relation between the data obtained from the two quite different 


types of tests is relatively good i appear therefore that 


for the values of strain common in utting the influence of 


hvdrostatie pressure on y ield stress is of a larg nough magnitude 
to be significant in the theory of machining 

Since the authors’ limiting chip-stress solution based on plasti 
actually results in a poorer 


itv theory ipproximation to exper 


ment than does the simple and approximate minimum-energ 


'*See Authors’ Bibliography, reference (18 
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TORSION & COMPRE SSION 
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COMPRESSIVE STRESS 
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2 30)~C« 50 
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NE 9445 STEEL 


SHEAR ANGLE, DEGREES 


TABLE 1 rTING AND 


COMPARATIVE RESULTS OF METAI { 
TORSION-COMPRESSION TEST DATA 


f slope 
i-stress 
~ torsion-com 
Steel wearit pression tests 
SAE 3115 : >. 87 78 
SAE 3150 
spheroidized 
SAE 3150 
pearlitic 
SAE 3450 


solution, even though on a theoretical basis it should be expected 
to be more exact, it would seem that the next desirable step for 
the authors to take 
assumptions into their solution to determine the extent to which 


would be the introduction of improved 


Fur- 
thermore, it would appear that the authors would benefit by 


these help to bring it into line with experimental findings. 


carrying out a few critical experiments to test certain points In 


their theory as a guide in choosing the assumptions to be used 


It is the writer's feeling that if these things are done a useful and 
practical theory can be developed from this approach which will 
mark a major advance in the science of machining 

ki. K. Henriksen. Metal-cutting operations have been in 


vestigated more or less continuously since the 1870's and, there- 


fore, it is not surprising that it is hard to find basically new 
thoughts in the constant flow of contemporary public ations on this 

The present, paper, however, is a very pleasant excep- 
this The 
by original methods, and they have laid a new cornerstone to the 


subject 
tion to rule. authors have approached the problem 
foundation of a great, vet unfinished building. 

However, many engineers and other workers in the metal-cut- 
ting field will fail to visualize the value and importance of this 
paper, simply because they do not feel themselves prepared to 
follow the this 
writer would like to present a simplified development, leading up 


mathematical methods used. For reason the 


to the principal equation 


Referring to Fig. 4 of this discussion, AC is assumed to be the 


plane where the shearing deformation occurs, and it follows that it 


is the plane of maximum shearing stress. B is the upper limit for 


ontact between chip and tool and, therefore, the upper limit 


or transmission of force between the two bodies It follows then 


There 


There is no stress on AB 


is no load on AB 


Specifically there is no shearing stre 


Hence AB is a plane of principal stress 


From the general rule that the maximum shear occurs under 45 


deg with the planes of principal! stress, it follows that 


XBAC = 45 deg 


With the shearing stress & and the normal stress p, it can be showr 
elementarily by means of a triangular free-body diagram p = / 

Now the shear stress is contant, equal to k, over Af it Is rea- 
sonable to assume that also the normal stress is constant, and 
equal to p But then the resultant R will go through the middlk 
of AC an angle of 45 deg and therefore, be parallel to AB 
and its point of intersection E with BC will be the middle of Bf 
Equilibrium of the body ABC 
forces on BC will fall in the same line, that is, Pass through E, but 
as E is the middle of BC it is reasonable 
to assume that the load on BC is contant 
Now 


, under 
requires that the resultant of the 


so far, at least, see later 


t{ DEC = —hry 


where A, is the friction angle between chip and tool, hence (from 


the triangles ABC and DEC) 


* Professor in Charge 
essing, Cornell University 


Materials Proc 


Head of Department of 
: ASMI 


Ithaca, N. ¥ Mem 
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velow indicates that the consensus of opinion on this contr 
versial issue appears to be against an appreciable hydrostati 
pressure effect 

The study of plasticity of metals has led to the general 


onclusion that while fracture conditions are highly influence 


hich is the bas quatic ; t Apel 4 . is ele- . 
— sic Equation : is ele hydrostatic pressure, this does not affect plastic-flow cor 
mentary an 

renta : , ' tions appreciably This is emphasised by the general accept 
s to t subject itself, th thors yased t ry 8 ‘ o- 
A ne Can aitecaphesandcsieaiee ave based their analysis ince of the Mises or Tresca yield conditions, both for initi 
upon the issumption of no load beyond point B, that is, ignoring ld. and 7 f rali i str train r tion 

ela, 2nd ast isis lor & generaize stress-strain reiath ove! 

stic stresses in the cl ‘or analysis of this character s : 
elasti . ves aa For aed . his character pens inge of strain. Both these limits are independent of the average 
ssumptior s nature anc ustihed, t consequences of this , } , 
. pre ‘ ms =e a : a drostatic pressure, depending only on the stress deviator The 
issumption being that AB is a straight line, and constant pressure 


over BC 


nsensitiveness of the flow stress in tension to hydrostatic pres 


ee sure up to large strains has been che¢ ked directly by Bridgeman 
he writer would like to show the effe of omitting this limit- . , 

rhis is in agreement with the concept of plastic potential (se« 
Ing assumMptior ° 
= ; Hill’) which has recently received both experimental®® and 


il?! justification. Thus we must answer the questio: 


In metal forming theory why is it generally accepted that hydro 
if 


tatic pressure affects fracture but not plastic flow, and ye 


whining theory an appreciable influence on plastic flow has 
een commonly accepted? A recent discussion of this questior 
connection with a paper by Chao and Bisacre™* shows that 
inv workers in the field of metal-cutting theory do not accept 
ble influence of hydrostatic pressure Discuss¢ 
M. C. Shaw, and bE. H, Lee independently empha 
Shaw states: “Calculations based upon atomi: 
tructure and inter-ionic forces revealed that normal stresses of 
ibout the magnitude that were found in metal cutting should have 
Fig. & no influence upon flow stress That the yield stress for ductile 
etals was found to be the same in tension and compression was 
there will be a region FB, Fig. 5 of this dis- further evidence that normal stress did not influence the ductile 
y tic stresses in the chip and there will be a sharp properties of metals. Schmidt and Boas had demonstrated tl 
s¢ in pressure from F to B. The plastic deformation begins at B_ normal stress did not affect the flow stress of ductile 
beyond 8B the rate of increase in pressure will be less, the pres- —« In the closure Chao and Bisacre agreed 
1aximum and will then decrease somewhat to- hav comments and mentione 1 independent che« 
Bridgeman’s compression-torsion-t yp st with a@ brass 
boundary AB between the elastic and the plastic region in = gave a zero hydrostatic-pressure effe n view of these 
f no longer will be a straight line, but will be curved ments, what then is the ificance of yapers quoted by Dr 
» side in the direction of the flow of the chip. Merchant? In the report,** of I per under discussio 
is a shortened version, the shortcomir f the notched torsior 
AuTHORS’ CLOSURE , 
compression te wi ‘ is test serves ad 
withors first wish to express their thanks for the interest I i o demonstt 1e inf of hvdre itic pressure 
hown | the discussers of this paper intuabit tr ‘ " t 1 tral ribution ma 
Chey appreciate the comments of Professors Chao and Trigger erp on i rms of flow stre ctremely difficult r} 
onfirming the use of an ideally plastic stress-strain relation, and est do not product under hydrostatic pressure 
emphasizing the smaller influence of hydrostatic pressure than ppreciable longitudir 
had been assumed in earlier work. The discussion of the in- 
fluence of temperature distribution will be important in the future 
development of the theory; in the analysis under discussion such 
influences were merely averaged and results compared which 
vould not be expected to be sensitive to such averaging 
We appreciate the comments of Dr. Hahn concerning the ex- 
perimental observation of a built-up edge as predicted in the 
paper We ag that a more detailed study of chip strains in- 
cts is called for. In this connection B. W 
SI as d oped a theory of the turning process with rela- 
he tool and the work, in contrast to the 
under discussion In this solution varving 
chip thickness occurs 
discussion is concer d mainly with the in- 
stress in shear. In 





Dis 


his 


we 


this test 


of the 


Rotner’s work 


precedence over results in Because < 


language difficulty were unable to examine 
in detail in the short time available 
Thus it seems that the 


Merchant’s theory inc 


agreement of experimental 
effect ca 
This point was made by J 
The 


stress solutions given in the paper under discussion offer 


results with 
lls for explana- 
M 


limiting 


luding a pressure 


tion on other grounds Lickley 


in the discussion mentioned previously chip 


such an 
explanation 

The 
Merchant correspond 


discussion of edge development given by Dr 


» intermittent nose usually classed 


is a type-3 chip lap mall built-up nose discussed in 


the paper behaves differently according to its connection with the 


base material lowever, it would be of interest to study the 
nfluence of material 


We do not 
the point F in Fig 
It 
smooth extension of the 
pw, = 1 
We agree that the 


under the 
the need 


shear nose base on the present 


tendency toward 
Merchant 


very well merely 


theory ippreciate¢ lor a 
12 of the 


vould seem that elast 


paper, as suggested by Dr 


s might involve 


constant beyond the bound- 


H 
In the ec 
it seems that it may not be possible 
We 


e of considering the bending effect in 


Professor n- 


simple derivation given by 


riksen may have greater appeal to some engineers ist 


of the built-up nose solution 
to avoid the mathematical considerations in this way 


with the 


agree 


importance i more 


} 
ilvsis 


detailed ar 


The Calculated Performance of Dy- 
namically Loaded Sleeve 
Bearings—IIT' 


experimental data whicl 


the 


loud case 


is 
short-bearing 
The 


Ocvirk is listed in 


to support 
the steady 
Professor 
12). 


W. Oecvirk and the writer are associated with a 


approximate me thod of solution for 


theoretical part of our project b 
the author's Bibliography as reference 

Prof. I 
sear hp ojyect at ( 


National 


project 


re- 
ornell University being sponsored and financed 
Advisory ( This 


was begun in it the present time is 


the ommittee for Aeronautics 


1948, and 


experimental investigation of eccentricity 


LAS 
reste ircl 
engaged in an ratio, 
friction, and oil flow to obtain actual data on this subject for 
The 


oad > Was 


the steady-load short-bearing approximation theory 


f 


ist 


or the steady extended recently by Professor 


Oecevirk and is being used as a convenient vehicle for plotting 


| Both Professor Oevirk and the 


the experimental results writer 
would like to offer some brief comments 


We shaft 
on projected area 900 psi, 


is 


n- 


at speeds up to 6000 rpm, loac 
with SAE 10 oi The ecce 


tricity of the shaft in the bearing is being measured by a system 


are using a 1*/,-in 


up to 
ilso are obtaining friction and 
At the present time 
the data cover length-diameter ratios of > ind | 
the ring ethod 


lines can be plotted to represent the eccentricity 1 


of levers and dial gages, and we 


oil-flow data in both directions of rotation 


Using short-he approximate mn theoretical 


illo, m versus 
the Sommerfeld number 8, and we obtain 
i of 1, 


ia 


fan curves, 
On 


a group of 


on which the lines for and be labeled 


an 


such a plot these lines are far enough apart for the l/d mentioned 


so that experimental data should show a trend quite readily 


1 By J 
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? Professor of Mech 
N. ¥ 


Deve 
ASMI 


her 


hed in the 
Mecnanics, Trans 


T, Burwell, publis 


APPLIED 


anical En 


SSION 


{ 


general agreement 
above 


u\ 


points show a slightly 


the theoretical lor the st 


load case will note within a 


few months 
Acc 


heoret 


ording 
1 


aking 
ne The 


moderate 


use 
expermne! 


spread slight num 


ber, originated by our pro 


12 
} 


turn out to be ol n 


author just 
We be 


practi 


number 


above his Equation 


may re | value -known 


f 
gs ol 


Sommerfeld number for bearu a short 
a litth 

These 
the 
profile, and 
which ean be 
ever 


Th 


flow 


theoretical curves are approxin n t include 


by a 


the 


illustrated 


ol 


which can be 
they 


illustrated | 


side flow parabolic velocit 


circumferential 


flow 
Hlow 
tri ingle 


the 


include art 


the | 
vy a triangular velocity profile 


the 


superimposed 
hand 
1 includes both parts of the circumferential flow 


the 


theory should be considered as the two halves of the same ty 


parabola on 
sick 


Thus 


uring 


Sommerfeld solution, on the other neglects 


ire 


it that the Sommerfeld theory and short-be 


, 
pe of 


solution, rather than as be ompetitive ach has a range of 


the 


deserves 


ing 


d in which it more useful In our opinion, this short 


ing theory wide recognition for its ability to in 


lude oil flow and give a useful approximation for short-bearing 


lengths commonly used 

The 
the 
others take 180 deg, neglecting the negative-pressure region 
that the the 


th 


in extent 


author also points out that some of solutions assume 


while 
He 


300-deg 


oi! film for load purposes to have of 360 deg 


ilso points out attitude is WO deg tor 


t ingle 


semicircle for the 1S0-deg or 


OO psi 


or 27 case, and ay 


yproximates a 


wT « Our experimental data for unit loading up to 


ase 
tend to follow the semicircle result 

experimental work reveals a number of 
s shaft deflection, v 
which 
importance ol 
This will be 


In general, the pr 


iriations in ¢ 


affect 


onsiderations such 
n oil-inlet pressure considerab! 


the 


uring approximation 


mid i 
minor 
of 


and tend to minimize error ip 


one the most 
teresting pl uses of our report 
the author states that 


1 of 


t If we use the same figure and enter with a 


ferring to Fig. 0 (a) of the paper 


0.8, the error in 1/8 is about 20 per cent for an 


», Which is corre« 


value of S oi, t 1, the two curves give Vv alues of n of 


ibout 79 per d 81 per cent 1 difference of about 2 per 


cent 
On 


i plot fn versus S, the two theories in question 


give 
between them in the 


hand 


distance 


ll. Om the 


together 


two lines close 
the n- other 


the 


direction of ixis is sma intercepts 


in the direction of S-axis are farther apart. 


B. L. Newkirk.’ Fig theor 
taining to it, throws light on the question of stability of 
In the 


study of the oil-film whirl‘ it was found that an unloaded journal 


2 of the paper, with the 


in its equilil rium position course of an ¢ xperimente 


vould whirl at of one half its rotation fre 


1 whirl frequency 
, 

how 
be 


a frequency equal t 


this whirl 
it did so 


» one half the shaft rotation speed 


quency when running at speed At higher speed 


to oscillate 


lisappeared, but when the 


at 


airing was free 


and 


ASM 


(re 


* Consultant, 17 Ros 1, Schenec 
‘Shaft WI 
Re 


tady, N.Y Mer 
Oi Action in Journal Bearings 


HS 


ig Due 
28, August, 1925, pp. 559 
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Hagg® found a similar result using a relatively large bearing and a 
light shaft 

This motion does not appear generally at higher speeds or at 
heavier bearing loads. In one case® in which the unit load was 
very small (of the order of 1 psi), the shaft ran steadily up to the 
top speed of 30,000 rpm 


at low speeds 


It might have occurred, and escaped 
notice, In one commercial machine having a 
7-in. journal with a cap load, the half-frequency whirl was a 
disturbing factor 

That a half-running-frequency stimulus is present at all speeds 
may be inferred from the fact that rotors running on journal 
whirl at when 


bearings develop a their resonant frequency 


running at twice that frequency The internal damping of 
i minute stimulus of the natural 
It would 
that journals, in general, do not follow such paths as the larger 
Motions of this magnitude 


rotors is generally low, and 


frequency calls forth a response seem, however, 
orbits shown in Fig. 2 of the paper 
could hardly have escaped detection in the studies reported in 
references mentioned.** Hummel’ in his study of the theory of 
the half-bearing concluded that all journals running with ec- 
His experiments confirmed 


this 


centricities less than 0.7 are unstable 


1 mild instability restlessness) in region ; however, his 
clearances were abnormally large (for example, 0.006 in. per in 

and his bearing was substantially without an upper half. He 
ess than 0.004 in 


states that with “‘small”’ journal clearances 


per in smooth operation 1s possible. 
The effect of this phenomenon on behavior of the rotating- 
plug viscometer ought to be studied. It is possible that with very 
small clearance the plug runs essentially in a centered position 
With normal bearing clearance an unloaded journal certainly 
moves out of the center and moves around 

The writer cannot subscribe fully to the statement just above 
Equation [25a] of the “This is the origin of the oil film- 
excited shaft whirl.” 
ence in the present state of knowledge whether one will agree that 
a whirl at, say, 1000 whirls per min when the shaft runs at 5000 


rpm, is due to a tendency of the journal to whirl at 2500 whirls 


paper. 
It is perhaps a matter of personal prefer- 


per min, slowed down by end leakage or other causes, or to 


ascribe the whirl in such cases to a different instability charac- 
teristic of the oil film. We are now studying a rotor that runs 
without any severe whirl up to 5 times its critical speed and then 
breaks into a vigorous whirl with the whirl frequency equal nu- 


merically to the critical speed 


F. W 


as given by various analytical expressions both approximate and 


Ocvirk.* This paper compares bearing-load capacities 


exact From these comparisons it may be seen that in so far as 
the bearing of useful width is concerned, most of the analyses are 
approximate ones with the exception of the results obtained by 
Cameron and Wood Although the Sommerfeld solution is ex- 


act for the infinitely wide bearing, it too is approximate in the 
realm of real bearings 
Of the 


capacity, it 


equations giving an approximate solution for load 
the NACA 


216 


short bearings under constant 


appears to us of Cornell bearing proj- 


ect, that Equations [21a] and of the paper are the most 


realistic for load. Experimental 


data for l/d of 1 y and , as obtained by the projec t show good 

Journal Bearings on the Stability of 
Hagg, JoURNA! ApPLiep Me- 
1946, pp. A-211-220 

\ Nonwhirling Bearing t B. L. Newkirk 

ASME, vol. 56, 1934, pp. 607-615 

Nachgiebigkeit des Schmuer- 


* “The Influence of Oil-Filn 
Rotating Machines,’ by A. ¢ 
cHANICS, Trans. ASME, vol. 68 

§ “Oil-Film Whirl 
and Lloyd P. Grobel 

Kritische Drehzahlen als Folge der 
mittels im Lager by Charles Hummel, VDJ, Forschungsarbeiten auf 
dem Gebiete des Ingeni Berlin Heft 287, 1926 

* Associate Professor of Mechanical Engineering, Cornell Univer 


sity, Ithaca, N.Y 
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wesens Germany 
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igreement with these equations as well as with the curves of 
Cameron and Wood as might be expected from the comparison 
in Fig. 9 of the paper. 

Experimental oil-film pressure distributions for bearings of l/d 
of 1 as obtained by this project and by McKee and McKee are 
ilso in good agreement with the pressure-distribution function of 
Michell and Cardullo from which Equations [2la] and [21b 
ire derived 

The author's illustration of the overlapping of load-capacity 


t 18a} and [19a 


curves given by is extremely in- 
the load-capacity 


Equations 
that 
Sommerfeld for the infinitely wide bearing intersects the curves 
by the bearings. At 
intersections either of these equations 
The 


of short bearings at which the intersections occur is given by 


teresting. He shows curve given by 


given approximate solution for short 


these may be used to 


predict bearing-load capacity length-to-diameter ratios 


A similar consideration of Equations [19a] and [2la] of the 


paper gives the following 


As indicated by the author, the plotting of //d against n results 
in a curve which is an aid in determining the length-to-diameter 
range over W hich approximate solutions will predict conserva- 
tive load capacities. The area under the curve gives the region 
in which the short-bearing approximate solutions should be used; 
which the Sommerfeld 

The 
be looked upon as defining the region of bearings of intermediate 
width 0, is 1.73 for 
Equation [1] of this discussion and 2.50 for Equation [2] indi- 
cating that the short-bearing approximate solution given by Equa- 


the area outside of the curve is one in 
curve itself may 


equations give the more realistic results 


The maximum value of //d, given at n = 


tions [21a] and [215] of the paper is applicable toa slightly greater 
and [18) It is interest- 
Wood 


show 


range of //d than Equations [18a 
that the 
Sommerfeld 


note Cameron and curves 
that 
are applicable to a than 
Equations [1] that 


reach well into the region of intermediate-width bearings 


give 
they 
that 
they 


ing to 


intersections with the curve to 


much greater range of l/d 


given by and [2], herewith, and 


KE. M. Simons.’* 
the Op/dx term instead of the Op/dz term in the Reynolds equation 


As the author points out, the idea ot dropping 
is not new He is to be commended on his ingenious and com- 
prehensive application of this narrow-h« iring theory to the case 
Also, his 
unidirectional 


of dynamically loaded bearings comparisons with 


existing theories for steady loads demonstrate 


that this treatment provides results which are sufficiently ac- 
curate to be of practical use to bearing designe rs for be arings 
with l/d up to '/, 

“How do the 


results of this study compare with experimental observations on 


One question which naturally comes to mind is 


real bearings?”’ The writer has conducted laboratory studies" 


on a bearing with I/d = Apparently, there was sufficient 
damping in the machine used for these tests to prevent orbital 
motion of the journal center in all experiments with steady uni- 


directional loads This means that comparisons can be drawn 


* Research Engineer, Battelle Memorial I Columbus, Ohi 
Mem. ASMI 
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author's bibliogray 





DISC 


only with the poles of families of orbits given by the author's 
Equation [16] for this loading condition 

Fig. 1 of this discussion is analogous to the author's Fig. 8 (a), 
with the addition of an experimental curve for steady loads with 
an unbroken oil film. Unfortunately, the theories predict load 
capacities which are higher than those found by test. Further- 
more, the discrepancies are greatest at high values of 7, making 
it risky to use this work for design purposes when //d is as large 
as ' This same situation exists in the case of ruptured films, 
although, as may be seen from F ig. 2, herewith, the agreement is 


somewhat better, especially at the higher values of n 


Fic. 1 Comparison 
FORMANCE 


or 
BEARING 
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OPERATING 


Per 
Piun 
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hic. 2 Comparison Tne 


FORMANCE OF 


or 
Be aRine 


Prer- 


Finn 


EXPERIMENTAI 
Rueturep O11 


ORETICAL AND 


Oreratinc Wrrn 


The writer has found experimentally’ that when the film is 


ruptured under a steady load, the eccentricity ratio for a given S 
can vary over a sizable range, depending on the oil-feed pressure 
The limiting curve, shown in Fig 
lubrication 


2, was obtained with flood 
Presumably, however, if the ambient and oil-feed 
pressures were below atmospheric, a new limiting curve would 
result. Such observations lead to the conclusion that the ec- 
centricity in a real bearing, operating under a constant unidirec- 
tional load with a given Sommerfeld number, is not uniquely de- 
termined when the oil film is discontinuous. In real bearings, 
incipient film rupture occurs in a highly localized region of peak 
negative If the external load is or the 
lubricant-feed pressure decreased, the area over which the film is 


pressure increased 
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ruptured increases correspondingly both the axial and « 


Any 
films must take into account this concept of varving degrees of 


in ir 


cumferential directions realistic theory for discontinuous 
film rupture 

Since a given 8 can produce a range of 7 when the film is in- 
omplete it seems logical that there also would be a range of 
Fig. 3 of 
Although 
the author’s approximation agrees fairly well with the curve of 
Wood Fig. 3, explains the 


broad range in the experimental values of ¢ for agiven 8S 


possible locations of the line of centers for each S 
this discussion illustrates that this is indeed the case, 
Cameron and in neither theory 

When the test bearing was subjected to a constant rotating 
load \ N, 


‘ For speed ratios less than '/,, the orbits 


the journal-center orbits were circular when was 
greater than about 
were more complic ated Unfortunately, there was not sufficient 
lla 


and y (t) for comparison purposes with 


time to carry out the integration of the author’s Equations 
ind [114] to obtain 9 (t 
these latter orbits. When the orbits were circular, ¢ was not r 2 


COMPARISON 
POR BEARING 


THEORETICAL AND Opserveno Paase 
Operatine With Ruperurep O11 Fits 


OF ANGLE 





Frequency of 00d apo 


bic. 4 Errect or Recative Frequency or Le 
Recentricity Ratio ror Constant I 
CONSTANT SPEED 


APPLICATION ON 
7.8 Pst AT 


VAD 


OAD OF {OTATING 


as predicted by Equation (24) However, experimental curves 
Me N, had the shape predicted by Equation 24a 

For example, Fig. 4, herewith, shows that the equation predicts 
values of 7 which are about 20 per cent higher than the test 
results 


ot 7 versus 


In the tests, the predicted eccentricity ratio of unity when 
N./N 
tiny misalignments and extraneous loads were reduced or elimi 
nated However, such changes affected the curve only in the 
immediate vicinity of V,/N 


= '/; was approached when such nonideal conditions as 


= 3 


'! Reference (4), author's bibliography 
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AuTHOR’s CLOSURE 


The foregoing discussions have been particularly valuable for 


the reason that they all report on experimental observations on 
} ' 
considered in the 


theory 


bearings under conditions like those present 


theory and this, of course, is the proof of any From these 


discussions one would conclude that there is a fairly good agree- 


ment, particularly for the case of a constant load, and that the 


theory can be used reasonably safely for calculating the perform- 


ance of bearings under this type of loading having //d ratios as 
large as unity This should be useful in simplifying bearing- 
pertormanet calculations since the application of a side-leakage 
factor to the result calculated or read from charts for the infinitely 
wide bearing has always been required in the past Here, one 


single calculation using a relatively simple equation gives the 
final answer 

Howe ver, the agreement even in the constant-load case is not 
as good as one would like and definitely suggests that one or more 
of the factors omitted from the theory for purposes of mathemat- 


ical simplicity is having a significant effect in an actual case 
In particular, Messrs. du Bois, Oevirk, and Simons al! observe 
that the assumption of a continuous oil film gives definitely poorer 
agreement than the assumption of a ruptured film 
That journal-bearing oil films are often incomplete has been 


known for a long time. Recent experiments providing direct 
visual observation of the film have been conducted in the M.1.T 
Lubrication Laboratory which substantiate this point. In these 
experiments a glass bearing having //d = 1 was used and the 
oil film was re ndered visible by adding a fluorescent dy eto the oil 
and then viewing the bearing with ultraviolet light The rup- 
tured portion of the oil film was plainly visible and its behavior 
was found to very clearly substantiate Mr. Simons’ statement 
that the 


change the extent of the ruptured portion 


lubricant feed pressure, among other factors, should 
The oil film was found 
to be entirely continuous if the supply pressure was sufficiently 
high. As this pressure was decreased keeping the speed and load 
constant, a point was reached where a very small ruptured or 
cavitated region appeared, located on the center line of the 
bearing and in the diverging region approximately where the 
predict the minimum absolute 


ontinuous-film theory would 


pressure. With further decrease in feed pressure, this ruptured 

region increased in extent both axially and circumferentially.'* 
spect it disagrees with the present narrow bearing theory 

1 require 


s the ruptured region to extend from @ rto@ = 2x 
oads 1 
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Obviously then, as Mr. Simons concludes, the loc 


shaft will be different for different supply pressures, as 
each of these different film dispositions 

Mr. Simons is quite correct in concluding that “the eccentri 
in a real bearing, operating under a constant unidirectional load 


vith a given Sommerfeld number, is not uniquely determin 
when the oil film is discontinuous.” 


inalytically It is seen most easily in the 


vide bearing where the expression for 
resulting from the firsc integration of the 


be written 


dp r\2 cos 6 cos 6,, 
- 120 ( ) \ 
dé r ” ” 1 + 7 cos 


vhere @,, is the From the equation i 


seen that #,, and @_ ,, correspond to the angular positions where the 


integration constant 


pressure gradient is zero. Furthermore, as stated on p. 395, th 
orrect boundary conditions for an incomplet« 
the pressure and the pressure gradient must b 
edge Hence the position of the trailing edge is determined 
the value of 6,, 


0 Therefore the boundary conditions for determi: 


m 


specifically, it is located at the angular positior 


ing the tw 


integration constants in the expression for the pressure resulting 
| 


from integration of Equation [1] above may be expressed ana 


ically: 


lear!y seen that the limits of the oil film are no 


From this it is clear 


uniquely determined but must be expressed in terms of 


unknown constant 6,, 
the film to obtain the 
taining @,, so that the load is not uniquely de 


Further integration of the pressure ove 


total load yields an expression still cor 


termined either 


Mr. Simons concludes. Obviously, another boundary conditior 


is needed and it is fixed very simply by the requirement that 


the position of the oil inlet the pressure in the oil film must equa 
the oil-inlet pressure. Applying this condition, it becomes pos 
sible to express @,, in terms of this pressure The same argument 
applies in the general case where side | is considered 
though it is more complicated since the 


llas 0 


depend on Z as we 


18 Strictly speaking, p equals the o 


region but this is substantially zero unc 
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Vibration and Shock ‘ olation and shock force is considered both for the cases of excessive and 


for intermediate deflections of the isolator. The roles of damping 
VIBRATION AND SHock Iso.ation. By Charle de. John and of nonlinearity, and the support of excessive static loads 
Wiley & Sons, l New York, N. Y., 195 x 9 in., 


. a form the subject of the following chapter There is also a brief 
xvii and 328 pp., illus., bibliography, tables, index, $6 


section on sound isolation. Properties of materials, design of 
Reviewep sy G. Horvay isolators are treated in chapter 5. The last chapter discusses 


rir a variety of industrial equipment in regard to particular vibra 
HIS timely and up-to-date book on a subject which ac quired ite © i 
tion or shock considerations that may be appropriate 
prominence only in recent years, by a man who actively “ : 
- The reviewer feels that in a parallel development of vibration 
participated in its development, is a welcome addition to any 
: . and shock an emphasis on contrasts, rather than on similarities 
design engineer's bookshelf One finds it hard to realize today ‘ 
‘ _ os ~~ may more readily penetrate a tired engineer's mind; he would 
that only 30 years have passed since 8S. Timoshenko’s arrival ‘ 
ae ; have also preferred a somewhat more extensive treatment of 
in our country rhe inspiration exercised by him, first on a : . 
; , advanced topics, particularly sound isolation, at the expense of 
handful of engineers at Westinghouse, and later, through his - 
, the more elementary considerations in the early chapters Phe 
teaching and books, on nearly every student of mechanical en- 
’ almost complete lack of foreign references was also somewhat of 
gineering, set off a tremendous upsurge in applied mechanics. 
,  & surprise. Last, but not least, he wished the author had gone 
Vibration analysis became an important subject, both in school ; 
. into greater detail in the description of shock stand characteris- 
and in industry I'wo dates seem to me particularly outstanding; i } 
1 tics, practice in military Gesigns, and similiar subjects, the de- 
the early date, 1928, when J. Ormondroyd and J. P. Den Hartog, I & 
, velopment of which Mr. Crede was so intimately connected with 
two young associates of Timoshenko, wrote their important paper ; ; 
» Evidently, security considerations did not permit more details 
on the “vibration absorber,” an idea that exercised profound — ,, ' , ' 
' , These are, however, minor complaints Lhe book is well written 
influence on the development of the subject, and the recent date, 
. > , and bound to acquire widespread use because it is informative 
1945, when R. D. Mindlin has, in another important paper, 
“¢- authoritative, and replete with practical information It will 
established a permanent pattern for many of the approaches ly { a : ' 
’ greatly wilitate the design engineer s ish t also fills nin 
used today in the study of vibration and shock. Crede’s book 
1 portant gap in tl literature 
deals with a branch of vibration analysis, that of vibration and < 


shock isolation The book is evidently not intended as a text- ° 2 ve ° . 
book; the subject is too specialized for that purpose. Rather, Aircraft \ ibration and I: lutter 
it is intended as a guide or handbook to the practical engineer 
whose job is ) rly design or install equipment that is ex- 
posed to vibration or shoch 

In addressing himself to a practicing engineer, Crede faces a 
deli ate questior Does the engineer, in the field, know more Reviswep sy Pav A. Lansr* 
or does he know st n his more studious counterpart in class- 


room? I on the job, forgotien part of his fundamental puat problems of aeroelasticity have occupied 





mathematics, or has he progressed well bevond his student days? of iutical engineers for the past 30 years is attested by 


ok, thinks for a moment that the re- th I volume of technical literature treating sucl 

lation si + = sit " 2 + ay cause blem ring im However, there has never beer 
difficulty to tl eader, a provides a referer the sub- » which ib idered satisfactory r the tes 
ject might be reviews Fortunately this is the only lapse of this ] co g these problems or which presented 
kind in the book. The text is on the level of the average book on — suc« 1e st velopment in the aeroelastic field at ar 
vibrations, and to some extent also overlaps in content. But the given time he r nt book by Professor Scanlan and Mr 
subject of vibratior id shock isolation is explored thoroughly Rosenbaum admirably fills these needs; it will be found t ny 
in a manner not ile elsewhere in book form. Most of the professors a useful teaching aid and by many practical dynami 
important formulas are presented in convenient charts;a de- engineers a convenient source of much useful information 
tailed discussion is given of the mechanical properties of rubber, The first two chapters of this book present the mathematical 
felt, and cork techniques which ar ploited subsequently Thus one finds a 

After an introductory chapter on b notions and the cal discussion of the mn of ordinary linear differential equations 
culation of stiffnesses in nonprinci directions Crede treats with constant « i ts, the solution of some entary partial 
in chapter 2 the vibrations of the one-degree-of-freedom and differential equatic ements of determinant and matrix theory, 
multi-degree-of-freedom systems due to excitation by applied and the use of the Lagrangean « juations to est juations 
periodic forces or by prescribed oscillations of the support. The of motion for dynamic systems 
notions of force and displacement transmissibilities are intro- The third and four nay be consid d i short 
duced I have found the notation on page 44, using the same but clear treatment of the material usually presented in under- 
letter both for length (positive quantit, and for co-ordinate graduate vibration courses, that is, the free and forced vibration 
(positive or negative quantity) somewhat confusing The of one and two degree of freedom systems 
method of decoupling modes is explained in great detail and with The remaining twelve chapters of the book are concerned witl 


great clarity. In the next chapter isolation of shock motion aeroelastic or related problems and methods of their solution 


Engineerin ision, Knolls Atomic Power Laboratory, General 2 Associate Professor of Aeronautical Engineering, Pol 


Electric Comp henectady, New York Institute of Brooklyn, Brooklyn, N. Y. 
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Thus several practically important methods of obtaining the un- 
oupled natural modes of bending and torsion oscillation of 
beams are presented. The modes so obtained are then applied 
to the three-dimensional wing and tail-empennage flutter prob- 
The use of coupled modes of vibration, as may be obtained 
The 
whose vanishing determines 
The 
seroelastic proble mis of wing reversal and aileron reversal and re- 
The book 


closes with a discussion of instrumentation and testing and a pre- 


ems 


by vibration test, is then discussed various methods of 


solving the stability determinant 


the critical flutter velocity and frequency are presented 


duction in aileron rolling effectiveness are discussed 


iminary treatment of the flutter of sweep wings. The treatment 


of the aerodynamic forees for the flutter proble mis presented in 
letail in an appendix; only the coefficient representation is, there- 
lore, given in the main body of the text 
In their preface the authors point out that the topies and pres- 
This 


that people using 


entation reflect their own experience and point of view 
is quite natural; it must be expected, however, 


this text will, for the same reason, prefer in certain sections their 


own points of view. Thus, for example, the reviewer feels that 
he presentation of the various methods of modal analysis of dy- 
imc prol lems, both steady state and transient, as an expansion 
n terms of various complete orthonormal] sets, that is, the un- 


oupled or coupled, free body or cantilever modes of vibration, 


brings considerable clarity and unity to the various techniques 


which are presented in this text and in the technical literature. 
Che possibility of other points of view does not detract from the 
hook; the 


ing avoids i the textbook literature 


suthors must be congratulated in so successfully fill 


Differential Geometry 


Dir METRY By Dirk J 


Matt 


( ASSICAI ERENTIA (; 


Addison-Wesle 
Mass., 1951 


ities Series 


und 221 pp 


Reviewep By H, L. LANGHAAR 


textbook on the classical differential 


Hes is an authoritative 
geometry of curves and surfaces in three-dimensional Luclid- 
:mphasis is placed on the historical development 
is it stemmed from the works of Monge and Gauss 

on to the 


In addit 


uded a surprising amount of supplementary material; the general 


standard basic theory, the author has in- 


theory of order of contact, theory of integration of natural equa- 


urves, theory of curves of constant slope, evolutes and 


thons of 
involutes, imaginary curves, a few theorems of differential geom- 


etry in the large, triply orthogonal systems of surfaces, theory 


of curviline co-ordinates in space, geodesic co-ordinates, varia- 


tional theerems on geodesics, surfaces of constant Gaussian 


irvature, the Gauss-Bonnet theorem, interpretation of the non- 


huclide 
f 


ve lopes ot 


in geometry of Lobachevsky and Bolyai, theory of en- 


families of surfaces, conformal mapping of surfaces 


sometri nd geodesic mapping of surfaces, theory of minimal 


surfaces, theories of ruled surfaces and de 


elopable surfaces, 
] il 


ind surface theory 


naginaries in 
The author 


reader 


s stvle is direct, but concise Needless to sav, the 


requires an understanding of the fundamentals of mathe- 


matical analvsis Some knowledge of Gibbs’s vector analysis 


* also advisable, since the brief explanation of this theory in 


the book is scarcely adequate. Gibbs's vector notation is used 
throughout, but tensor notation is not emploved 


The 


ones from 


text contains many drawings, including some notable 


ther books. Also, there are some illustrative photo- 
Theoretical and 


versity of Dlinois, Urbana, Ill 


’ Professor Applied 


Mem 


Mechanics De partment, | 
ASMI 


APPLIED 
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graphs of plaster models. Numerous examples are worked out, 


ind many problems for students are given, Biographical sketches 
of the lives of the principal contributors to classical differential 
novel feature of the book. In all, 


teachable graduate 


geometry are an attractive 


the requirements for an instructive and 


fulfilled 


mathematics text are 


Algebra of Vectors and Matrices 


Wade 
ambridge 


x an Y py llus., bibliog- 


Matrices. B 
Mather 


Thomas L, 


Serie ( 


VecTORS AND 
Press, In 


loth. 5 x SS 


THe ALGEBRA OF 
Addison-Wesley 
Mass., 1951. ¢ 
aphy $4.50 


atics 


Reviewep By R. F. Tinenarr* 


"THERE has lack, cited by the 
this book’s preface, of a good textbook for a first course in 
There This 


neither the mathematiciar 


existed a serious withor in 


the theory and application of matrices still does 


book seems likely to find favor with 
nor the physic al scientist or engineer 
The stated aim of the book is to present in elementary ex 
position of the algebra of vectors and matrices, that exposition 
being articulated with the basic concepts of modern algebra 
The author achieves his purpose to the extent that the exposition 
is lucid in most of the places where his mathematics is correct 
The book opens with a brief exposition of the concepts of group, 
ring, integral domain and field, followed by an elementary treat- 
ment of the algebra and geometry of two and three-dimensional 


A chapter on the algel 


vectors algebra of n-dimensional vectors 


setting for a itment of matrices trom 


After the un 
established, 


( ompletes the stage 


damental properties of 


the vector standpoint 


matrices have been linear transformation theory, 


linear 


The 


the characteristic equation, rank equations, 


ind bilinear and quadratic forms are t1 exercises are 


not very numerous or meaty 


In contrast to mathematical uracy, the selection of 


material and the order of its presentation subject to or 


remembers that the is designed for 


only No 


The reviewer believes that the “law 


minor criticism, if one 


a beginning course nongeomet applications are 


given of nullity” is impor 


presentation 


tant and elementary enough to have been included pri 


pal defect with regard to order of material is the 
of the characteristic equation before s 
This robs the 
the 


ind barren 


linear equations 


been discussed equation of 


makes 


somewhat awkward 


have 


motivation and discussion of i teris vectors 


A number of theorems in the bo« without proof 


with references to proofs provided. Proofs of some of these 


could have readily been included without violating the elementar 
spirit of the book and without perceptible extension of the text 

On hand, 
might better have been stated without pr 


Proofs 


the other some of the theorems hich are proved 


‘ool 


view 


containing maj« theoren 


therein 
on the equality of the row rank erminant 


140 


regarding a fundamental set of 


rank of a in the implies the 
154 
homogeneous linear equations; in the theorem transformation 


of a 165 
Other confusions are 


matrix (page orem (page 


solutions of a system of non 


svmmetric bilinear form (page 


ind 


are of such a nature as to perplex ar 


mathematical errors plentiful and 


many wut the best students 


For example 
two theorems ‘ 
is called 1 field F, when 


are any elements whatever of S, then p + 9, 


these 


1 Page 10—‘As a 


set of two or more elements 


consequence ot 
numbers 
p and q, (q #90 


* Profe 


Ohio 


r of Mathematics, Case Institut rechnology, Cleve 


land 





Pq, P-9 iP 
tributive laws are appare ntly needed 

2 An effort is made 
matrices by a 


q are also elements of S No 


assoc1atiy 


page 56) to 


points.” 


motivate the sum 
sum ol 
confusing 

3 After defining a matrix, addition of matrices, multip 
of a matrix by a scalar, and proving the theorems (six 


B+A,(A 4 +C =A4+(B+C),rA + 
rs)A = B) =rA + /rB, (rand s scalars), the 
then decides it is time to discuss equality ot matrices 
4 Having defined a scal 
summarizes the 


Theorem VII 


sA 
sA),r(A 


' 


ir as an element of a field, 
last four theorems mentioned in the st 
Linear combinations of matrices 


coefficients obey the laws of ordinary scalar algebra.” 
5 On page 70 the author 


which all 


wwona 


with 


defines a diagonal matrix 


the elements 


T his 


square matrix 


the 


in are zero except 
long principal di ads to 


contusion 
the 
6 


diagonal 
On 


Invers¢ 


page 84 the author alludes to a right inverse as 
of a matrix 


verse is also a left inverse and that it is unique 
7 On page 168 


field of characteristic zero, the author then restates the 
for ther 


proof he regards as equired in this case 


In several instan 


not stated (« 2 


exer 


£.. 
stances the neces 
“the 
through the book 


ise page 96, page 164 


is not stated 
sufficient 


In other 
sary ‘“‘squareness 
The phrase necessary and condition 


The reviewer | 
graduate 


s just used the 


havi 


book as a text in a course 
engineers 


f belated publi 


m the 


book did not 
healtl 


to read the text wit! 


basis of its table of 


work badly The 


conten 
Strangely enough, tl 


the 


errors 
book have the i 
gre are All in all, 


it was good fun 


Applied Mechanics 
1pD MECHANICS 


App 


bridge Universit 
1951. Clot 


N 
Lond 


q in 


RS By Sir Charles Inglis 
England, ar 


xu and 404 py 


Pre on d New York, N. Y 


437 figs., $7.50 


Reviewep sy bk. STerNBer 


"T’HIS book is based on lectures given by the 
working toward an honor’s degre« 


rhe level the 


ourse classification traditional in 


In engineering at Car 
bridge Universit of book is intermediat 


liate 
its scope defies the this country 
The first two chapters are devoted primarily to 


ind graphical treatment of plane rigid-body statics 
nterest 


in analytics] 


The author's 
in structures is mace 


apparent quite early 
o bridge design. There 


itenary and two chapters on the 
itions in ideal trusses, 


Vv excursions 
Into questions re tei t 


ible and 


and deform 


follow a discus- 
sion of ¢ problems, 


stresses 
cally CASE 


the treatment of statics 


including some stati- 


indeterminate 4 chapter on dry friction concludes 


The remainder of the book deals with dynamics, and includes 
several chapters on the kinematics and kinetics of particle motion 
is well as a discussion of plane rigid-body motion 
vibrations are given considerable 
the subject matter of five chapters 


Mechanica 
constitute 
f one and 


prominence and 
Systems ¢ two 

Professor of Mechanics, Illinois Institute of Technology, Chicago 
Ill. Mem. ASM} ° 


This seems meaningle ss and 


lication 


{ 
he : 


when it becomes clear that he intends to permit zero elements on 


then subsequently shows that a right in- 


having established the theorem on the re- 
ducibility of a quadratic form to a diagonal form for an arbitrary 


theorem 
eal field after supplying the three additional sentences of 


es the necessary nonsingularity of a matrix is 


e.g., page 43 
abounds 


for 
g chosen the book sight unseen, because 


effect of stimulating the good student 
Cam 


suthor to students 
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e or dis degrees of freedom, as well as simple examples of one-dimensional 
continuous systems, are covered A separate chapter on non- 
ot two 


WAVES, 
noteworthy 
discussion of gy roscopic motions, 


linear restoring forces, pressure and 
lations is particularly 


The 


self-excited oscil- 


book with a 


It is clear from the beginning that the 


+ s)A 


wuthor 


withor is anxious to 
ind specific applications, 
This desire is probably responsible 
hurried treatment accorded to basic 
reader 


reach techniques of solving problems 
with a minimum of delay 
the occasionally 


Thus this 


lor 
ples 
tha 


princi- 
wonders whether the 


two-dimensional 


statement on 
withor 


number of 


page 


force 


in a system a 
in always be balanced by a single force 
interpreted. The however, 
in accomplishing his primary 
with genuine applic 
1 host of illustrative 
are, 


irtling =< 


s 

might not | 
sealar succeeded 
The 
in elementary 
es which ar 


«* TDIS- 
author, has idmirably 

objective hook 
as “a 


level 
e presented with 


and an uncommon attention to numerical detail 
mechanics 


ibounds 
tions even at it contams 
those vam unusual 
later ( 


Teachers 


exercises 


of ire apt to find the large supply of 
helpful 


rhe 


most 
the emphasis on approximate methods is a novel feature 
Numerical and graphical integration of differential equations 
us well as harmonic analysis, are used, at a relatively early stage, 
to bring withir f the book stimulating problems of 
is refreshing to note the conspicuous 

ibsence of the myriads of unhappy pulleys which have long popu 
lated the pages of traditional texts on mechanics 
The point of vie 


resolution 


the s« ope 
engineering significance It 


w adopted by the 
employ the simplest 


preference to those which may | 


author is 
to mathematical methods 
in ” more ele gan bui less obvious 
This attitude has led the author to place a great deal of emphasis 
on intuitive insight into the mathematical methods used and upot 
physi al interpretation of the 


results obtained 
hand, it has induced him t 


On the 
methods 
the “mysticism 
to allude to the 
with his treatment 


other 


to dis 


clear of ve« 


tor 
miss Lagrange’s equations with a refers 


nee to 
ts 


of an ingenious mathematical technique 
“red herring of convergency 
of Fourier More 
the » make 


to become a master 


ind 
in 


* in connection 


he 
it servant 
narks of this kind 
hey are also apt 

For, as the 
wa ol 


| escaping mastery by 


series general] announces hi 


and 
likely 


to cre 


s ideal in 
pretace, “‘t« mathematics 


> 
te 


hows ver 


not allow it 


are to appea 
to many readers 


ite 


an impres 
sion which was not intended 


suthor has often demor 
strated, one 


mathematics is to 
to master it 


This is not just another text on mechanics; it is a ve per 
sonal 


nd eloquent account of the subject matter which through 
wit bears the mark of its author's long and rich experiet 
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exercise 
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Reviewep sy C 
‘THs book forms a part of the 


Series under the editorship of I 


elastic solids where the 


W. MacGrecor 


Wiley 


s 


Apple <l 
Sokolnikoff 
le that 


such that the 
products of the strain components are negligible 


Mathematics 


and treat 


assumption is not 


sufficiently small 


raat leformation 
magnitudes are 


squares and 
It is a subject 


on which the author has published frequently and 


spread acceptance 


wile- 


with 
* Vice-President in cha of 


engine 
University of Pennsylvania, I 
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Loffective use of the theory of matrices has been made through- 


out the beok. As in the modern treatment of many other engi 
neering problems, matrix theory is so useful here as to be beyond 


An «¢ 


ments in the text 


xamination of the general proofs and develop- 


this 


ontention 
For those wishing a re- 
the first 


will substantiate 


matrix theory or even confronting it for time, 


chapter l provides idequate treatment for the material which 


ollows 


Chapters 2 and 3 are devoted to the geometries of strain and 


stress and the relationships between them. In the treatment of 


isotropic elastic media in chapter 4, the influence of initial states 


isotro are considered The increased number of 


of stress on 


elastic constants for isotropic bodies under higher-order strain 
effects are also outlined, ompressibility relations are developed 


between hydrostatic pressure and volume, both for the linear and 
made of 


strain theories Comparisons are 


the finite high-elastic 
these with Bridgman’s experimental results 


The 


dis- 


Nonisotropic elastic media are treated briefly in chapter 5 


elastic sensitiy of the medium to special rotations is 


cussed 


Ihe separate problems of simple shear and tension are de- 


veloped in chapter 6 where both the elementary linear theory and 
The last chapter 


the effect of in- 


ond-order corrections are determined 


the sec 


iins solutions of three problems, name 


rnal and external pressure On & spheric al shell, on a circular 


‘ ind the torsion of a circular vilinder In each case the 


itions based on elementary considerations are compared with 
ose considering finite elastic strains 
lhe reviewer found this book to be a stimulating treatment of a 
lected While clearly written 


added greatly 


the use of more 
to the 


0 long neg 
lew hgures would have 
and other varia- 


a 
entation rhe notations for stress, strain, 


do not follow current practice ountry, and especially 


this reason the inclusion of : f nomenclature would 
tive The latter ar I ! or crit 


ciable 


icisms and 


strued ul ippre manner 
er hopes that the author 


inc that 


The 


ve VOrkK 


of this book 


further 


Linear Camputations 


‘8 By 1 wyer John Wiley 
Y., 1951 loth, 6 9 in iand 


& Sons, 
$44 pp., ap 
lexe 5 


Reviewep py Georce bk. Forsyrut 


putations” the the solution 


iuthor means 


quations, the inversion of a 


' 
erminant il 


inear algebrai 


ition of a de 1 similar calculations 
ion of refined Gaussian 


putations on desk 


devoted to the « 


forming these ir com 


ire all elimination method 


but few omputers are aware of the 
ide possible by over a century of development 


rked out Doolittle, Cholesky, Banachiewicz, 


take full advantage of the ability of 


umulate products easily, and reduce the slow 


tminimun 


is a profe ssor of statistics at the 


University of 


In the problems of experimental statistics he is con- 


it iterative methods are less successful than pivotal 


methods in earrving out linear computations 


ind Southwell’s school that the conclusion 


wo of Hardy Cross 


Los Angeles 


APPLIED 


I gather from the 
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does not carry over to all engineering calculations. Even so 


elimination methods are important, and the present book must 


surely be a standard reference on these techniques. There is 


nothing comparable in the literature 
Four or five chapters are devoted to calculating machines and 


to the principles of designing a computation for them. One 


chapter deals numbers as 


with computation with approximate 


This material is drawn upon in d 


etermining the errors in 


such 
linear computations. A number of chapters deal with such 
problems as obtaining the characteristic polynomial of a matrix 
evaluating a linear form whose variables are solutions of a linea 
system, and getting an inverse matrix. 
Professor Dwyer’s book is written at an unusually elementary 
level Nothing but easy high-school algebra is needed to follow 
the de velopment of the efficient elimination schemes for solving 
Then 
and their elementary properties are described and used 
thereafter The therefore be a 


computers with littl 


linear equations determinants and matrices are intro 


duced, 
book will useful reference to 
The elemet 


thar 


The reader who desires to develop 


mathematical background 


tary level is partly a result of presenting techniques more 


tools of numerical analysis 


analytical tools would do better to start with a good grounding 


in elementary matrix theory, including the triangular decomposi 


tion and the use of various norms 


The numerous carefully prepared tables include layouts and 
! 


one or more numerical examples of almost all the methods dis- 


The 


| nglish langu ge 


cussed lists of references are unusually complete im the 
rather modest on European work, and quite 
devoid of Russian entries. The book is well printed, with few 
misprints or other errors 


The chapter titles are l 2) Computa 


Introductory remarks 
tion with approximate numbers, (3) The principles of computa 
The solution of simultaneous equations with 
the method of multiplication and subtraction, (5) The method of 
t) division—method of 
with 


methods, (7) Relations between the coefficients, (8) The solutior 


tional design, (4 


multiplication and subtraction with (exa 


determinants, (6) The solution of equations approxima 


of related and associated equations, (9) Introduction to deter 
10) The evaluation of determinants and determinanta 


12) An intro 


and its cal 


minants, 
The evaluation of linear forms, 

(13) The inverse 
(14) The caleulation of the 


thed of 


ratios, (11 


to the algebra of matrices, matrix 


culation with methods, 
and 
involving the ch 
17) The 


statistics 


ipproximate 

determinants, (15 
Othe 

Applic 


adjoint inverse with the n 


Problems iracteristic equation, (16 


methods, errors of linear computations, (18 
tion to 19) Application to nonlinear problen 


cluding remarks 


Introduction to Acoustics 


An IntTRopU« Lobe H. Ra 
Wesley Press s., 1951 
tables, ref me ind prot s, xii and 


Reviewep By Samvuet Levy’ 


t students sics and engineering who plan to special 


, agp + 
in communications, this book will form a valuable basi 


understandable form, the phvsical 


ind the 


It presents, in easily 


principles underlying acoustic phenomena methods of 


ipplying these principles to engineering problems. Thus one 


finds both the development of the wave equation for the propa- 
ippli- 


gation of sound in a gas and the engineering equations 
4 1 


transmission of sound in 





BOOK 


The author the 


of other 


has emphasized, wherever possible inalogy 


between the equations of acoustics and those branches 


of physics. His book is thus of added value in giving the student 
a broad view of physic 8 I xamples of his ise Of analogies are 
found in his comparison of acoustic with electromagnetic waves 
Th 
mpedance is 
ot 


and of the electrical transmission line with an air conduit 
use of complex notation in discussing acoustic 
presented with frequent comparison to the analogous use 
complex notation for impedance of electrical circuits 

The mechanics of vibrating strings, plates, and bars is pre- 
sented in detail This is followed by a discussion of musical 
The 


to point out to the student the great difficulty of setting down 


instruments withor takes advantage of this discussion 
the boundary conditions of most actual acoustical problems in 
precise mathematical form and the need, in many cases, of mak- 
He feels 
that the student must be impressed at an early stage in his learn- 
\ 


review is given of the available knowledge regarding speech and 


ing approximations to obtain an answer at all wisely 


ing with the practical importance of approximate methods 
hearing both from the subjective and mechanical points of view 


und the 
tages and disadvantages of the different types are 


The methods of designing loudspeakers are given idvan- 
compared 
Numerous other examples of the applic ation of acoustic theory 
are presented 

The large number of problems included in the text for solution 
by the student have been carefully selected to emphasize the 
The 
prepared and add to an understanding of the subject matter 

This book be to the of 
undergraduate students in physics, and communications engineer- 


ing 


principles presented many illustrations have been well 


will welcome teacher acoustics, to 
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REVIEWED BY SAMUEL Levy® 


"THIS is a book on advanced acoustics written primarily for 
graduate students in either physics or engineering electronics 


It should 
also, however, find wide use as a reference for practicing engineers 


The book is intended as a textbook for classroom use 


because of the broad coverage, the detailed consideration of the 
basic assumptions, and the derivation of each important equa- 

The 
a“ good check 


tion from the fundamental laws of physics comprehensive 


set of excellent problems should afford on the 
student's grasp of the subject 

The first nine chapters, comprising 264 pages, present clearly 
the 
solid bodies 
and of fluid 
The presentation of this material is illuminated by the mention 
The 


and concisely the fundamentals of theoretical acoustics in 


form of analyses of the various types of vibration of 


the propagation of sound waves through media 


of practical applications for the equations that are derived 


use of vectors to express harmonic variation with time 4s 
commonly done in electrical engineering, is presented logically 
The definitions 


of motional impedance and transfer admittance proceed from a 


in discussing the motion of a simple oscillator 


discussion of the dynamic loudspeaker The concept of wave 
impedance is presented in connection with the motion of strings 
The presentation of vibration theory includes dis« ussions of the 
1xial and transverse vibrations of bars, the transverse vibration of 
The treat- 


membranes, and the transverse vibrations of plates 
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ment the derivation of the 
wave equation and proceeds logically with discussion o 
of 


reflection at interfaces, and absorptior 


of transmission starts with 


sound 


gener: 


spherical acoustic waves, directivity waves, radiation imped 


ince, resonators, filters 


The 


cerned with applications 


con 
cle 


loud- 


remaining seven chapters, comprising 235 pages, are 


The topics considered include the 
baffles 


A chapter on p 


sign of dynamic loudspeakers, cabinets horn-t ype 
ill kinds cho 


speakers, and microphones of 8 
acoustics describes the behavior of the ear, the subjective cl 
teristics of loudness, pitch and timbre, and binaural localizatior 
\ chapter on architectural acoustics provides design information 
nd the 
covered in unusual detail with the pres 


on reverberation time damping of standing waves 
Underwater acoustics is 
ol 


and temperature 


entation much valuable information regarding attenuatior 


refraction gradient effects in sea water well 


as 
ind discussions of echo ranging 
The last chapte r, dealing 


as design information for sonar 


and underwater communication vith 


ultrasonics, discusses the generation of ultrasonic waves, the ce 
tection of flaws and cracks in solids, and propagation phenor 

in gases, liquids, and solids 

this book 
vetween theoretical and applied acoustics 
he forefront 
should il 


d engineers 


In the opinion of the reviewer is well balanced as to 


content t Its clear and 


rigorous style puts it in t smong graduate textbooks 


on acoustics This text ao be of considerable value as a 


reference work for designers at 
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6 =x xii and 542 pp., 


LEVIEWED By Eric Retssner"® 


present book is concerned with the theory 
of Bessel 
esse] integrals, Mellin integrals, and Laplac« 


and appli 


trigonometric series, function series 
Fourier integrals, 
suthor 


transforms. In three chapters of some ninety pages the 


discusses the theory of the various procedures. The remainder 
of the book is devoted to various applications which in general 
are preceded by a brief introduction to the fields from which these 
applications are taken. We find problems of vibrating string» 
membranes and pl ites proble ms from the theory of heat conduc 
dime nshons 
such 
und 


tion, from hydrodynamics, from elasticity in two 


and with 
timely subjects as the slowing down of neutrons in matter 


axial symmetry Separate chapters deal wit! 


with applications to atomic and nuclear physic 8 
Some sixty pages of the book are devoted to problems of hydro 
the 


dynamics Incompressible potential flow is illustrated by 


problems of flow through an aperture (without free boundaries 
and linearized problems of gravity waves This is followed by 
the solution of some linearized problems in incompressible vis 


cous flow ame 


the 


There is one nonlinear problem, concerning pl 


viscous flow, where the Fourier tec hnique can be used, since 
differential equation 1s Of sec ond degree, to reduce the problem 
to a nonlinear integrodifferential equation for the Fourier trans 
form of the stream function. In the main the author has con 
tented himself with discussing classical results of the type found 
One 


applications of the Fourier te« hnique to the 


in Lamb’s Treatise regrets the omission of present-day 


fields of steady 
to the field 
in this connection, of the subject of generalized 


and nonsteady linearized flow and of 


turbulence and, 


supersonic 


harmonic analysis 
In elasticity the treatment is primarily concerned with bound 
the 


ary value problems for the half space, the infinite plate, and 
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infinite circular cylinder \ departure from the usual treatments 
in this field concerns the solution of mixed boundary-value prob- 
lems, such as the problem of the rigid punch on the surface of a 
half space, by Fourier methods. Such problems can be reduced 
by the use of the Fourier or Bessel integral to so-called dual in- 
tegral The author shows examples where these 
dual integral equations can be solved in closed form. 


equations 
However, 
the problems in question have previously been solved by other 
more direct techniques and where these more direct techniques 
break down the dual integral-equation technique likewise does 
not seem to supply a solution 

Perhaps not enough stress is laid by the author on the problem 
of evaluating a Fourier transform solution once it has been ob- 
tained In this connection one thinks of such problems as the 
forced oscillations of an elastic half space as treated by Lamb in 
a classical paper or of the problem of transient beam vibrations 
with transverse shear deformation taken into account which re- 
cently was found to be of considerable technical importance. 
book is 
the 


Aside from limitations such as these the very well 


written and represents a useful addition to literature on 


applied mathematics 


Bou nda rv-Layer Th eoryv 


Hermann Schlichting. Verlag und 
Druck G. Braun, Karlsruhe, Baden, Germany, 1951 Cloth, 
6/5 x 9 in., xv and 483 pp., figs., tables, bibliography, author 
and subject indexes, $12.50 U.S.A 


GRENZSCHICHT-THEORIE By 


Reviewep By Ascuer H,. SHartro" 


ij {! RE is a book which is indispensable to advanced designers, 
research workers, and students dealing with problems of 
fluid motion in the fields of aeronautical engineering, mechanical 
engineering, chemical engineering, and applied physics 
The appreciation of the value of the boundary-layer concept 
und of the place of boundary-layer theory in fluid dynamics has 
grown during the past 20 years from a point where comparatively 
few understood its implications to the present condition where it 
Yet the 
principal treatises (Prandtl’s article in volume 3 of Durand’s 
lerodynamic Theory; 
Fluid Dynemics) prior to Dr. Sehlichting’s book are now 15 years 
old. The present book, besides incorporating recent develop- 
ments in the theory of 
of transition from laminar to turbulent flow, and in the turbu- 
lent the superior to 


Prandtl’s and Goldstein's treatises in respect to choice, presenta- 


is matter-of-factly accepted as an indispensable tool. 


and Goldstein’s Modern Developments in 


compressible boundary layers, in the theory 


boundary laver, is in reviewer's opinion 


tion, and organization of subject matter 

In his foreword, Dr. Schlichting points out that the book is 
written primarily for engineers, This aim, together with the 
of the book, seems to have guided the author in his choice 


Those theoretical aspects of boundary-layer theory 


tithe 
of material 
which are relevant to and useful in engineering design form the 
main subject matter of the book. Experimental data are intro- 
of validity of 


the 


to illustrate the accuracy or range 


not all 


duced mainly 


the theoretical considerations. Of course topics 


treated are of direct applicability to the quantitative solution 


of engineering problems—many give rather a general under- 


standing of the mode of boundary-layer behavior in typical 
situations and provide some basis for qualitative decisions on the 
In accordance with the practical aims 


part of the designer 
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of the book, there are numerous charts and tables which help 
to eliminate lengthy computations 

The book is divided into four main parts. In Part 
forth the fundamental concepts and laws for dealing with the 
fluids: The concepts of viscosity, dynami 
similitude, and of the boundary 
Navier-Stokes equations and the assumptions implied therein 
the general behavior of viscous fluids as indicated by the Navier 
solutions of the Navier-Stokes 
Part B covers the laminar bound 


A are set 


motion of viscous 


layer; the formulation of the 


Stokes equation; some exact 


equation; and creeping flow 
ary layer: The development of the boundary-layer equations 


the general behavior ol laminar exact 


boundary layers; some 


unsteady flows; integral 
the thermal 

Part ¢ 
flow, and 


solutions for laminar boundary layers 


methods of solution; boundary-layer control; 


boundary layer, and the compressible boundary layer 
deals with the transition from laminar to turbulent 
covers mainly the stability of laminar boundary layers to disturb 
inces of very small amplitude. Part D is concerned with turbu 
The apparent 
the semiempirical mixing-length theory 


lent boundary layers: turbulent stresses owing 
to momentum exchange 
of Prandtl, vorticity-transport theory of Taylor, and similarity 
theory of Karman ; pipes; turbulent 


boundary layers on flat plates and rotating disks; effects of pres 


von turbulent flow in 


sure gradients on turbulent boundary layers; problems of fre: 


turbulence; and some miscellaneous information relating to 
profile drag. 

The format, typography, and illustrations live up to the mag 
nificent prewar standards of German technical books, and set 
an example for American publishers 

It is to be hoped that some American or British publisher will 


book translated and 


take the initiative of having this valuable 
published in an English text 


Theory of Control 


THeory or Contro. in MecHANICA 
Macmillan Cambridge Universit 
New York, N. Y 1951 Clotl 


illus., examples, tables, bibliograp! 


AN INTRODUCTION TO THE 
EnGineertnc. By R. H 
Press, London, England, and 
7 X 10'/:in., xiii and 195 pp 
appendixes, author and subject indexes, $6 


Reviewep BY Joun A. Hrones™ 


"THIS book of less than two-hundred pages is one of the first 

devoted to the theory of automatic control and written from 
the viewpoint of a mechanical engineer. The chapter titles 
indicate the nature of the ground covered by the author, namely 
the principles of control, control-system elements, servo-system 
operation, servo-system performance, general servo analysis, 
unalytical methods, and graphical methods 

Also included brief 
numbers and the Laplace transform method 


The book is well written; is adequately illustrated with clear 


are appendixes dealing with complex 


sketches and diagrams; and is supplemented with numerous 


problems. It should prove to be more easily read by novices 


in the field than most existing texts 

The transient and steady-state behavior of simple single loop 
systems in response to square wave and sinusoidal disturbances 
and analytical methods for deter- 


is discussed. Graphical 


mining such responses are presented. Certain simple nonlinear 
systems are briefly The 
scope, will be of considerable interest to those who wish a brief 


covered book, although limited in 


ut competent introduction to the subject 
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